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Shipping of Radio-[sotopes* 


KARL Z. MORGAN 
Oak Ridge National Laboratory, Oak Ridge, Tennessee 


(Received December 19, 1947) 


A. GENERAL DISCUSSION OF THE PROBLEM 


HE need for suitably revised regulations 

governing the shipment of radioactive ma- 
terials has been recognized for several years. 
Present regulations prohibit the shipment of 
radioactive materials by mail! and specify the 
conditions under which 100 me or less of radium 
and radon can be shipped by railway express 
and by air.’ A few years before the past war the 


*Dr. K. Z. Morgan, Director of Health Physics at 
Oak Ridge National Laboratory, gave this paper before the 
Conference on Radiobiology and Medicine held at Brook- 
haven National Laboratory on October 18, 1947. 

‘One gram of luminous paint made from pure polonium, 
free of other radioactive members, may be shipped by 
mail if properly packaged. Dials painted with pure po- 
lonium are mailable. Otherwise, the following statement 
of the Postmaster General describes the case. ‘‘Any radio- 
active material, including radium, thorium, or other radio- 
active substances, or any materials containing radioactive 
substances, such as powders containing radium or thorium, 
liquids containing radium emanation, radium salts, or 
radioactive minerals, is prohibited transmission in the 
mails.”’ 

* The regulations at present (referring to the old regula- 
tions that were superseded January 26, 1948) prohibit the 
shipment of more than 100 mg of radium unless special 
arrangements are made with the local shipping agent. A 
table is available (L. F. Curtiss, J. App. Phys. 12, 346 
(1941)) which gives the lead shielding thickness required 
for up to 600 mg of radium and for different times of 
transit. The radioactive shipment must be placed 15 ft. 
or more from x-ray film and 10 ft. or more from ordinary 
film. The radioactive material must be securely packaged 
with labels indicating the radium or ‘“‘equivalent radium”’ 
content. According to Mr. R. B. Wilsey these regulations 
were intended to limit on the average the exposure to 11 
mr to x-ray film and 25 mr to other types of film. Supple- 


shipment of radioactive materials from Van de 
Graff machines and cyclotrons began, and the 
magnitude and number of such shipments has 
increased rapidly since because of the develop- 
ment of other high voltage accelerators and the 
widespread expanding use of radio-isotopes in 
research, medicine, and industry. It was not 
until the large scale production of radio-isotopes 
by uranium piles began and until the relaxing 
of security regulations since the war that it 
became possible to furnish fission products and 
pile-induced radio-isotopes on a large scale so 
that serious efforts were made to formulate a 
completely new set of shipping regulations. It 
was recognized that suitable regulations should 
be developed for radioactive shipménts by air 
as well as by rail because many of the radioactive 
materials are short-lived and often the most 
practical and economical method of shipment is 
the fastest method available. Also, regulations 
dealing specifically with radioactive materials 
other than radium, polonium, radon, and thorium 
were badly needed. All the prewar regulations 
were developed with particular emphasis on the 
protection of photographic and x-ray films from 
radiation exposure during shipment. Now the 
large increase in frequency of radio-isotope ship- 
ment 5 to the Railway Express Agency, Air Express 


Division, Tariff No. 8 (May 31, 1946) lists regulations for 
shipment of radio-isotopes by air. 
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ments and the fact that (potentially, at least) 
single sources can be made up rather cheaply for 
shipment, which produce a gamma-radiation 
equivalent to thousands of curies of radium, 
make it urgent that regulations be provided 
which will protect man as well as films. 

During this interim period the Health Physics 
Divisions of the Plutonium Projects (at the 
University of Chicago, Hanford, Washington, 
and Oak Ridge, Tennessee) have set up their own 
more rigid and complete protective measures 
for the shipment of radio-isotopes. They have 
maintained what is considered to be safe radia- 
tion protection standards for these shipments. 
For example, Oak Ridge National Laboratories 
solves this problem by specifying the type of 
shipping container, by limiting the radiation in 
the cab of a truck to a tolerance dose per day, by 
limiting the radiation level at the surface of a 
package to 50 mr/hr. for truck and rail shipment 
and 15 mr/hr. for plane shipment, by forbid- 
ding contamination on the surface of a package, 
etc. 

Ideally, perhaps, one would like to require that 
radio-isotopes be so packaged that no radiation 
escapes from or through the surface, but this 
criteria must be ruled out as impractical because 
it would make the packages cumbersome and 
heavy and would lead to excessive shipping 
expense. 

Recognizing the desirability of establishing 
shipping regulations that were safe and suitable 
to all interested parties, H. A. Campbell of the 
Bureau of Explosives called together, on May 
22, 1946, representatives from the film industry, 
the radium industry, the shipping industry, the 
Bureau of Standards, the Atomic Energy Proj- 
ects, the universities and hospitals, the Air 
Transport Association, the National Research 
Council, and the National Research Council of 
Canada. Many exchanges of letters between these 
representatives during the following months led 
to a further understanding of what the final 
shipping regulations should be like. 

All during this period the Committee on 
Radioactivity of the National Research Council 
(now the Committee on Nuclear Science of the 
NRC) had indicated an interest in this problem 
and on February 1, 1947, the Subcommittee on 
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Shipment of Radioactive Materials** was called 
together by the chairman, R. D. Evans of M.I.T. 
A number of other meetings of this committee 
have been held and the present form of the 
Regulations for the Shipment of Radio-Isotopes 
has been arrived at after many ramifications and 
a healthy exchange and development of ideas, 
and through the assistance of the Bureau of 
Explosives and with the valuable advice of such 
men as R. C. Gibbs, G. Failla, J. E. Rose, H. M. 
Parker, J. J. Quinn, T. C. George, and others. 
An act of Congress authorizes the Interstate 
Commerce Commission to formulate regulations 
governing the transportation of dangerous arti- 
cles, and the I.C.C. is authorized to call upon 
the Bureau of Explosives to maintain safe trans- 
portation of explosives and other dangerous 
articles. It is the responsibility of the Bureau of 
Explosives to assist in preparing and enforcing 
regulations governing the transportation of dan- 
gerous articles such as radioactive materials by 
rail freight, rail express, and by common con- 
tract and private carriers by motor vehicle. If 
penalization of a shipper is necessary, the Bureau 
of Explosives reports the case to the I.C.C. 
which directs it to the proper federal authorities 
for prosecution. The Post Office Department is 
responsible for setting up its own regulations 
governing (or prohibiting) radioactive shipments. 
It is not certain at present when the Civil Aero- 
nautics Board will publish regulations for radio- 
active shipments by air. For the present some of 
the airlines have agreed to carry shipments of 
radio-isotopes packaged in such a manner that 
the gamma-radiation does not exceed 14 mr at 


30 ft. Other airlines will not accept such ship- 


ments unless the surface radiation of the package 
does not exceed 15 mr/hr., and unless other 
presently. proposed regulations for rail and motor 
vehicle shipments are adhered to. Mr. R. D. 
Speas of the American Airlines and John Law- 
rence and Harold Crouch of the Eastman Kodak 


- Company have been working with R. D. Evans 


of M.J.T. since 1943, laying the ground work for 





** Members consist of R. D. Evans of M.I.T., L. F. 
Curtiss of the National Bureau of Standards, H. S. Youngs 
of Air Transport Association, R. B. Wilsey of Eastman 
Kodak Company, A. Morrison of National Research 
Council of Canada, W. G. McKenna of the Bureau of 
Explosives, and K. Z. Morgan of Oak Ridge National 
Laboratories. 
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the regulations for shipment of radio-isotopes by 
air. The final draft is being completed at present 
by H. S. Youngs of the Air Transport Associa- 
tion.* The United States Coast Guard is respon- 
sible for the transportation of radioactive ma- 
terials by water, and the regulations it establishes 
must not conflict with the Interstate Commerce 
Commission’s regulations when shipments are 
transferred to rail or motor vehicle. It is to be 
hoped that if and when these three sets of regula- 
tions are adopted, there will be maintained a 
minimum of conflict so as to minimize confusion 
when shipments are transferred from air to rail 
or to water. 

The regulations for the shipment of radio- 
active materials by rail and by highway have 
been prepared by this Shipping Committee of 
the National Research Council and submitted 
by the Bureau of Explosives to the Interstate 
Commerce Commission for final approval.‘ 
They can and should be amended from time to 
time when experience and further developments 
indicate the desirability of such changes. It 
would take a long time to give a complete dis- 


* This proposed Airline Tariff Rule for the Shipment of 
Radioactive Materials was distributed April 21, 1948 to 
the air carriers by Air Cargo, Inc. with the joint recom- 
mendation by the Air Transport Association of America 
that it be adopted. As soon as several air carriers have 
indicated approval, it will be submitted to the Civil Aero- 
nautics Board as the official airlines tariff. 

*These regulations were approved by the I.C.C. on 
October 24, 1947. 
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cussion of these regulations, so they are outlined 
briefly below, with particular emphasis on those 
regulations that apply to the shipment of radio- 
isotopes for research use. 


B. NEW REGULATIONS FOR THE SHIP- 
MENT OF RADIO-ISOTOPES BY 
RAIL AND MOTOR VEHICLE 


There are three groups of radioactive materials 
as follows: 

Group I.—Radioactive materials that emit 
gamma-rays only or both gamma- and elec- 
trically charged corpuscular rays. The radiation 
from these packages must be limited to 10 
units*** by the use of shileding and there must 
be no primary a- or §-radiation escaping from 
the shipping container. 

Group II.—Radioactive materials that emit 
neutrons and either or both types of radiation 
characteristic of group I. In this group the 
radiation limits are: 


1. y-radiation of 10 units. 
2. Electrically charged secondary corpuscular radiation 
that is the physical equivalentt of 10 units of y-radia- 


***1 unit=1 mr/hr. at 1 meter for hard y-radiation, 
or that amount of radiation which will have the same effect 
on sensitive film as 1 mr/hr. at a meter of hard y from 
radium filtered through 4 inch of lead. 

t The physical equivalent of a roentgen is that amount 
of radiation that would be absorbed in tissue to the extent 
of 83 ergs per gram. It is sometimes called the ‘Roentgen 
equivalent physical” or abbreviated as “‘rep.” 
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tion. None of the primary a@- and £-radiation may escape 
from the package. 


. Neutron radiation which is physically equivalent of 2 


units of gamma-radiation. 


. Combination of the above must not exceed the equiva- 


lent of paragraphs 1, 2, or 3. 


Group ITI. 


Radioactive materials that emit 


electrically charged corpuscular rays only, i.e., 


a or B, etc. None of the primary a- and @-radia- 


tion may escape from the package, and the 


secondary radiation must not exceed 10 mr/day 
at the surface of the package. 


Some of the more important specific regulations 


a. 


d. 


are as follows: 


If radioactive material of groups I, II, or III 
is liquid, it must be packaged with the use of 
a metal container and gaskets or with ab- 
sorbent material so that can leak to 
surface in case of an accident. (The former 


none 


part of this rule is understood to apply also 
to radioactive gas.) 

No loaded shipping container may have a 
y-radiation exceeding 200 mr/hr. or equiva- 
lent at any point of the readily accessible 
surface. 

The smallest permissible outside dimension 
of the shipping container is 4 inches. 

All packing must be such that the degree of 
fogging of undeveloped film with a 24-hour 
exposure at 15 feet is not greater than that 
produced by 11.5 mr of y from Ra filtered by 3 
inch of Pb. (The energy of such filtered Ra 
gamma is about 1.8 Mev.) 

Not more than 2 curies of Ra, Po, Pu, or 
other members of the heavy radio active series, 
or more of other materials which emit greater 
than 10" dis/sec., may be shipped in a single 
package. 


TABLE I, 


Distance in feet Distance in feet 


Potential 
(Mev) 


The package must have no significanttt sur- 


face contamination, and the shipping con- 

tainer must be of sturdy construction which 

meets certain specifications of the Bureau of 

Explosives. 

Packages are exempt from these regulations 

provided they meet the following specifica- 

tions: 

(1) There can be no leakage of radioactive 
material from the package. 

(2) The package does not contain more than 

0.1 me of Ra, Po, Pu, or 510° d/sec. of 

Sr*®, Sr®, and Ba'**, or 50 10° d/sec. of 

other less harmful radio-isotopes. 

There is no a, 8, or neutron radiation 

reaching the surface of the package, and 

the y-radiation is less than 10 mr/day at 

the surface. 


(3) 


Note: Empty shipping containers must conform to 
item f and to the three regulations above. 


The packages contain manufactured arti- 
cles, other than liquids, such as clock 
dials that are properly packed, and the 
radiation at the surface of the package is 
less than 10 mr/day. 

Special regulations are set up for the shipment 
of ores. They are exempt from these regula- 
tions, provided they are of low activity and 
are securely packed in leak-proof containers 
for rail freight shipment. The gamma-radia- 
tion must not exceed 10 mr/hr. at a distance 
of 5 ft. from any surface of the freight car. 
Materials that tend to remain fixed in the 
human body are classed as especially hazard- 
ous. Examples are Ra, Pu, Sr**, etc. These 
materials must be packed in special inside 
TABLE II.* 


Potential 


Multiplication 
c (Mev) 


Multiplication 
factor (K) : 


factor (K) 





0.035 21 
0.040 
0.045 


0.090 ¢ 
0.10 
0.15 


= 


20 
20 











to area that may 
be continuously 
occupied up to 8 
hrs. by passengers 
or employees. 


fotal Minimum distance 
number in feet to nearest 
of units* undeveloped film 


1 to 10 15 3 
11 to 20 20 4 
21 to 30 25 5 
31 to 40 30 6 


* See reference ttt. 
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to area that may 
be continuously 
occupied over 8 
hrs. by passengers 
or employees. 


0.050 
0.060 
0.070 
0.080 


* The values were obtained by R. B. Wilsey of Eastman Kodak 
Company, 


~_— ou 


20 0.20 
16 0.3 
12 0.5 
9 2.0 


pmo pt pt om om gy) 





and are probably the best available at present. 


tt “No significant” contamination is defined as <500 
d/min. measured over 100 cm? of surface for a and <0.1 
mr/hr. of y, or 0.1 mrep/hr. for 8. 
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metal containers or other containers approved 
by the Bureau of Explosives (as required by 
paragraph a). 

i. In case of fire, wreck, or breakage, great care 
should be taken to prevent radioactive ma- 
terials from getting on or inside the body. 
The Bureau of Explosives must be notified 
immediately. 

j. Special labels as obtainable from the Bureau 
of Explosives and shown in Figs. 1 and 2 must 
be properly filled in and placed on the outside 
of the package. The label in Fig. 1 is printed 
in red and the one in Fig. 2 is printed in blue. 

k. Table | indicates the minimum distances of 
radioactive materials permitted to films and 
persons for shipment by motor vehicle. In 
case of rail shipments, column 4 in Table | 
is omitted and column 3 is applied to distances 
from the dividing partition of a combination 
car as well as to distances to passengers and 
employees. The distance in Table | must be 
measured from the nearest point of the radio- 
active container or containers. Not more 
than forty unitsttt of radioactive material 
(red label) shall be transported in any car or 
stored in a depot at one time. If for any reason 
a package containing radioactive material 
(red label) remains in the same building for a 
period longer than 24 hours, it must be moved 
to a different location after each 24 hours. 


ti? 1 unit=1 milliroentgen per hour at one meter for 
hard gamma-radiation, or that amount of radiation which 
will have the same effect on film as 1 milliroentgen per hour 
at one meter of gamma-radiation from radium filtered 
through } inch of lead. 
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Experimentally, there is no perfect instrument 
for making these measurements. The Landsverk- 
Wollan electrometer shown in Fig. 3 is the best 
available in that it has the necessary sensitivity, 
it is not very energy dependent, and it will 
measure a-, §-, and y-radiation. The unit of 
radiation may be determined by making a 
measurement in mr/hr. at a meter from the 
nearest surface of the shipping box and then 


units of radiation = K X(mr/hr. at 1 meter) 


in which K is given approximately for gamma- 
radiation in Table II]. The value of K is unity 
for neutron and beta ray shipments. 

In spite of the large values of K in Table II 
for low energies, it is of interest to note that 
one need not use more than a ;g-in. sheet of lead 
to ship 2000 mc if the energy is less than 0.1 
Mev because of the large coefficient of absorp- 
tion of lead for this energy. In other words, 
Table II and the constant K can be ignored if 
one will always use a small amount of lead 
shielding for low energy sources. 

Table III for cobalt and radium was prepared 
by J. H. Roberson of Oak Ridge National Labo- 
ratories from experimental data he has collected 
and is typical of those tables being used for the 
shipment of radio-isotopes. This table is for 
cylindrical lead containers with a cylindrical 
cavity having a diameter equal to the height for 
holding the radio-isotope. The corners of the 
cylinders were rounded off to reduce the weight. 

The shield thickness required for beta-emitting 
radio-isotopes is given approximately by the 
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TABLE III. Wall thickness and weight of lead containers 
for shipment of radium and cobalt®. 


Maximum millicures Weight of lead container (Ibs) 








Wall : - ; ‘ 

thickness Ra in 0.5 diameter of cavity (inches) 

(inches) Co mm Pt. 5 1 2 
0.01 7.7 12 0.005 0.019 0.077 
0.1 8.5 14 0.066 0.22 0.83 
0.25 10 18 0.25 0.71 2.34 
0.5 13 26 0.86 1.98 5.65 
1.0 25 50 3.8 6.88 15.8 
1.5 48 92 10.1 15.9 31.8 
2.0 104 165 21.0 30.4 54.9 


| 


equation X =(0.21E,,—0.065)/d, in which X 
=thickness of material in inches to stop all 
B-rays, Em=maximum energy or f-rays in Mev, 
d=density relative to water (=g/cc). 

This equation cannot be used for energies less 
than 0.3 Mev, but in this case, as with a-emitting 
substances, a few layers of tin foil are sufficient. 
In shielding large B-ray sources, low atomic num- 
ber material should be used for the inner wrapper 
to reduce bremsstrahlen. 

Neutron sources should be shielded with ma- 
terial containing hydrogen (and perhaps iron 
and boron). Paraffin and water are not very 
satisfactory for shields because of their liquid 
state. No completely satisfactory meter is avail- 
able for these neutron measurements, but the 
double ionization chamber called Chang and Eng 
is best for fast neutron measurements and the 
boron lined electroscope is suitable for thermal 
neutron measurements. Figure 3 shows three 
instruments that might be recommended. 

All the above regulations for shipment by 
rail apply equally as well to motor vehicles. In 
addition, motor vehicles carrying radioactive 
material must be marked with a placard reading 
“‘Dangerous—radioactive material.’’ Also, no 
radioactive material may be carried on a bus 
when passengers are being transported unless no 
other practicable means of transportation is 
’ available. 
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TABLE IV. 





Minimum distance to 
films* (feet) 


Minimum distance to 
humans (feet) 





Total No. under 6 to over under 6to continuous 
of units 6 hr. 12 hr. 12 hr. 6hr. 12hr. exposure 

Up to 2 3 5 7 1 2 3 
3to5 5 8 10 2 3 4 
6 to 10 8 10 15 3 4 6 

11 to 20 10 15 20 4 6 8 

21 to 30 13 18 25 5 7 9 

31 to 40 15 20 30 6 8 11 








**Films” is interpreted to include mail, baggage and unidentified 
cargo. 


C. PROPOSED REGULATIONS FOR THE 
SHIPMENT OF RADIO-ISOTOPES BY AIR 


These proposed regulations are very similar to 
those given above for shipments by rail and motor 
vehicle. There are a few differences as follows: 


a. The limits for especially hazardous materials 
unless packaged so they cannot spill in case of 
an accident are: 

1. 100 me of Ra as a soluble salt or 300 mc when 
hermetically sealed in a small metallic capsul, 

2. 100 mc of Pu, Po, and similar hazardous a-emitters, 

3. Quantities of Sr, Sr”, Ba™, etc. that emit 410° 
d/sec. 

b. The limit of radiation exposure to animals is 
100 mrep/day. 

c. No photographic film should be shipped on 
planes carrying plane-load cargoes of radio- 
active materials. In any case, the exposure 
to films should be limited to 11.5 
equivalent) total exposure.f{ 


mr (or 


d. The gauge pressure of any gaseous shipment 
must not exceed 25 Ibs. per square inch at 
70°F. 

e. A different distance table is used 
shipments and is shown in Table IV. 


for air 


t It should be noted that no provision is made to protect 
film shipments that will be en route by plane, rail, or motor 
express for periods exceeding 24 hours. When such film 
shipments are to be made, special labels should be placed 
on the film package requesting that it not be in shipment 
with radioactive material. 
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Irradiation of Biological Materials by High Energy Roentgen Rays and Cathode Rays 


J. G. Trump, Electrical Engineering Department, AND R. J. VAN DE GRaaFF, Physics Department, 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received February 16, 1948) 


In this work roentgen rays and cathode rays of several 
million-volts energy have been applied to an investigation 
of their biological, photo-chemical, and germicidal effects, 
particularly as they are related to the processing of foods 
and biological materials. A constant-potential electrostatic 
generator, together with an acceleration tube, was used 
to produce continuous streams of electrons with homo- 
geneous and controllable energy. [R. J. Van de Graaff, 
K. T. Compton, and L. C. Van Atta, Phys. Rev. 43, 149 
(1943).] These high energy electrons were utilized both for 
the production of penetrating roentgen rays and for the 
direct irradiation of materials. The mechanism of the bio- 
logical action of both roentgen rays and cathode rays is 
discussed, as well as the energy considerations in their 
application to various absorbers. The companion paper 
[C. G. Dunn, W. L. Campbell, H. Fram, and A. Hutchins, 
J. App. Phys. 19, 605 (1948)] reports on measurements of 
the lethal action of these radiations on a wide variety of 
micro-organisms and also on the effect of the radiations on 
enzymes, vitamins, and certain whole food products. Both 
investigations have been conducted cooperatively by the 


Departments of Electrical Engineering and Food Tech- 
nology at the Massachusetts Institute of Technology. 
The systematic and quantitative study of the effects of 
roentgen and cathode rays on various elementary forms of 
living matter is regarded as an essential prelude to their 
possible widespread application to the preservation of 
foods and the sterilization of various biological materials. 
It was realized that the mechanisms of biological action 
of high energy roentgen rays and cathode rays are similar 
and depend closely upon the ionization energy absorbed 
per gram of material. The studies were begun with roentgen 
rays since these were immediately available and permitted 
quantitative measurement of the energy absorbed in the 
biological materials. Later, the water-cooled gold x-ray 
target was replaced with an aluminum cathode-ray window, 
and corresponding studies were begun on the direct appli- 
cation of high energy electrons to these materials. During 
the period of these biological studies the x-ray source was 
also used in a daily clinical program of deep cancer therapy 
under the medical direction of Dr. Richard Dresser and 
also in the investigation of the physical properties of high 
energy radiations. [R. Dresser, Radiology, in publication. ] 





ELECTROSTATIC SOURCE OF 
IONIZING RADIATIONS 


HE source of high energy radiation was a 

pressure-insulated electrostatic generator 
of the belt type at the Massachusetts Institute of 
Technology rated at 3-million volts constant 
potential.! This generator, shown in Fig. 1, is 
completely contained in a steel pressure tank 43 
feet in inside diameter and about 12 feet high, 
containing air at 13 atmospheres pressure. 

The high voltage terminal of the generator is a 
stainless steel hemispherical spinning about 24 
inches in diameter supported by an 88-inch high 
insulating column of metallic equipotential mem- 
bers separated by glass spacers. A single 3-ply 
butyl rubber cotton fabric belt 12 inches wide 
travels at 4000 feet per minute within the column 


‘J. G. Trump and R. W. Cloud, Am. J. Roentgenology 
and Rad. Therapy 49, 531 (1943). 
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and transfers negative electric charge continu- 
ously from ground to terminal. Although both 
sides may be utilized, it has been more convenient 
to spray charge merely on the upward run of belt. 
Under these conditions the continuous current 
available for the acceleration tube has been about 
500 microamperes at full voltage. 

The generated voltage may be controlled from 
values of a few hundred thousand volts to a maxi- 
mum of about 4 million by regulating the amount 
of charge sprayed on the belt for a given fixed 
setting of the electron current in the tube. The 
negative voltage of the terminal is divided uni- 
formly lengthwise the column by means of a 
variable corona gap between each of the sepa- 
rated equipotential planes. 

The acceleration tube is mounted within the 
generator column and in a position parallel to the 
charge conveyor belt. The source of electrons is 
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Fic. 1. 3-million volt pressure-insulated electrostatic 
accelerator with tank removed. 


a hot tungsten filament supplied with power from 
a small permanent-magnet generator contained 
within the high voltage terminal and driven by 
the upper pulley. The insulating column of the 
tube is a rigid self-supporting structure made of 
alternate glass rings and flat stainless steel disks 
sealed together by a thermoplastic-film technique 
developed in this laboratory. Each disk has an 
axial 3-inch diameter hole and is electrically con- 
nected to the equipotential plane of the generator 
column at that level. The electrons are thus 
progressively accelerated as they pass from the 
filament axially downward through the planes of 
the diaphragms. The ground potential end of the 
tube extends on through the tank base into the 
room below and terminates in the target arrange- 
ment shown in Fig. 2. X-rays are produced when 
the electrons impinge on the water-cooled gold 
‘disk, which is 0.25 inch thick. The electrostatic 
focusing action of the tube diaphragms brings the 
electrons to the target in a focal spot about 10 
millimeters in diameter. This may be brought to 
less than 0.25 millimeter by passing d.c. current 
through an axial magnet coil at the grounded end 
of the tube. By withdrawal of the gold disk from 
the path of the stream the electrons may be 
allowed to pass through the 0.002-inch aluminum 
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window into the atmosphere as a high energy 
cathode-ray stream, as shown in Fig. 3. This 
target arrangement has permitted great flexibility 
in working alternately with x-rays and cathode 
rays since it eliminates the necessity for breaking 
the vacuum and thereafter reconditioning the 
electrode surfaces. The acceleration tube is con- 
tinuously evacuated by a 
pumping system. 


mercury diffusion 

In the roentgen-ray studies the biological 
material was placed in a thin metal or glass con- 
tainer located an average distance of 3 centimeters 
below the x-ray source. The dosage rate at this 
position could be adjustéd to a maximum of 
about 5000 roentgens per second. A small special 
ionization chamber was mounted directly under 
the sample and connected to an integrating cir- 
cuit to measure the total ionization dose received 
by the material. The total doses applied to 
samples ranged from 10-thousand to 2-million 
roentgens, a roentgen being defined as that 
amount of radiant energy which, when passing 
through 1 cubic centimeter of air at standard 
temperature and pressure, will release 1 electro- 
static unit of charge. Figure 4 shows containers 
of 0.005-inch aluminum and 0.01-inch stainless 
steel used in this work, together with the carrier 
for sample and ionization chamber which could 
be automatically moved into or out of the beam 
of radiant energy. 


NATURE OF THE BIOLOGICAL REACTION 
TO RADIATION 


The photo-chemical and biological reactions 
produced by the absorption of x-rays and cathode 
rays are similar in their physical nature and re- 
quire closely equivalent energies. When matter is 
traversed by x-rays, energy is absorbed from the 
x-ray photons by the photoelectric, Compton 
scattering, or pair-production mechanisms, and 
electrons are accelerated with related amounts of 
energy. These energetic secondary electrons ex- 
cite and ionize many other atoms in their path 
and may thereby produce low energy tertiary 
electrons still capable of excitation and ioniza- 
tion. Chemical changes in an absorber are 
brought about primarily by the ionization and to 
a lesser extent by the excitation of the constituent 
atoms. These higher energy states generally favor 
the dissociation of complex molecules, though 


JOURNAL OF APPLIED PHYSICS 














| ELECTRON BEAM 





/ | ED 
, y \ BEAM - DEFINING 
VA L, AL STEEL BLOCK 





° ~ . \ x Y 
Fic. 2. Combined x-ray target NN \ , /, ' WOQQQAR ada 
and cathode-ray window. ‘ N J, Af . . /s 



















































t 
: \WATER LINES 
SS E 


\ 025" GOLD ay A % > RETRACTABLE WATER-COOLED 
\ % GOLD TARGET 
INSULATED STEEL O1SK = 


> a7 «| «.002" ALUMINUM WINDOW 
5 CC SAMPLE HOLDER 
RETRACTABLE X-RAY IONIZATION CHAMBER 
SAMPLE HOLDER 


COMBINED X-RAY TARGET and CATHODE RAY WINDOW 








combination is also often made possible, as in the ionization of other atoms. The biological effects 
formation of hydrogen peroxide by the irradiation are produced by these excitation and ionization 
of water. The ionizing energy is imparted so events. 
directly to the absorber that profound chemical 
changes may often be accomplished with negli- 
gible gross heating effects. In the case of living 
tissue, injury—and in the limit death—of a cell 
is the result of physiochemical changes brought 
about by direct absorption of ionizing energy or 
indirectly by the reaction of surrounding irradi- 
ated tissue. Much is now known concerning the 
influence of dosage rate, the linear ion density 
along the track of the high energy particle, and 
the influence of the surrounding medium on the 
complex response to radiation of molecules, cells, 
and microorganisms. Ultimately, the physio- 
chemical and biological reactions depend on the 
number of ionizing electrons produced in the 
passage of the radiation through the material. 
In cathode-ray irradiation electrons are di- 
rectly projected into the material and produce 
ionization similar in its general characteristics to 
that produced by x-rays. Each high energy 
electron in the cathode-ray stream proceeds into 
the material, losing energy by collision and 
causing the excitation or ionization of thousands 
of atoms in its path. These primary electrons 
thus distribute the energy of the cathode-ray 
stream through the volume of the absorber. 
Many of the secondary electrons produced in 
these encounters may themselves possess suffi- 
cient energy to act as biological agents by the 


In general the primary difference between the 
action of supervoltage x-rays and cathode rays on 





Fic. 3. 3-million volt cathode-ray beam. 
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Fic, 4. Retractable car- 
riage with ionization cham- 
ber and sample holders for 
x-irradiation. 








living matter is in the gross distribution of the 
ionization energy within the absorber. Whereas 
cathode rays may penetrate into water a maxi- 
mum distance of about 1 centimeter for each 
2-million volts of energy and produce no effect 
beyond this maximum range, x-rays are relatively 
penetrating and diminish exponentially in pass- 
ing through an absorber. Thus the intensity of a 
2-million volt roentgen-ray beam diminishes 
almost one-half for each 10 centimeters of water 
through which it passes.” In the case of initially 
parallel cathode rays of homogeneous energy, the 
maximum ionization density produced in an 
absorber occurs at about one-third the maximum 
range.* This location of the region of maximum 
ionization toward the forward part of the maxi- 
mum range for any given voltage is due to the 
naturally high scattering tendencies of electrons 
and would be accentuated if the cathode-ray 
beam were accelerated by pulsed or varying 
voltages. The insert of Fig. 5 compares the gross 
distribution of ionization produced by a beam of 
3-million volt x-rays traversing water with that 
produced by cathode rays passing through the 
same material. 


ENERGY ABSORBED AND TEMPERATURE RISE OF 
: IRRADIATED MATERIAL 


The amount of energy absorbed in a gram of 
air when 1 million roentgens of x-radiation pass 
through it can be measured and calculated. 
Irradiation of air with a dose of 1 million 
roentgens is equivalent to the absorption of 8.5 


2D. E. Lea, Actions of Radiations on Living Cells, The 
Macmillan Company, New York, 1947). 

3J. G. Trump, R. J. Van de Graaff, and R. W. Cloud, 
Am. J. Roentgenology and Rad. Therapy 43, 728 (1940). 
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joules of ionizing energy per gram of the ab- 
sorber.* Over a wide range of wave-lengths this 
relation between energy absorbed per gram and 
x-ray dose in roentgens also holds fairly closely 
for water which has constituents of atomic 
number close to those of air. Dependence upon 
the atomic number of the absorber becomes pro- 
nounced for radiation of less than 100 kilovolts 
energy because of the photoelectric absorption 
process and above several million volts where 
pair production becomes an influential absorption 
mechanism. For the energy region between, in 
which the Compton scattering process predomi- 
nates, the energy absorption per gram is nearly 
independent of atomic number. 

In another investigation* we have found, as 
might be expected, that this relation likewise 
holds when cathode rays are directed at an 
absorber. In our physical measurements close 
equivalence has been found between the number 
of ion pairs produced by the absorption of a given 
amount of roentgen-ray energy and the same 
amount of cathode-ray energy. It is therefore 
reasonable to expect that experimental studies of 
the biological and photo-chemical effects of x- 
rays conducted under conditions of high physical 
control can usually be correlated to closely 
corresponding effects obtained with equal 
amounts of cathode-ray energy. 

The energy absorbed by the passage of 
roentgen-ray or cathode-ray energy through a 
material ultimately appears as heat and is mani- 
fested by a corresponding temperature rise. The 
temperature rise in water produced by an x-ray 
dose of 1 million roentgens or by the absorption 


4L. H. Gray, Proc. Roy. Soc. A156, 578 (1936). 
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of 8.5 joules per gram is 2.04°C. This temperature 
rise is directly proportional to the x-ray dose and 
inversely proportional to the specific heat of the 
absorbing material. Since 6000 roentgens is a high 
cancericidal dose and 1 million roentgens is a 100 
percent bactericidal dose, it is evident that pro- 
nounced biological and photo-chemical effects are 
realized with amounts of radiation which produce 
a relatively unimportant temperature rise in the 
absorbing material. The energy of these pene- 
trating radiations is not expended directly in 
raising the average thermal energy of the mole- 
cules as in the case of heat sterilization, but is 
used more efficiently to induce by the process of 
excitation and ionization discrete changes in the 
atomic and molecular structure of the irradiated 
material. 


CHARACTERISTICS OF CATHODE RAYS 
FOR IRRADIATION 


The advantage of cathode rays over x-rays in 
the irradiation of living materials lies primarily in 
the high inherent efficiency of this process and 
in the controllability of the cathode-ray dosage 
under proper conditions. The efficiency of util- 
izing the electron beam energy can ordinarily be 
several hundred times higher in cathode-ray 
irradiation than when x-rays are used. In utiliza- 
tion of the cathode rays accelerated in a constant- 
potential vacuum tube, an over-all efficiency 
approaching 100 percent is realizable. All of the 
emergent cathode-ray beam energy may be di- 
rected into the absorber and expended in ioniza- 
tion processes within the well defined maximum 
electron range, which increases linearly with 
voltage. The efficiency and controllability of the 
electrostatic cathode-ray process makes possible 
the continuous delivery of ionization doses 
equivalent to millions of roentgens in a small 
fraction of a second. 

With x-radiation, on the other hand, most of 
the kinetic energy of the electrons within the 
acceleration tube is lost as heat at the target. At 
100 kilovolts only 0.2 percent of the energy of the 
electron stream is converted into x-radiation. At 
2-million volts about 5 percent is transferred into 
radiant energy. Moreover, since x-rays are pene- 
trating and proceed in all directions, only a small 
fraction of the x-ray energy can readily be ab- 
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Fic. 5. Distribution of ionization in water by 0.5-, 1-, 2-, 
and 3-million volt cathode rays and by 3-million volt 
x-rays. 


sorbed uniformly within the desired volume of 
the selected absorber.® 

The cathode-ray dosage falling on an absorber 
can be completely specified with great exactness 
as it is determined by the number of impinging 
electrons per second and their accelerating volt- 
age. With the constant-potential electrostatic 











Fic. 6. Compact 2-million volt x-ray or cathode-ray source 
using pressure-insulated electrostatic generator. 


rr J. G. Trump, C. R. Moster, and R. W. Cloud, Am. J. 
Roentgenology and Rad. Therapy 57, 703 (1947). 
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TABLE I. 
104r dose 10% dose 
Rate per (0.85 joules /gram ) (8.5 joules /gram) 
Second 500 grams 50 grams 
Minute 30 kilograms 3 kilograms 
Hour 1800 kilograms 180 kilograms 
Day $3,000 kilograms 4300 kilograms 


300 24-hr. days 13,000,000 


kilograms 


1,300,000 
kilograms 


accelerator the ionizing power of the incident 
cathode rays can be accurately determined by 
direct measurement of the voltage and current. 
Thus the ionizing power of the gompact 2-million 
volt generator shown in Fig. 6 operating with an 
electron beam current of 250 microamperes is 500 
joules per second, or 500 watts. Applied to the 
sterilization of water-equivalent materials, such a 
cathode-ray source can produce a sterilizing dose 
of 8.5 joules per gram (equivalent to 1 million 
roentgens) in about 50 grams of material per 
second. Table | gives the calculated sterilizing 
rate in continuous-processing applications for the 
cathode-ray source shown in Fig. 6 and operating 
with 500 watts output. 

The distribution of ionization produced in 
water by a cathode-ray beam of 1-, 2-, and 3- 
million volt energy is shown in Fig. 5. For the 
2-million volt electrons, for example, the maxi- 
mum range in water is seen to be close to 10 
millimeters, and the maximum ionization density 
occurs in a broad-thickness region centering at a 
depth of 3 millimeters. For many continuous 
processes a 2-million volt stream of cathode rays 
would be capable of irradiating effectively a 
stream thickness of about 6 millimeters. Two 
such cathode-ray sources irradiating from oppo- 
site sides would deliver a fairly uniform dose to 
approximately twice this thickness. Correspond- 
ingly higher voltages would be required for the 
irradiation of greater thicknesses, and allowance 
. must be made for absorption by an intervening 
wall materials. 

When the cathode-ray beam emerges into the 
atmosphere from the acceleration tube, some 
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divergence of the beam results from the normal 


scattering by the window diaphragm. Hence the 


area which is irradiated by the beam increases 
with distance from the cathode-ray window. The 
cathode-ray stream may also be magnetically 
diverged. Thus the depth of penetration of 
cathode rays is determined by the voitage and the 
density of the absorbing material, while the area 
which is irradiated can be controlled by selection 
of window-absorber distances, by magnetic de- 
flection, and by suitable scattering arrangements. 


CONCLUSION 


High energy cathode rays and x-rays now offer 
a new and effective means of producing biological, 
bactericidal, and even chemical changes on a 
practical scale. These changes are brought about 
by excitation and ionization of constituent atoms 
of the absorber and need relatively small amounts 
of ionizing energy. Cathode rays produce similar 
reactions and can be utilized with several hundred 
times the efficiency of x-rays. The temperature 
rise evoked by cancericidal doses of these radia- 
tions is about 0.01°C, while complete sterilization 
of bacteria can be accomplished with ionization 
doses resulting in a 2°C rise. The constant- 
potential output of electrostatic accelerators is 
well suited to quantitative radiation studies and 
to large-scale sterilization and_ inactivation 
processes. 
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Biological and Photo-Chemical Effects of High Energy, Electrostatically Produced 
Roentgen Rays and Cathode Rays 


Cecit G. DuNN, Wit_iam L. CAMPBELL, HARVEY FRAM, AND ARDELIA HUTCHINS 
Department of Food. Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received February 16, 1948) 


INTRODUCTION 


N a companion article by Trump and Van de 

Graaff the physical characteristics of electro- 
static accelerators for the production of high 
energy roentgen rays and electrons are discussed, 
together with factors influencing their application 
to the irradiation of biological materials. The 
principal portion of the research reported herein 
is concerned with the lethal action of x-rays pro- 
duced at high voltages on pure cultures of 
microorganisms and those found in spices, apple 
juice, milk, soil, water, and catgut sutures. Some 
preliminary data on cathode rays are included. 


MICROBIOLOGICAL INVESTIGATIONS 
WITH X-RAYS 


Literature 


A number of scientists have studied the lethal 
and sublethal effects of x-rays (particularly those 
produced at low or comparatively low voltages) 
on microorganisms. References to some of the 
more important of these papers are given at the 
end of this report. 

Proctor, Van de Graaff, and Fram (1942) 
carried out research concerning the action of 
x-rays produced at high voltage on finely ground 
beef and observed that the bacterial counts were 
considerably reduced asa result of the irradiation. 
The present research was suggested as a result of 
these findings. 


Microorganisms Used 


Of microorganisms investigated to date the 
main emphasis has been on bacteria. Those ex- 
amined consisted of food spoilage, heat-resistant, 
Gram-positive, Gram-negative, spore forming, 
non-spore forming, and pathogenic organisms. 
The yeasts included both active and dormant 
cultures; the molds, two of the more commonly 
occurring species. 

The microorganisms used were obtained from 
the stock collection of the Department of Food 
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Technology, M.I.T.; the American Can Com- 
pany ; and the National Canners’ Association. 


Procedure 


The lethal effects of high voltage x-rays on pure 
cultures of microorganisms were determined es- 
sentially as follows. The organism was grown in a 
medium (liquid or solid) and at a temperature 
conducive to the development of a large number 
of cells or spores. A uniform suspension of the 
organism was prepared either in a broth medium 
or in sterile distilled water. When necessary, the 
culture was shaken with sterile glass beads in 
order to obtain a uniform suspension. Occasion- 
ally the supernatant portion of the suspension 
was transferred aseptically to a second sterile 
tube before use. The suspension thus prepared 
was pipetted into the sterile glass or stainless 
steel containers which were used for irradiation 
purposes. The containers with their contents 
were placed successively in the x-ray field and 
subjected to various roentgen doses, varying 
from as low as 10,000 to as high as 8,000,000 
roentgens. At the time of irradiation, data were 
recorded concerning the dosage in roentgens, the 
voltage in megavolts, the electron beam current 
in microamperes, and the exposure time in 
seconds or minutes. The irradiated samples were 
immediately plated in a solid medium suitable for 
the growth of the organism. The plates were 
incubated at temperatures favorable for the 
growth of the survivors and the visible colonies 
were counted with a Quebec Colony counter 
after a suitable incubation period. Details con- 
cerning some of the experimental data are pre- 
sented in Tabie | of this report. 

Suspensions containing viable bacterial spores 
were prepared by heating suspensions of the 
organisms being investigated at 80°C for 10 
minutes or at the boiling temperature for five or 
more minutes. The suspensions thus contained 
the viable bacterial spores and the heat-destroyed 
vegetative cells and weaker spores. 
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The spices examined were purchased on the _ taining sterile glass beads and dilutions prepared 
open market. One-gram samples were used for therefrom. Aliquot portions of various dilutions 
irradiation purposes. Control and irradiated were plated in duplicate in Bacto nutrient agar 
samples of the spices were shaken with sterile and incubated at 30°C for three to six days before 
distilled water in wide-mouth glass bottles con- final observations were made. 


TABLE I. Some experimental data. 














Tem : 
pera- Age of 
ture at cul- Composi- Medium in Incuba- _Incuba- 
Medium which ture tion of which bacteria tion tem- tion 
Source of in which grown, (in irradiation suspended for Plating perature period 
Culture culture grown a days)! container irradiation medium —_ (in days) 
Aerobacter Stock Bacto 37 1 Glass Nutrient Bacto 37 1 
aerogenes collection* nutrient broth nutrient 
broth agar 
Escherichia Stock Bacto 37 1 Glass Nutrient Bacto 37 1 
coli collection nutrient broth nutrient 
broth agar 
Pseudomonas Stock Bacto 37 2» Stainless Nutrient Bacto 37 2 
aeruginosa collection nutrient 3 steel broth nutrient 
broth agar 
Pseudomonas Stock Bacto 24-26 1» Stainless Nutrient Bacto 24-26 3-44 
fluorescens collection nutrient 3 steel broth nutrient 
broth agar 
Serratia Stock Bacto 30 Zz Stainless Nutrient Bacto 26 2 
marcescens collection nutrient steel broth nutrient 
agar agar 
Sarcina flava Stock Bacto 30 1® Stainless Nutrient Bacto 30 4.5 
collection nutrient 2 steel broth nutrient 
broth agar 
Staphylococcus Stock Bacto 37 1 Glass Nutrient Bacto 37 1 
aureus collection nutrient broth nutrient 
broth agar 
Bacillus Manioca Bacto Stainless Nutrient Bacto 26 3 
mesentericus, sirup nutrient steel broth nutrient 
strain broth agar 
B. sterothermo- American Dextrose 55 1 Stainless Dextrose Bacto 55 2 
philus Can tryptone steel tryptone dextrose 
Company broth broth tryptone 
agar 
B. thermoaci- National Special 37 5® Glass Distilled Special 52-54 1-44 
durans Canners proteose 7 water proteose 
(No. 43-P) Association peptone peptone 
agar® agar® 
Canco 6B American Bacto 37 3® Glass Distilled Bacto 37 2 
Can nutrient 5 water nutrient 
Company agar ‘Nutrient agar 
broth 
Flat sour National Bacto 53 3 Glass Distilled Bacto 53 3-44 
" (No. 1518) Canners dextrose 4 water dextrose 
Association tryptone tryptone 
agar agar 
Spore former Catgut Bacto 37 2° = Stainless Nutrient Bacto 37 3 
suture nutrient 4 steel broth nutrient 
. agar agar 


Department of Food Technology, Massachusetts Institute of Technology. 

Some tests carried out using cultures of one age, some using cultures of another age. 
Refer to Food Research 7, 186 (1942). 

The incubation period varied between these limits, depending on the test. 

Held at room temperature for several days in addition. 

At time of irradiation. 
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Molds Spices Miscellaneous items 


Fic. 1. Lethal effects of x-rays on microorganisms. (Roentgen doses versus percentages destroyed. ) 
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Organism 


Aerobacter aerogenes 
Escherichia coli 


Pseudomonas 
aeruginosa 


Pseudomonas 
fluorescens 


Serratia marcescens 


Sarcina flava 
Staphylococcus 
aureus 


Bacillus mesentericus, 
strain 


Bacillus 


sterothermophilus 


Bacillus subtilis 


Bacillus 


thermoacidurans 


Canco No. 6B 


‘Flat sour (No. 1518) 


Spore former from 
catgut suture 


Mor- 
phology 


Rods 


Rods 


Rods 


Rods 


Rods 


Cocci 


Cocci 


Rods 


Rods 


Rods 


Rods 


Rods 


Rods 


Rods 


Gram 
reaction 


Negative 
Negative 


Negative 


Negative 


Negative 


Positive 


Positive 


—s 
I ositive 


Negative 


Positive 


Positive 


Positive 


Positive 


Positive 








* The sample contained 40.32 percent of spores. 


» A mixture of spores and vegetative cells. 
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Relation to 
oxygen 





Aerobic, 
facultative 


Aerobic, 
facultative 


Aerobic, 
facultative 


Aerobic 


Aerobic 
facultative 


Aer« bic 
Aer« bic 
facultative 
Aerobic 


facultative 


\erobic 
facultative 


Aerobic 
facultative 


Aerobic 
facultative 


Aerobic, 
facultative 


Aerobic, 
facultative 


Aerobic, 
facultative 


Sporu- 
lation 


None 
None 


None 


None 


None 


None 


None 


Spc Tes 


Sp. res 


Spores 


Spores 


Spores 


Spores 


Spores 


Pathogenicity 





None 
None 


Pathogenic 


to some small 


animals 


None 


None 


None 


Pathogenic 


None 


None 


None 


None 


None 


None 


None 


Number ot 
organisms 
per ml 
238,000,000 
700,000,000 


760,000,000 


400,000,000 


. 


2,090,000 ,000 


337,000,000 


990,000,000 


380,000,000 


2,400,000 


1,000,000,000 


6,200,000* 


65,000,000 


75,000,000 


4,800,000" 


240,000,000" 


44,500,000° 


400,000,000" 


420,000,000" 


2,200° 


1,100,000° 


22,000,000» 


13,000,000 


TasBie Il. Summary of data on lethal 


Percent- 
age de- 
stroyed by 
25,000 


99.44 


99,999 + 
99.39 


98.68 


94.60 4 


26.7 


92.91 


95.60 
81.8 
96.15 


46.4 


roentgens roentgens 





Percent 
age de 
stroyed by 
35,000 


92.86+ 


99.999 + 


99,999 + 


99.97 


64.58 
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Percent- 
age de- 
stroyed by 
100,000 


roentgens 
99.99 


99.97 + 


100.00 
99.98 +- 


99.999 + 


99.998 


99.0 


99.6 


98.1 

99.90 
95.41 
99.97 


71.0 





effects of x-rays on selected bacteria. 


Percent 
age de 
stroyed by 
140,000 
roentgens 


100.00 


100.00 


99.999 +- 


99.88 + 


Percent 
age de- 
stroyed by 
500,000 


roentgens 
100.00 


100.00 


100.00 


100.00 


100.00 
99.999 +4 


100.00 


100.00 


99.97 


99.99 


99.98 


99.997 + 
99,999 + 
99.996 
99.999 + 
100.00 
99.999 +- 
98.0 


99.6 


Percent 
age de 


stroyed by 
700,000 


roentgens 


Number of roentgens required 
for destruction 


100.00 Less than 250,000 _ 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


but more than 100,000 


Less than 150,000 
but more than 90,000 


Less than 140,000 
but more than 35,000 


Less than 140,000 
but more than 35,000 


Less than 350,000 
but more than 140,000 


Less than 100,000 
but more than 25,000 


Less than 500,000 
but more than 250,000 


Less than 350,000 
but more than 140,000 


Less than 250,000 
but more than 100,000 


Less than 1,500,000 

but more than 1,000,000 
Not destroyed by 
1,000,000 


Destroyed by 1,500,000 
but not by 1,000,000 


Less than 1,000,000 
Less than 1,000,000 
but more than 500,000 


Less than 1,000,000 
but more than 500,000 


Not destroyed by 
500,000 


Less than 1,000,000 
but more than 500,000 


Less than 500,000 
but more than 250,000 


Less than 1,000,000 
but more than 500,000 


Not destroyed by 
1,500,000 


Less than 2,000,000 
but more than 1,500,000 


Age of 
culture 


1 day 
1 day 


2 days 


3 days 
1 day 

3 days 
2 days 
1 day 
2 days 


1 day 


1 day 

7 days 

7 days 

3 days 

5 days 

4 days 
3+ days 
4 days 


5+ days 





Special comments 


Of sanitary significance 


Produces Fluorescence 


Produces Fluorescence 


Produces red pigment 


Produces yellow pigment 
Produces orange pigment 


Strain especially resistant 
to heat; isolated from 
manioca starch 


Produces flat-sour 
spoilage of tomato juice 


Causes spoilage of canned 
foods 


Source—American Can 
Company 


Causes flat sour spoilage 
Source—National 
Canners Association 


Resistant organism iso- 
lated from catgut suture 





* All spores-—no viable vegetative cells. 


“Only one organism survived. In three other tests, the organism failed to survive the application of 500,000 roentgens. 
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TABLE III. Summary of data concerning lethal effects of x-rays on yeasts. 











Percentages of yeasts destroyed by 








25,000 100,000 250,000 500,000 1,000,000 1,500,000 Number of roentgens 
Designation of yeast roentgens roentgens roentgens roentgens roentgens roentgens required for destruction 
Saccharomyces cerevisiae 85.70 99.98+ 99.999 + 100.00 100.00 Less than 500,000 but 
strain 1 more than 250,000 
S. cerevisiae 86.20 98.62 99.966 100.00 100.00 Less than 500,000 but 
strain 2 more than 250,000 
S. cerevisiae 58.98 99.45 99.998 + 100.00 100.00 Less than 500,000 but 
strain 3 more than 250,000 
S. cerevisiae 93.75 99.824 99.998 + 100.00 100.00 Less than 500,000 but 
strain 4 more than 250,000 
S. cerevisiae 87.70 9995 99.997 99.999 + 100.00 Less than 1,000,000 but 
strain 5 more than 500,000 
S. cerevisiae 96.53 97.92 99.999+ 100.00 Less than 1,500,000 but 


(dry state) strain 4 


Torulopsis pulcherrima 98.75" 99.99+  99.999+ 


Torulopsis rosea 67.30 97.98 99.90 
17.64 96.41 99.95 








* Received 28,500 roentgens. 


Irradiated and control samples of raw and 
pasteurized milks, apple juice, soil, and waters, 
after appropriate dilution with sterile distilled 
water, were plated in Bacto nutrient agar and 
incubated at approximately 25°C for two to five 
days before making final observations. 

Vials containing catgut sutures rolled on a 
spindle and partly submerged in liquid were sub- 
jected to irradiation by the high voltage x-rays. 
The tubes were dipped in alcohol and flamed 
following irradiation and broken at scored points. 
The catgut sutures were removed aseptically and 
placed in large culture tubes containing nutrient 
broth enriched by the addition of 0.5 percent 
sodium chloride and 0.1 percent of dextrose. 


TABLE IV. Lethal effects of x-rays on Aspergillus niger.* 








Average 
percentage 
Dose in destruction Range of percentage 
roentgens of molds destruction of molds 
25,000 95.96 91.29 to 98.25 
50,000 99.63 98.93 to 99.96 
100,000 99.98 99.954 to 99.999 
250,000 99.999 99.999 to 100.00 
500,000 100.00 


100.00 to 100.00 


* The results are based on a minimum of five separate tests. The 
ages of the cultures irradiated varied from 2 to 28 days. The molds 
were grown on Bacto Sabouraud’s dextrose agar and suspended in 
Bacto malt extract broth for irradiation. 
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100.00 100.00 


more than 1,000,000 


Less than 500,000 but 
more than 250,000 


99.999 More than 500,000 


99.999 + 100.00 Less than 1,000,000 but 
more than 500,000 


Observations were made after incubation for one 
week at 37°C. 


Observations and Results 


Many of the observations and results con- 
cerning the lethal effects of x-rays produced at 
high voltage are summarized in Tables II to VI, 
inclusive, and in Fig. 1. 


Discussion of Results 


An examination of Table II reveals a number 
of interesting facts. The non-spore forming bac- 
teria, both Gram-negative and Gram-positive, 
were destroyed by doses of less than 500,000 
roentgens and from 64.5 to 99.999 percent of them 
were destroyed by 35,000 roentgens. Pseudomonas 
aeruginosa and Ps. fluorescens were particularly 
susceptible to the lethal action of x-rays, the 
former being destroyed by less than 140,000 
roentgens. Bacterial spores, however, were con- 
siderably more resistant to the action of x-rays 
produced at high voltage than the vegetative 
cells and generally required more than 500,000 
roentgens for destruction of all cells. The most 
resistant spore-former investigated required more 
than 1,500,000 but less than 2,000,000 roentgens 
for destruction. However, under the conditions of 
the tests, most spores were destroyed by the 
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TABLE V. Summary of data concerning lethal effects of x-rays on the microorganisms in some spices. 











Number of Percentage destroyed by 
organisms 35,000 140,000 350,000 700,000 1,400,000 Number of roentgens 
Spice per gram roentgens roentgens roentgens roentgens roentgens required for destruction 
Allspice, ground 740,000 39.18 93.37 99.45 99.998 100.00 Less than 1,400,000 but 
more than 700,000 
Pepper, white, ground 8,600 45.34 97.79 99.88 99.96 100.00 Less than 1,400,000 but 
more than 700,000 
Sage, ground 84,000 42.85 83.33 99.70 99.994 100.00 Less than 1,400,000 but 
, more than 700,000 
Percentage destroyed by 
25,000 100,000 250,000 500,000 1,000,000 
roentgens roentgens roentgens roentgens roentgens 
Mace 2,200 27.27 93.18 97.72 99.77 100.00 Less than 1,000,000 but 


more than 500,000 





application of less than 1,000,000 roentgens and 
from 15 to 96 percent were destroyed by 25,000 
roentgens. 

The yeasts (Table I11) examined appeared to 
be slightly more resistant on the average to x- 
radiation than the average non-spore forming 
bacteria. Destruction of all cells of four strains of 
Saccharomyces cerevisiae occurred with the appli- 
cation of less than 500,000 roentgens and from 
58.9 to 93.7 percent of the cells of five strains of S. 
cerevisiae were destroyed by 25,000 roentgens. It 
required more than 500,000 but less than 
1,000,000 roentgens to effect complete destruc- 
tion of a strain of S. cerevisiae in the dry state. 
Two species of Torulopsis were destroyed by less 
than 1,000,000 roentgens. 7. rosea appeared to be 
somewhat more resistant to x-rays than 7. 
pulcherrima. 





Research carried out on two species of molds 
indicated that considerable variation in sus- 
ceptibility to x-rays may be expected amongst 
different species. In four out of five tests 
(Table IV), A. niger was destroyed by the appli- 
cation of 250,000 roentgens. Both the youngest 
(2 days old) and the oldest (28 days old) cultures 
of A. niger were destroyed by this dosage. In five 
tests carried out with a mold of a species of the 
genus Mucor, using essentially similar procedures 
as for A. niger, a dose of 1,000,000 roentgens 
destroyed all of the individual molds, whereas 
500,000 roentgens destroyed over 99 percent of 
them. The species of Mucor was thus considerably 
more resistant to x-rays produced at high voltage 
than the strain of A. niger used. 

The microorganisms associated with ground 
allspice, mace, white pepper, and sage (Table V) 


TABLE VI. Summary of data concerning lethal effects of x-rays on microorganisms found in miscellaneous items. 











Number of 
organisms per 
ml or per gram 





Percentage destroyed by 











of unirradiated 25,000 100,000 250,000 500,000 1,000,000 Number of roentgens 

Item sample roentgens roentgens  roentgens roentgens roentgens required for destruction 
Apple juice 440,000" 99.47 99.96 99.999 +- 
Milk, raw 1,550,000 96.67 99.67 99.998 99.999 +- 100.00 Less than 1,000,000 but 

' more than 500,000 
Milk, pasteurized 15,000 47.33 99.82 99.98 99.986 100.00 Less than 1,000,000 but 
more than 500,000 
Soil 8,000,000> 91.7 98.9 99.4 100.00 Less than 1,000,000 but 
14,000,000 87.0 98.0 99.9 100 00 more than 500,000 

Water, tap 146 100.00 Less than 250,000 
Water, river 12,000 99.91 100.00 100.00 100.00 Less than 250,000 but 


more than 100,000 








* The apple juice contained 200,000 molds per ml and 240,000 bacteria and yeasts per ml. No molds survived the application of 1,000,000 


roentgens. 
> Organisms per gram of sample. 
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TaBLe VII. Summary of data concerning the lethal effects on some microorganisms of cathode rays 





Tem- 
pera- Age ot 
ture at cul- 
which ture 
Medium in grown, (in 
Culture which grown* sn days)* 
Staphylococcus aureus Bacto 37 1 
nutrient 
agar 
Escherichia coli Bacto 37 1 
nutrient 
agar 
Bacillus cereus Bacto 37 5 
nutrient 
agar 
Serratia marcescens Bacto 30 1 
nutrient 
agar 
Torulopsis rosea Bacto malt 30 1 
extract agar 
Saccharomyces cerevisiae Malt 30 1 
extract 
broth 
Aspergillus niger Sabouraud’s 30 
dextrose 
agar 


* Prior to irradiation. 
> Survived equivalent of 180,000 roentgens, the largest dose used. 
© Maximum period used was 5 seconds. 


were destroyed by doses of less than 1,400,000 
roentgens, but in all cases there were survivors 
after the application of 500,000 roentgens. From 
27 to 45 percent of the organisms were destroyed 
by 35,000 roentgens, which indicates that de- 
struction was less rapid during the initial stages 
than when pure cultures of microorganisms were 
irradiated. 

Preliminary tests carried out with raw and 
pasteurized milks, soil, and waters (Table V1) 
indicate that sterilization may be accomplished 
by doses of less than 1,000,000 roentgens. It re- 
quired less than 250,000 roentgens to destroy the 
_ organisms associated with the tap and river 
waters. Apple juice was not rendered sterile by 
the use of 1,000,000 roentgens, although over 
99.999 percent of the microorganisms were de- 
stroyed by this dosage (based on one experiment 
only). 


BIOCHEMICAL INVESTIGATIONS WITH X-RAYS 


A limited number of preliminary tests have 
been carried out to determine the effects of x-rays 


612 


produced at high voltage. 





Equivalent 
Number of roentgen value Time required 

Plating organisms required for for destruction 

medium per ml* destruction (in seconds) 
Bacto 56,000,000 1,400,000 Less than 5 
nutrient 
agar 
Bacto 381,000,000 1,000,000 Less than 5 
nutrient 
agar 
Bacto 28,000,000 1,900,000 Less than 20 
nutrient 
agar 
Bacto 480,000,000 275,000 Less than 6 
nutrient 
agar 
Bacto malt 7,800,000 610,000 Less than 29 
extract agar 
Bacto 19,700,000 b ¢ 
Wort 
agar 
Bacto 52,000 52,000 Less than 3 
Sabouraud’s 
dextrose 
agar 


produced at high voltages on fats, vitamins, and 
enzymes. 


Fats 


Butter and olive oil were subjected to roentgen 
doses varying from 0 to 1,000,000 roentgens and 
their peroxide oxygen values determined by the 
method described by Lea.*® There were slight in- 
creases in the peroxide oxygen values (ml 0.002 
thiosulphate per g of fat) for both the butter and 
olive oil. The rate of increase did not appear to be 
linear. It was observed that the orange yellow 
color of the butter was progressively destroyed 
with increased radiation. 


Vitamins 


Preliminary tests were carried out to determine 
the effect of x-rays produced at high voltage on 
ascorbic acid in pure solution and in fresh fruit 
juices. The ascorbic acid used was U.S.P. Refer- 
ence Standard Ascorbic Acid (20 mg made up to 
100 ml with 0.5 percent oxalic acid). The fruit 
juices were grapefruit and orange juices, freshly 
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expressed and strained through several layers of 
cotton gauze. Two-milliliter amounts were placed 
in each of the stainless steel dishes for irradiation. 
The x-ray doses were varied from 10,000 to 
500,000 roentgens. Reduced and total ascorbic 
acids on control and irradiated samples were 
determined by the Indophenol Colorometric 
Method.‘ 

The application of 500,000 roentgens resulted 
in rather high losses in the reduced ascorbic acid 
contents of the pure ascorbic acid solution and 
the fruit juices. However, the losses in total 
ascorbic acid were considerably less—23.4 and 18 
percent for the ascorbic acid solution and the 
orange juice, respectively. 


Enzymes 


Preliminary tests were carried out to determine 
whether x-rays produced at high voltage had the 
ability to destroy or inactivate those enxymes in 
vegetables which may be inactivated by the heat 
process known as blanching. The tests commonly 
used to check the efficacy of heat blanching by 
steam or hot water were used as an index of the 
value of x-rays as a blanching agent. Tincture 
guaiacum (1 percent in ethyl alcohol) and tincture 
benzidine (1 percent in 50 percent ethyl alcohol) 
were used to determine the presence of peroxidase 
in -inch slices of white potatoes, before and after 
irradiation. Hydrogen peroxide (1 percent) was 
used to accelerate the tests and to detect the 
presence of catalase. 

The slices of potato showed positive tests for 
peroxidase and catalase even after the application 
of 8,000,000 roentgens. 

Treatment of comminuted pieces of chicken 
liver with as high as 5,000,000 roentgens (which 
required 25 minutes) failed to destroy all of the 
catalase in the liver, although boiling for six 
minutes in water accomplished destruction of 
this enzyme. 

Preliminary qualitative tests indicated that x- 
rays produced at high voltage destroyed or 
inactivated much of the amylolytic activity of 
the enzymes in dilute (1 to 10 percent) malt 
syrup preparations when applied in doses of 
1,500,000 to 3,000,000 roentgens. 

The actions of x-rays produced at several 
million volts on crude enzyme preparations 
appear to be in agreement with those found by 
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other investigators when using x-rays produced 
at low or comparatively low voltages. Scott®® re- 
ported that enzymes are inactivated by x-rays 
only after enormous doses. Dale," as a result 
of researches with hard x-rays on crystalline 
carboxypeptidase and partly purified polyphenol- 
oxidase, showed that the percentage inactivation 
of both enzymes was a function of its concentra- 
tion for a given dose of radiation and also that the 
rate of inactivation of the enzyme preparation 
depended on its purity. Similar findings were 
made by Tytell and Kersten®® while investigating 
the effects of soft x-rays on catalase and urease. 
They demonstrated that crude preparations of 
catalase did not show inactivation in dilute or 
concentrated solution after long irradiation. 
However, they showed that crystalline prepara- 
tions were partially inactivated by large doses of 
soft x-rays. Forssberg'*!® showed that catalase 
was extremely radioresistant in vivo, whereas 
dilute water solutions of the pure enzyme proved 
to be rather sensitive. Catalase appeared to be 
highly protected from roentgen rays when in its 
natural cell environment. 


MICROBIOLOGICAL STUDIES WITH 
CATHODE RAYS 


Only a limited number of tests have been 
carried out on the uses of cathode rays produced 
at high voltages. The results of a few such tests 
are shown in Table VII. 


THE TREATMENT OF CONTINUOUSLY FLOWING 
RAW MILK WITH CATHODE RAYS 


The effect of cathode rays produced at high 
voltage on raw milk continuously passing through 
the irradiation field was determined. 

The raw milk was allowed to flow by gravity 
from a glass carboy reservoir through the irradia- 
tion field to a receiving glass carboy through glass 
tubing with rubber connectors. The carboy reser- 
voir, located on the floor of the room directly 
above the irradiation room, was equipped with 
inlet and outlet tubes. The inlet tube, which was 
plugged with cotton, entered. the top of the 
carboy through a rubber stopper. Its purpose was 
to allow sterilized air to enter the carboy to re- 
place milk withdrawn. A glass outlet tube, 
located at the bottom of the carboy, led to the 
exterior of the container through a rubber stopper 
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to a short length of rubber tubing. A screw clamp 
was attached to the rubber tubing and was used 
to regulate the flow of milk through the field of 
irradiation. 

Glass tubing, with an external diameter of six 
mm and an internal diameter of about four mm, 
was used as the means of conducting the milk 
through the field of irradiation. The area occu- 
pied by the wound portion of tubing exposed to 
irradiation (3.5 turns) was about four inches 
square. This portion of the tubing, which was to 
receive the irradiation, was coated with graphite 
to prevent damage by static effects to the glass 
during the irradiation. 

Below the irradiation unit was another short 
piece of rubber tubing which served as a connec- 
tion between the former and a piece of glass 
tubing which led through a rubber stopper to 
approximately the middle of the receiving glass 
carboy. The carboy was also provided with an 
air-inlet tube which allowed air to be displaced 
aseptically as the irradiated milk was collected. 

Before initial use, both carboys and all of the 
tubing and connections were sterilized in a steam 
autoclave. 

The raw milk was irradiated with cathode rays 
produced at 2,000,000 volts. 

Bacterial counts were made on the raw and 
cathode ray treated samples of milk, using Bacto 
nutrient agar as the plating medium and incuba- 
tion for 2 days at 30°C. There were 37,000,000 
bacteria per ml in the raw milk and two bacteria 
per ml in the irradiated milk. 

Cottage cheese was prepared from the raw 
milk and from the milk treated with cathode 
rays, using lactic acid ‘‘starters’’ obtained from a 
local source. Outside of the fact that the raw milk 
appeared to curd a little more rapidly (un- 
doubtedly associated with the large number of 
viable bacteria present before inoculation with 
the starter), there appeared to be no important 
organoleptic differences between the cheeses pro- 
duced from the raw and cathode ray treated 
milks. 


CONCLUSIONS 


1. X-rays produced at high voltage (approxi- 
mately 3 megavolts) readily destroyed bacteria, 
yeasts, and molds in massive concentrations in 
pure culture and when associated with liquid or 
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solid substrates such as milk, water, apple juice, 
ground spices, soil, and catgut sutures. 

2. Spoilage types, spore formers, non-spore 
formers, Gram-positive and Gram-negative, pig- 
mented and fluorescent bacteria were destroyed 
by high voltage x-rays. The spore formers were 
most resistant and the fluorescent bacteria the 
least resistant of the bacteria examined to date. 

3. The application of 25,000 roentgens de- 
stroyed from 94 to 99.999 percent of the non-spore 
forming bacteria. However, it occasionally re- 
quired 10 to 20 times as many roentgens to achieve 
complete destruction of the residual bacteria. 

4. Doses of 500,000 roentgens usually de- 
stroyed all non-spore formers and 99.9 percent or 
more of the spore formers. 

5. In only one case did a spore former survive 
the application of 1,500,000 roentgens. Some 
spore formers were destroyed by 1,000,000 
roentgens. 

6. The dose required to destroy bacteria with 
cathode rays appears to be similar to that neces- 
sary with x-rays. However, it was possible to 
obtain the result with cathode rays much more 
rapidly than with x-rays because the available 
output in roentgens per minute of cathode rays is 
much greater than when using the same equip- 
ment for x-rays. Preliminary tests indicate the 
possibility of sterilizing fluids by a continuous 
process. 

7. Potatoes cannot be blanched successfully 
with x-rays, although preliminary investigations 
have shown that there is some destruction of 
oxidative and amylolytic enzymes, especially on 
prolonged radiation. 

8. The peroxide oxygen values of butter and 
olive oil are increased slightly by x-radiation. 

9. Preliminary tests indicate there is a gradual 
decrease in the reduced form of ascorbic acid with 
increasing doses of x-rays, especially when a pure 
solution (in oxalic acid) of the acid is irradiated. 
Similar trends, but at a slower rate, were shown 
when orange and grapefruit juices were x- 
radiated. However, it was indicated that most of 
the ascorbic acid was oxidized only to the 
dehydroascorbic acid stage. 
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This paper presents a mathematical discussion of retarded control systems. The analysis 
presented is applicable to general control systems and is based on the solution of a linear 
differential equation of infinite order. The mathematical investigation is performed by opera- 
tional methods, and the question of stability is studied by the use of the complex variable. 


INTRODUCTION 


HE development of modern electronic ap- 
paratus has increased the range of applica- 
tion of control devices to industrial processes 
and military equipment. The mathematical 
analysis of control problems has received con- 
siderable attention in the last few years.'~* It 
has been found that the methods of electrical 
circuit synthesis and analysis are very useful in 
the design and the study of stability of linear 
control systems.* ® 
In the general literature of control systems** 
the effect of the retarded action of the controller 
is discussed. A critical analysis of the effect of 
time lag in control systems reveals the fact that 
erroneous conclusions may sometimes be reached 
if the problem is not analyzed from a rigorous 
mathematical point of view. 


It is the purpose of this discussion to analyze 
the effects of time lags in simple control systems 
by the use of the modern operational or La- 
placian transform theory that is now used so 
extensively in electrical circuit analysis and syn- 


1H. L. Hazen, ‘‘ Theory of servomechanisms,”’ J. Frank. 
Inst. 218, No. 3, 279-331, 534-580 (1934). 

*N. Minorsky, “ Directional stability of automatically 
steered bodies,” J. Am. Soc. Naval Eng. 34, 280 (1922). 

3L. A. MacColl, Fundamental Theory of Servomecha- 
nisms (D. Van Nostrand Company, Inc., New York, 1945), 
contains an excellent bibliography. 

*A.C. Hall, The Analysis and Synthesis of Linear Servo- 
mechanisms (The Technology Press M. I. T., Cambridge, 
1943), contains an excellent bibliography. 

5H. W. Bode, Network Analysis and Feedback Amplifier 
— (D. Van Nostrand Company, Inc., New York, 
1945). 

®*S. D. Mitereff, ‘Principles underlying the rational 
solution of automatic control problems,” Trans. A.S.M.E. 
57, 156 (1935). 

? Callender, Hartree, and Porter, “‘ Time lag in a control 
system,” Phil. Trans. Roy. Soc. 235A, 415 (1936). 

* N. Minorsky, “Control problems,” J. Frank. Inst. 232, 
451, 519 (1941). 
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thesis.*'® The analysis is extended to a discus- 
sion of some simple non-linear servomechanisms. 


I. FUNDAMENTAL EQUATION OF SIMPLE 
CONTROLLED SYSTEMS 


In the usual control problem, a certain quan- 
tity, ¢(t), a function of the time ¢, is a measure 
of the departure of the system from a preselected 
state which may be designated as ¢=0. The 
quantity that ¢(¢) denotes varies in various types 
of control problems. 

For example, the function ¢(¢t) may represent 
the following typical quantities: 


¢(t) =the angular deviation of a ship or an air- 
plane from its predetermined value which 
the control system is supposed to main- 
tain, or 

¢(t) =the deviation of the angular velocity of an 
engine from a certain desired angular 
velocity. . 

¢(t)=the deviation of pressure or temperature 
from a predetermined value in a pressure 
or thermal control system. 


In general, most control systems deal with 
motions of a system in the vicinity of an equi- 
librium position. The stability and general be- 
havior of such systems may be studied by the 
application of the principles of the theory of 
small oscillations. This method of analysis leads 
to linear differential equations with constant 
coefficients. The justification of such a procedure 
is based on the fact that if the small oscillation 
analysis shows that the small deviation from the 


*M. F. Gardner and J. L. Barnes, Transients in Linear 
Systems (John Wiley and Sons, Inc., New York, 1942), 
Vol. I. 

LL. A. Pipes, Applied Mathematics for Engineers and 
Physicists (McGraw-Hill Book Company, Inc., New York, 
1946). 
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predetermined value vanishes as time increases, 
then the use of this analysis is justified and 
stability of the system is assured. 

The analysis of most dynamical systems that 
contain inertia terms leads, in general, to linear 
equations of the second order with constant co- 
efficients. If the purpose of the system is to con- 
trol temperature, pressure, degree of illumina- 
tion, etc., the differential equation of the system 
is usually one of the first order. 

It is possible, of course, that the control 
system may be a “‘coupled”’ system; in this case 
its analysis leads to a system of linear differential 
equations with constant coefficients. Elimination 
between the system of equations leads to a 
linear equation with constant coefficients of 
higher order than the second. In this paper only 
systems whose equations is of the second order 
will be considered. The analysis developed is 
applicable to more complicated control systems. 

The analysis of most controlled dynamical 
systems may be based on the following equation: 


Ad(t) + Bolt) + Colt) = Fa(t) — F(t), (1.1) 


where ¢(/)=the departure of the system from 
the preselected state; A, B, C=constant coeffi- 
cients determined by the nature of the system; 
F,(t)=the external ‘“‘disturbing’’ force; F.(t) 
=the “control action” of the system. 


Il. RETARDED CONTROL ACTION 


Let us consider Eq. (1.1) in the case where the 
external disturbing force, F4(t)=0. In this case 
the equation reduces to 


Ad(t)+ Bd(t) + Co(t) = — F.(t). (2.1) 


In order for the control system to be analyzed 
it is necessary to specify the functional form of 
the “control action” of the system, F.(t). From 
a theoretical standpoint, there are an endless 


- number of possible specifications for F.(t). One 


of the simplest forms that the control action 
may be given is the following one: 


F(t) = nop(t) +mid(t) +n29(t), (2.2) 


where mo, m, and m2 are constants. It is possible 
to realize this type of control practically. How- 
ever, in any control system, whether electrical, 
mechanical, or thermal, there are always time 
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lags in the action of the system so that the con- 
trol action takes place at a retarded time. 

For simplicity, let the time lag for all the 
three types of control terms of Eq. (2.2) be T. 
Then instead of the various terms in (2.2) being 
evaluated at the instant of time, ¢, they must 
be evaluated at a past instant, t-7, so that we 
must write: 


F(t) =noo(t—T)+nio(t—T)+n.@t-—T). (2.3) 
If we now substitute F.(f) into (2.1), we obtain: 


A $(t) + Bolt) + Co(t) +noo(t— T) 
+mo(t—T)+nd(t—T)=0. (2.4) 


This is the equation of the controlled system 
in the absence of an external disturbing function. 
Equation (2.4) is a mixed difference-differential 
equation. It may be reduced to a differential 
equation by the use of the following symbolic 
form of Taylor’s expansion :'° 


e~T9(t)=o(t—T), (2.5) 


where D is the derivative operator d/dt and e is 
the base of natural logarithms. Hence we may 
also write: 

o(t—T) =De-T’9(t), (2.6) 
and 

o(t—T) = D*e-T9(t). (2.7) 


Substituting these expressions in (2.4), we 


[AD?+BD+C+me~7?+n,De—?? 
+n2D*e-T \o(t)=0. (2.8) 


This is a linear differential equation of infinite 
order with constant coefficients. In the mathe- 
matical literature it is designated as a “hystero- 
differential’ equation." This equation arises in 
the study of dynamical systems whose behavior 
depends on the past history of the system and 
also in certain mathematical investigations in 
economics." 

The characteristic equation of the differential 
equation (2.8) is given by 


(AZ?+BZ+C) 
+e-7?(n2Z?+niZ +m) =0= F(Z). (2.9) 





uL. Silberstein, ‘‘On a _ hystero-differential equation 
het) in a probability problem,”’ Phil. Mag. May 29 
(1940). 

2H. T. Davis, The Theory of Linear Operators (The 
Principia Press, 1936), Chapter VI, contains an excellent 
bibliography. 


JOURNAL OF APPLIED PHYSICS 











Stability of the system 1s assured if Eq. (2.9) 
has no roots with positive real parts. An ap- 
proximate analysis of stability may be based on 
the expansion 


eT =1—TZ/1!4+T?Z?/2!—T3Z3/3!+---. (2.10) 


If the time lag is small, then an approximate 
analysis of stability can be carried out by 
neglecting the 7? term in the expansion (2.10), 
so that Eq. (2.9) becomes 


(AZ?+BZ+C)+[(1-—(TZ/1!) ] 


X (m2oZ?+,Z + no) =0, (2.11) 
or 
— Tno.7*+Z?(n2+A —mT) 
+Z(n,+B—noT)+(C+m0)=9. (2.12) 
This cubic equation is of the form 
A,Z?+A iZ°+A 2 +A3=0. (2.13) 
The Hurwitzian criteria of stability" are 
Ay>0, Ai>0, Az>0, As>O (2.14) 


and 


A\A2>AovA3. 


It must be noticed that if more terms of the 
expansion (2.10) are taken, the Hurwitzian cri- 
teria of stability are not satisfied. This indicates 
the need of checking the approximate theory by 
exact analysis, since the results of the approxi- 
mate theory are inconclusive and hence of no 
value. 


III. OPERATIONAL ANALYSIS OF RETARDED 
CONTROL ACTION 


The physical significance of retarded control 
action may be seen if Eq. (2.4) is analyzed by 
the modern operational or Laplace transform 
method.*!® Let us consider Eq. (2.4) that in- 
cludes the effect of the retarded control action 
and also an external disturbing function of the 
“unit function”’ type which will be designated 
by the symbol 1(t) and defined by 


0 for t<0 


(= 1 for t>0° 


(3.1) 


In this case, the hystero-difference equation 


3 A. Hurwitz, “Ueber die Bedingungen unter welchen 
eine Gleichung nur Wurzulen mir negativen rellen Theilen 
besitzt,”” Math. Ann. 46, 273 (1895). 
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(2.4) may be written in the form 


Ag(t) + Bolt) + (4) +noo(t—T) 
+mo(t—T)+n2d(t—T)=1(t). (3.2) 


Following the usual method used in solving 
differential equations with constant coefficients 
by the Laplacian transformation,®!® we intro- 
duce the Laplace transform 


L(t) = 6(p). (3.3) 


The system will be supposed to be at rest at 
t=0 so that ¢(¢) and its first derivative ¢(t) are 
both zero at t=0. We then have from the theory 
of the Laplacian transformation,® !° 





Lo(t) = pe, (3.4) 
Lo(t) = p*6, (3.5) 
Lo(t—T)=e-78, (3.6) 
Ld(t—T) = pe-™"8, (3.7) 
Lo(t—T) = p’e-76, (3.8) 
L1(t)=1. (3.9) 
As a consequence of these transforms, Eq. (3.2) 
transforms into the equation 
(Ap?+Bp+C)0+e-??(nop?+nip+mo)0=1. (3.3) 
If we solve for 6, we have 
1 
6 = --—— ——- -——. —- . (3.4) 
(Ap?+Bp+C)+e-7?(nop*+nip+no) 
Now let 
Zo(p) = (Ap?+Bp+C) (3.5) 
and 
Zi(p) = (m2p?+mp+no), (3.6) 
then (3.4) may be written in the form: 
1 
ao f 
Zo(p) te 7Z1(p) , 
(3.7) 





~ Zo(p)\1+e-7”(Z/Zo)} 


619 











Fipl- Punt 


was 














Fic. 2. 


As in the usual method used in the study of 
wave propagation by the Laplacian transforma- 
tion, we expand the term in braces in the following 
manner : 








1 Z\ 
: =1l1—e Tp 
f1+e°7?(Z1/Zo)} Zo 
Zi 2 Z\ 3 
+e on( ‘) —¢ a : ) +--+. (3.8) 
Zo Zo 
Hence (3.7) becomes, 
1 Zi 
§an———_ — g- TR 
Zo(p) Zi 
ZY Zi 
+e-?7?—__ —¢-3TP__+..., (3.9) 
Ze Zo 


The series may be shown to be uniformly 
convergent. Now, the required deviation of the 
system is given by 


o(t)=L~"0(p), 


the inverse transform of @(p). The physical sig- 
nificance of the various terms of Eq. (3.9) are 


(3.10) 


evident. In the range of time between 0<t<T 


the inverse transform of the second term is zero, 
so that the first term gives the response of the 
system to the external disturbing force before 
the controller has had time to act. In the range 
‘T <t<2T the controller has acted on the sys- 
tem, and the value of ¢(¢) is obtained by taking 
the inverse transforms of the first two terms of 
(3.9) since the inverse transform of the third 
- term is zero until t=27 and the controller has 
acted again. Each succeeding transform is of 
increasing complexity. 

As a simple example of this method of solu- 
tion, consider a system with no appreciable 
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Pic.. 3. 
inertia (such as a thermal control system) so 
that A =0, and let the control action be given by 
Fo=nyp, (3.11) 


so that we have derivative control. In this case, 
the functions Z» and Z, reduce to 


Zy=Bp+C=B[p+(C/B) | (3.12) 
and 
Z\ = mp. (3.13) 
To simplify the writing, let 
C/B=a. (3.14) 


For this case, the general series (3.9) reduces to 





1 np 
0(p) =——_—_—_ — e~T?______ 
B(p+a) B*(p+a): 
2p2 3A3 
pg tre vt? —e-3Tp ‘ae (3.15) 
B*(p+a) B*(p+a)* 


¢(t) may now be calculated by calculating the 
individual transforms of (3.15). We thus obtain 


o(t)=1/C(1—e-*") O<t<T, (3.16) 
o(t) =1/C(1—e~") 
—n,/BXt—T)e"-?) T<t<2T, (3.17) 
(t) =1/C(1—e-*") —m,/B2(t— T)e-#—) 
+-n;?/2B%{2(t—2T)e-*t27 
—(t—2T)'ae*""27)} IT <t<3T. (3.18) 


If the time lag, 7, is large, the method of solu- 


tion based on Eq. (3.9), although laborious, may 


be carried out as many terms as desired. In this 
manner the behavior of the system may be 
studied for as long a value of time as is required. 

The usual approximate theory valid for small 
time lags may be based on Eq. (3.7). The func- 
tion e~7? may be expanded in a power series in 
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the form 


e-T? =1—Tp/1!4- T2p?/2!—T8p3/3!4---, 


(3.8) 


and as many terms retained as is desired. The 
characteristic equation of the system is 


F(p) =Zo(p) +e-7?Zi(p) =0. 


If its roots are known, then the inverse trans- 
form of 6(p) may be completed by the theory of 
residues or an application of the Heaviside ex- 
pansion theorem. This procedure is, however, 
most laborious and does not give a closed form 
solution. 


(3.9) 


IV. TIME LAG IN A SIMPLE SERVOMECHANISM 


As a simple application of the general theory 
discussed above, consider the system of Fig. 1. 
Let the following notation be introduced: 


6,(p) =the Laplace transform of the input elec- 
tromotive force impressed on the am- 
plifier. 

6.(p) =the Laplace transform of the angular 
displacement of the motor shaft. 
64(p)=the Laplace transform of the output 

potential of the output of the amplifier. 
6,(p)=the Laplace transform of the output 
potential of the potentiometer. 


In order not to obscure the analysis, and to 
illustrate the salient features of time lag, the 
amplifier will be supposed ideal, the inertia of 
the rotating parts of the motor will be assumed 
negligible, and the potentiometer action will be 
supposed to be ideal but will be assumed to 
introduce a time lag of magnitude 7. As a 
consequence of the above assumptions, we will 
suppose that the component parts of the system 
will have the following transfer functions: 


P,=yu=the transfer function of the 
amplifier. 

Py=G/Rp=the transfer function of 
the motor, where G is a positive 
constant and R is a positive constant 
that is a measure of the friction of 
the rotating parts of the motor. 

Q=kye-T? = the transfer function of the 
potentiometer feed-back loop where 
ko is a positive constant and T is a 
measure of the time lag of the feed- 
back loop. 


(4.1) 


(4.2) 


(4.3) 
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Hence the transfer function of the forward 
loop is given by 
P=P4Py=uG/Rp. (4.4) 


The relation between 49 and 6; is given by the 
usual equation 





69=[P/(1—PQ) 6, (4.5) 
and in this case this becomes 

uG uG 1 
6) = —__—__——_ = —_—————.. (4.6) 

Rp+puGkoe7? Rp uG 

|! +—he | 
Rp 
If we let 

b=ynG/R, (4.7) 


Eq. (4.6) may be expanded in the following series 
form: 


b b3 
09 =———kye- 7? +—k?e 77” 
2 p® 
5! 
——hkye*T?+4---, (4.8) 
p' 


The angular displacement of the motor shaft, 
6o(t), is the inverse transform of 4» and is given by 


go(t)=bt for O<t<T, 
go(t) = bt—(b*ko/2)(t-T)? for T<t<2T, 
do(t) = bt — b*ko/2(t—T)? 

+ (b*k?/6)(t—T)3 


for 27 <t<3T. 


(4.9) 


As many terms may be computed as desired. 
This method gives a closed form solution for 
any value of time. An alternative method of de- 
termining the inverse transform of 0 9 depends 
on the determination of the roots of the char- 
acteristic equations of the system." This equa- 
tion may be written in the following form: 


F(p) = p—ae-T?=0, (4.10) 
where 
a=pGko/R. (4.10a) 
If 7 =0, then p= —a is solution of the system 
and @¢o(t) is given by 
go(t) =1/ko(1—e-**). (4.11) 


In this case the angular displacement of the 
motor follows the unit input signal rapidly, and 


14 N. Minorsky, “Self excited oscillations in dynamical 
systems possessing retarded action,”’ Trans. A.S.M.E. 64, 
A65-A71 (1942). 
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the system is stable. An approximate analysis 
that given qualitative results for small values of 
the time lag, 7, is based on approximating the 
exponential function e~7” by the first three 
terms of its Maclaurin series expansion in the 
form: 


F(p) =p+a[1—Tp+(T%p?/2)]=0. (4.12) 


In this case the characteristic equation is 
p?+[2(1—aT)p/T*a}+(2/T*)=0. (4.13) 


To this degree of approximation it may be 
seen that if 
T>1/a (4.14) 
the characteristic equation will have roots whose 
real parts are positive and the system will be 
unstable in that the effect of transient dis- 
turbances will be magnified indefinitely. If 


T=1/a, (4.15) 


then 

p=j(V2/T) 
and the output will oscillate with an angular 
frequency of w=v2/T. A precise analysis re- 
quires an investigation of the roots of Eq. (4.10) 
by assuming a solution of the form 


(4.16) 


p=atjuw, (4.17) 


where w and a are real and are to be determined 
from Eq. (4.10). Substituting (4.17) into (4.10), 
we obtain 


tanwl’=w/a, e??*=a?/(a*—w*). (4.18) 


If we eliminate a between these two equations, 


we obtain 
2Tw 1 ww 
exo| - -~|-— +a’. 


tanw7 J} a*® tan’wT 





(4.19) 


This transcendental equation may be solved 
‘graphically for the possible values of w if nu- 
merical values of T and a are given. 

Having determined the values of w, then the 
corresponding values of a may be determined by 


- means of the equation 
a=—w/(tanw7). (4.20) 


If the values of a and w, and hence the zeros 


of F(p) given by Eq. (4.10), have been deter- 
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mined, then the response of the system, ¢o(?), 
is given by®*!° 


go(t) = residues of (b/pF(p))e?' (4.21) 


at the zeros of F(p). 

This procedure is extremely laborious. The 
wave solution of Eq. (4.9) is much simpler and 
has greater physical significance. 


V. STABILITY ANALYSIS 


If only a knowledge of the stability of the 
system is desired, then it is only necessary to 
investigate whether or not the function F(p) has 
zeros in the positive half of the p-plane. This 
may be done by an application of the following 
well-known theorem of Cauchy. 

“Given a function, F(p), that is analytic 
except for pole singularities within a closed con- 
tour in the complex p-plane. If the point p 
traverses the given contour once in the clock- 
wise direction, then the number of times that 
the representative point in the F(p)-plane en- 
circles the origin in the clockwise sense is equal 
to (No—N,), where No is the number of zeros 
and N, is the number of poles of the function 
F(p) that are within the given contour.” 

A suitable contour that may be used to de- 
termine whether or not the function F(p) has 
zeros in the positive half of the complex p-plane 
is shown in Fig. 2. 

The radius of the circle R tends to infinity. 
In the case under consideration, F(p) = p+ae-7?. 
Now along the axis of imaginaries, let p=jw, and 
consider 


F(jw) =jw+ae-T?*, (5.1) 


As w runs from — « to + ~ the point p traverses 
the imaginary axis. The point 0 of the p-plane 
corresponds to the point a of the F(p) plane; 
as w increases, the plot of ae~7’* in the F(p) 
plane is a circle of radius a as shown in Fig. 3. 

F(jw) may be regarded as the addition of the 
vector jw and the vector ae~7’*. As long as p 
continues to increase along the axis of imagi- 
naries, the vector ae~7’* revolves around the 
circle of Fig. 3. The plot of F(jw) may be con- 
structed graphically. Since limit | F(p)|—>~, the 
portion of the F(p) plot corresponding to the 
infinite semicircle is of no concern. It is only 
necessary to plot F(jw) for values of 0<w<~, 
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since the plot for — « <w<0 may be obtained 
as the mirror image of the plot for positive fre- 
quencies. It will be seen that whether or not 
F(jw) encloses the origin will depend on the 
relative magnitudes of a and T. If a and T are 
small, F(jw) will not enclose the origin and the 
system will be stable. If T is large and a>1, then 
the origin will be enclosed and the system will 
be unstable. This procedure may be applied to 
more complicated systems, and it is essentially 
an application of the familiar Nyquist criterion® 
to systems having time lag and whose charac- 
teristic equation is of the transcendental type. 


In conclusion, it may be mentioned that in 
many control mechanisms certain component 
parts of the system are non-linear. However, in 
order to facilitate the analysis the effects of the 
non-linear components is regarded as introduc- 
ing time lags. The more rigorous approach to 
problems of this type is to consider the actual 
non-linear differential equations of the system. 
This procedure has been carried out by Bugakov 
in a recent paper.!® 


16 B. V. Bugakov, ‘‘ The method of van der Pol and its 
application to non-linear control problems,” J. Frank. 
Inst. 240, 1 (1946). 








Erratum: Response of Linear Resonant Systems to Excitation of a 
Frequency Varying Linearly with Time 
[J. App. Phys. 19, 242 (1948) ] 


GUNNAR HOK 
Wesleyan University, Middletown, Connecticut* 


HE author is indebted to Mr. R. R. Bennett 

of the California Institute of Technology 

for pointing out .some errors of sign in Eq. (39) 

in the Appendix to the above paper. Since the 

computations for the graphs provided a satisfac- 

tory check of the consistency of the numerical 

methods used, the errors have evidently crept in 

during the final arrangement of the subject matter 
for publication. 

The correct form of the equation is: 


1 1 3 3°5 





Prat — + $e: 
2y 47° 8y* 167’ 327° 
1 1 3 
>— cosé+-—- sin3@—-—— cos5@—-:- 
2r 4r’ 8r° 
1 1 a 
+7} —— sind+— cos36+— sin5é—--- (39) 
2r 4r* 8r® 


Furthermore, Eq. (35) should be replaced by 





* Now at University of Michigan, Ann Arbor, Michigan. 
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the following expansion, also obtained by inte- 
gration by parts: 


t 
e* f e*dT 


—2 


1 Gf") Gf"  5f*®O-3L/"OF 
+ == ~ 





3 4 5 








v1 1 $1 1 
If) — Lf" OF" 
- - 
0 TASH OF" O— LOCH") P— 1S HF" P 


¢1" 


—-++ (35A) 


where ¢,; =d¢/dt. The asymptotic quasi steady- 
state solution for slow frequency modulation is 
obtained when all subsequent terms are negligible 
compared to the first one. (Compare Eq. (23).) 

If f(t) = et?, #1; =wo and both members are multi- 
plied by je’, Eq. (35A) gives Eq. (39). 
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Rising Sun Magnetrons with Large Numbers of Cavities* 


A. V. HOLLENBERG,** N. KRoLL, AND S. MILLMAN** 
Columbia Radiation Laboratory, Columbia University, New York, New York 


(Received February 2, 1948) 


Studies of the rising sun design for magnetron oscillators have been extended to large numbers 
of cavity resonators. Successful operation of pulsed magnetrons having as many as 38 cavity 
resonators has been achieved in the 1.25-cm wave-length region. The open cavity resonators 
heretofore used in rising sun magnetrons were found unsuited for oscillators containing a large 
number of cavities. The extension to higher numbers was made possible by the use of closed-end 
resonators. This modification in cavity geometry produces a favorable redistribution of the 
spectrum of the resonant modes of the rising sun structure, thereby making possible the opera- 
tion of the magnetron in the desired x-mode without disturbance from other modes. The in- 
crease in the number of cavity resonators for magnetrons should serve to facilitate the extension 
of present magnetron techniques to the millimeter range and to increase the radiofrequency 
power obtainable from magnetrons at 1 cm and at longer wave-lengths. 


I. INTRODUCTION 

HE rising sun magnetron’? was developed 

in the course of a search for stable mag- 
netron operation at wave-lengths close to 1 cm. 
With this design, stable operation is achieved in 
the desirable x-mode without the use of straps. 
The resulting simplification in anode block con- 
struction, together with some electrical ad- 
vantages, makes this type of magnetron particu- 
larly suitable for the generation of radiation in 
the short wave-length end of the microwave 
region. The early work'* was confined almost 
entirely to the development of magnetrons having 
18 cavity resonators in the anode block, a cross 
section of which is shown in Fig. 1. This number 
was quite adequate for the design of pulsed mag- 
netrons in the neighborhood of 1.25 cm. With the 
development of interest in millimeter generators, 
however, it was realized that a larger number of 
resonators would be required. This is a direct 
consequence of the scaling laws in magnetrons* 
which require that if a scaled copy be made of a 
known workable magnetron by reducing all linear 
dimensions by a factor K without changing the 

’ number of resonators, the wave-length decreases 
by the same factor while the anode voltage and 


* This paper is based on work done for the OSRD under 
Contract OE Msr-485. 

** Now at Bell Telephone Laboratories, New York, 
New York. 

1S. Millman and A. T. Nordsieck, J. App. Phys. 19, 
156 (1948). 

?N.Krolland W. E. Lamb, Jr., J. App. Phys. 19, 166 (1948). 

3 NDRC Div. 14, Report No. 588. 

‘J. B. Fisk, H. D. Hagstrum, and P. L. Hartman, Bell 
Sys. Tech. J. 25, 202 (1946). 
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anode current remain unchanged despite the re- 
duction in cathode area by a factor K*. This im- 
poses severe limitations on the cathode both with 
regard to electron emission and to heat dissipa- 
tion. A similar need for a larger number of 
resonators arises when one wants to design a 
short wave-length CW magnetron from a known 
workable pulsed magnetron at the same wave- 
length. It was therefore considered desirable to 
attempt to extend the rising sun design to anode 
blocks having a number of resonators appreciably 
greater than 18. 

In the following pages we shall first describe 
(Section I1) some of the difficulties encountered 
in attempts to extend the rising sun design to 
anode blocks having a large number of cavity 
resonators. This will be followed (Section II1) by 
a discussion of the use of closed-end cavity 
resonators, together with an explanation of how 
this modification in the anode block design helps 
to remove some of the difficulties encountered 
with *magnetrons using open-end resonators. In 
Section IV some features incident to the con- 
struction of magnetrons having as many as 38 
cavities will be discussed. In Section V a com- 
parison will be given of the performance of the 
various tube types studied in these experiments. 
In the concluding section a general discussion of 
the results will be presented. 


Il. OPEN CAVITY RESONATORS 


The early work consisted principally of at- 
tempts to develop rising sun magnetrons, for 
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pulsed operation in the 1.25-cm wave-length 
region, having anode blocks of the same general 
geometry as the satisfactory 18-vane magnetron, 
but with a larger number of cavity resonators. 
While no difficulty was experienced in con- 
structing a successful oscillator having 22 cavity 
resonators,’ considerable effort was expended in 
developing a moderately satisfactory 26-cavity 
magnetron for operation in about the same volt- 
age range (~15 kv) as the 18-cavity magnetron. 
lt was found that the performance of the 26-vane 
tubes was too readily affected by small changes 
in the value of 7;, the ratio of the resonator 
depths of the two different cavities (Fig. 1), and 
in the value of o, the ratio of the cathode diame- 
ter (D.) to the anode diameter (D,). For the low 
values of r;, mode competition was observed from 
the m = 12 resonance, which is next to the 7-mode 
(n=13) on the short wave-length side, with a 
mode separation of a few percent. For high values 
of r;, the competing resonances were those mem- 
bers of the long wave-length multiplet! of the 
rising sun spectrum for which 5 3, and opera- 
tion in these unwanted modes occurred at nearly 
the same voltage as that for the 2-mode. 

Attempts to design an open cavity magnetron 
having 38 resonators have met with complete 
failure. No value of 7; could be found which 
would eliminate the interference from both the 
n=18 mode and from some of the long wave- 
length resonances. The interference was so severe 
as to prevent any effective operation in the 
x-mode. 

Similar mode competition was subsequently 
observed in 18-cavity magnetrons, except in 
much milder form. The performance of these was 
carefully studied*® in the range of r; values from 
1.35 to 2.40. It was found that for values of r; 
below 1.6 interference from the n=8 mode be- 
comes troublesome, while for values above 2.0 
competition is observed from members of the 
upper multiplet. In the range of 7; values of 1.6 to 
2.0, entirely satisfactory 7-mode operation can be 
obtained. 

From the results of the experiments with open- 
ended cavity magnetrons, the following con- 
clusions were arrived at: 

(a) When r; is too small, competition from the n =((N/2) 


—1] mode interferes with stable x-mode operation of the 
magnetrons. 
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(b) When r; is too large, mode competition from mem- 
bers of the long wave-length multiplet having n=3 pre- 
vents satisfactory operation. 

(c) The range of satisfactory values of r; decreases as the 
number of resonators is increased. For magnetrons in the 
voltage range under discussion, satisfactory open-cavity 
vane type magnetrons with N>26 can probably not be 
designed. 


In order to understand the nature of the mode 
competition we review the general features of a 
typical mode spectrum for the rising sun resonant 
system. These have been fully discussed in pre- 
ceding papers.'? In Fig. 2 is reproduced a plot of 
the computed values of the resonant wave- 
lengths of the various modes as a function of 7; 
for N=18. The long wave-length multiplet in- 
cludes the resonances with mode numbers 1 to 4. 
The mode spectrum for a resonant system with a 
larger number of cavities is similar to the one 
shown. The long wave-length multiplet contains 
additional resonances close to the N=4 mode 
with mode numbers up to V/4 or (N—2)/4, de- 
pending upon whether NV/2 is even or odd. On the 
short wave-length side of the #-mode (n= N/2) 
there are the resonances with mode numbers 
(N/2—1), (N/2—2), ---N/4 or (N—2)/4 in or- 
der of decreasing wave-length. From the known 
principles of magnetron electronics the fact that 
the operating voltages for the competing reso- 
nances nearly coincided with that for the 7-mode 



































ri =(Di—Da)/(Ds —Da) o =D./Da 








Fic. 1. Rising sun anode block and cathode. 
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Fic. 2. Resonant wave-lengths of rising sun anode modes 
as a function of the ratio r;(.N = 18). 


implied that 
N 


| Y | Ana—As, 
2 


where \, is the wave-length of the competing 
resonance of mode number n, and 7 is the effective 
mode number,*** defined by the relation 


y=n+mN for a symmetrical anode block, 
and by 
y=n+m(N/2) for a rising sun resonant system, 


with m=0, +1+2:---. 

From the measured values of \, and \, and the 
above relation it was apparent that the effective 
‘ mode numbers for the competing long wave- 
length resonances did not correspond to y =” the 
fundamental, but to y=n—(N/2), one of the 
Hartree harmonics‘ of the mode number n. For 
the competing mode n=[(N/2)—1], however, 
the closeness of \, to \, makes it difficult to 
decide between y=n and y=(N—n) for the 


*** ~ as used here has the same meaning as the symbol 
k in reference 4. 
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effective mode number of the competing reso- 
nance. 

The general behavior of the spectrum as a 
function of r; leads one to conclude that competi- 
tion from the nm =[(N/2)—1] mode occurs when 
the separation between this mode and the 
m-mode [n=(N/2)] becomes too small, while 
competition from members of the long wave- 
length multiplet occurs when its separation from 
the x-mode becomes too large. From the meas- 
ured wave-lengths of competing modes it has 


‘been possible to show that the minimum separa- 


tion permissible between the w and [(N/2)—1 | 
modes does not change appreciably with N. On 
the other hand, the maximum allowable separa- 
tion between the z-mode and the long wave- 
length multiplet decreases considerably as N 
increases. It is this latter circumstance which re- 
duces the available range of r; as N is increased. 

The decrease in maximum allowable separation 
between the z-mode and upper multiplet modes 
is readily understandable if one takes as a condi- 
tion for mode free operation that the starting 
voltages for the members of the upper multiplet 
with m > 3 operating in the [ (V/2) —n ] harmonics 
exceed the starting voltage of the -mode. This 
leads to the conditionf that A,{(N/2)—3] 
<(N/2)dm or A3<Ae[N/(N—6) ]. 

The above observations suggest that a modi- 
fication of the resonant system which would 
depress the long wave-length multiplet relative to 
the z-mode might eliminate mode competition 
from the upper multiplet. Of course such a 
modification would elevate the starting voltages 
of the harmonics of the = 1 and m = 2 modes into 
the neighborhood of the z-mode. Since, however, 
the open-end 18-resonator magnetrons had oper- 
ated satisfactorily with such a distribution of 
starting voltages, one might hope that this would 
not be a source of trouble. 


Ill. CLOSED-ENDED RESONATORS 


One method of depressing the long wave-length 
multiplet consists of brazing cover plates to the 
end surfaces of the anode block, as shown in 
Fig. 3. The closing of the resonators depresses the 


t This condition is somewhat stronger than necessary as 
the »=3 mode omes a progressively weaker com- 
petitor as N is increased. Some reasons for this will be 
given later. 
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wave-length of each mode in accordance with the 
approximate formula 


1 1 


1 
eee = + tT 
* closed open (2h)? 





where h is the anode block height. The resultant 
depression of the long wave-length modes relative 
to the rest of the spectrum is exhibited in Fig. 4 
for the case N = 18, from which it is clear that for 
a given separation between the z-mode and the 
n=8 mode the separation between the 2-mode 
and the upper multiplet resonances is con- 
siderably less for the closed-end anode blocks. 

To test this idea a magnetron having 26 closed- 
end cavities was constructed. The anode and 
cathode diameters were approximately the same 
as the corresponding diameters in the open cavity 
26-vane tube. The height of the anode block had 
to be increased to exceed one half-wave-length. 
The resonator depths were also increased so as to 
compensate for the wave-length shortening 
effects of closing the ends. It was observed that 
competition from the upper multiplet modes 
could be entirely eliminated, even in tubes with 
an r; value as large as 3.3. These encouraging re- 
sults were followed up by the construction of 38- 
cavity closed-end magnetrons, where r; values as 
large as 2.7 were used, without any appreciable 
competition from the long wave-length resonant 
modes. The value of r; below which competition 
from the shorter wave-length (NV /2) —1 mode be- 
gins to appear is greater for the closed-end anode 
blocks. Since, however, the upper multiplet 
moding problem is eliminated, there remains no 
objection to large values of r;. 

No competition has been observed from the 
harmonics of the m= 1 and n=2 resonances even 
though the operating voltages for these harmonics 
are considerably lower than that for the 7-mode 
and therefore are passed through in every appli- 
cation of the voltage pulse to the magnetron. 
Furthermore, the competition from the n=3 
mode has been rather light, particularly for mag- 
netrons with the larger number of resonators in 
which the operating voltage for the [(N/2)—3] 
harmonic has been quite close to that of the 





it This expression is based on the assumption that the 
open cavity wave-length corresponds to an infinite anode 
height and that the brazed-on covers include the inter- 
action space region. 
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n-mode. These observations imply that conditions 
are unfavorable for the starting of harmonics of 
low number modes even though the voltage 
satisfies the resonance condition for operation. 

A partial explanation for the absence of inter- 
ference from harmonics of low number modes 
may be derived from a study of the dependence 
of the tangential and radial components of the r-f 
electric field on the radial coordinate for the vari- 
ous Hartree harmonics. Approximate expressions 
for these are? 


‘s (—)" ' siny§ (1—(De, ar)*i'} 
~ yO {1—(D./D,)2'"\}" 





. ty vi- ' siny@ {1+(D./2r)?!7'} 
" Xp, v8 {1—(D./D,)2"}" 


where @ is half the angle subtended at the origin 
by the mouth of the resonant cavity. These ap- 
proximations improve with increasing values 
of |v}. 

It may be seen from the above expressions that 
the field amplitude falls off rapidly as r decreases 
from D,/2 to D./2, with the rate of decrease 
increasing with |y|. Thus near the cathode sur- 
face the low |y| components, and in particular 
the fundamental, generally predominate. In con- 
sidering the possibility of a given mode competing 
with the z-mode it seems reasonable to suppose 
that the chance of x-mode operation improves as 
|y| increases. For as |y| increases the field 
amplitude near the cathode becomes weaker and 
it becomes more difficult for the oscillation in the 
competing mode to build up. In the case that the 





Fic. 3. Rising sun anode block with cover plates. 
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competing mode is being driven by a harmonic 
one must further consider the effect of the funda- 
mental component, which in general will be much 
stronger than the harmonic in question near the 
cathode. The disparity in strength will increase 
with increasing |y|—m, which in the case 

y| =[(N/2)—n] is [(N/2)—2n]. The main 
effect of the presence of the fundamental proba- 
bly consists of distorting the electron orbits and 
thus preventing proper coupling with the har- 
monic and thereby improving the probability of 
starting in the r-mode. 

According to the above considerations one 
would expect the lower m members of the upper 
multiplet to be the weaker competitors, a con- 
clusion which is consistent with the observed 
non-interference from the (n=1) and (n=2) 
modes. ‘The fact that the »=3 component be- 
comes a weaker competitor as N increases is con- 
sistent with the facts that (a) the disparity 
between the [(N/2)—3] and N/2 components 
decreases as N is increased due to the increased 
value of D./D, which always accompanies an 
increase in N, and (b) the disparity between 
the y=3 (fundamental) component and the 
\y| =(N/2)—3 increases as N increases, which 
increase is only partially compensated by the in- 
crease in D,/D,. 

The above considerations can also be applied 
directly to the question of the effect of the ratio 
D./D, wpon moding. It has been quite generally 
observed that increases in D./D, always reduce 
competition from the long wave-length multiplet. 
This would follow from the fact that any such 
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Fic. 4. Ratio of wave-length of the lowest member of the 
upper multiplet to that of the #-mode as a function of 
the ratio of the wave-length of the x-mode to that of the 
highest member of the lower multiplet, for open- and 
closed-end anodes (N = 18). 
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increase increases the N/2 component amplitude 
relative to that of the competing |y| component. 

We now examine the effect produced by the 
closing of the resonator ends on the zero com- 
ponent contamination of the 7-mode interaction 
field.? Within fairly wide limits this contamination 
has very little effect on the moding situation, but 
it affects the operating efficiency of the magnetron 
for some region of magnetic field. It is worth 
noticing, however, that the situation is no worse 
for the closed-end resonant system than for the 
open-end cavity magnetrons. This is illustrated in 
Fig. 5 in which the ratio of zero to N/2 com- 
ponent amplitudes at the anode surface is plotted 
as a function of \,/As for the same two anode 
blocks used in Fig. 4. 


IV. CONSTRUCTION 


In this section we shall review briefly some of 
the features of construction of the magnetrons 
discussed above, emphasizing some of the prob- 
lems that arose in connection with an increased 
number of resonators. The study of the variation 
of the performance with the ratio of resonator 
depths required the construction of one or more 
magnetrons for each value of 7;. The effect of 
varying the ratio of cathode-to-anode diameters 
could frequently be accomplished by replacing 
cathodes of different diameters in the same tube. 
Counting each cathode replacement as leading to 
a new magnetron, over 90 magnetrons were built 
in the course of these experiments. 

The construction of these magnetrons was 
similar to that of 18-vane, 1.25-cm rising sun 
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Fic. 5. Ratio of zero mode to x-mode amplitude at the 
anode surface as a function of the ratio of the wave-length 
of the r-mode to that of the highest member of the lower 
multiplet for open- and closed-end anodes (N = 18). 
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magnetrons which were, during the period of 
these experiments, in an advanced state of de- 
velopment in this laboratory and elsewhere.* 4 
While some of the tubes were made with radially 
mounted cathodes, more attention was paid to 
the tube design in which the cathode is mounted 
axially and pole pieces are built into the vacuum 
envelope. The latter design yields a more efficient 
magnetic circuit and greater accuracy in centering 
the cathode with respect to the anode. Both 
features are particularly desirable for short wave- 
lengths. For these reasons it has seemed desirable 
to meet the problems which might arise in this 
design when the number of resonators is in- 
creased. 

The essential features of a typical magnetron 
with inserted pole pieces and axially mounted 
cathode are shown in Fig. 6. The wave-guide 
output shown was used in all cases. Some 
variation in the dimensions of the output quarter- 
wave transformer was required in the various 
tube types. 

In some of the closed-end anode blocks, steps 
were introduced at the ends of the cylindrical 
anode surface in order to provide adequate volt- 
age clearance between the anode and the end hat 
of the cathode in scaled down versions of these 
designs. A sectional view of a portion of the 
cathode and anode block involved in this modifi- 
cation is shown in Fig. 7. Successful operation 
both with and without the step has been ob- 
served. The built-in pole pieces had cylindrical 
holes in the center, large enough to pass the end 
hats of the cathode. This design feature is used to 
facilitate the alignment of cathodes and their 
replacement. Since the end hats were only slightly 
smaller than the anode diameter, magnetrons 
with a large number of cavity resonators required 
large holes in the pole pieces with some resulting 
complications. 

One consequence of the large holes in the pole 
pieces was an unfavorable distribution of mag- 
netic field in the interaction space near the pole 
pieces, giving rise to excessive current leakage in 
some of the 38-vane magnetrons containing the 
large pole piece holes. Permendurfff cathode 
hats, shaped to produce uniform field throughout 
the interaction space, according to electrolytic 


ttt Reference 4, p. 328. 
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Fic. 6. Typical magnetron assembly. 


bath model tests, were employed in a later design, 
and the current leakage was materially reduced, 
to an amount not appreciably more than is ob- 
served with magnetrons having no holes in the 
pole pieces. 

Another consequence of the large pole piece 
hole was serious leakage of power out of the 
magnetron through the coaxial transmission line 
formed by the cathode mounting stem in this 
hole. A number of facts concerning the leakage 
have been observed, and will be listed in the next 
paragraph. This will be followed by the interpre- 
tation of the observed behavior and by sugges- 
tions for the elimination of the leakage without 
introducing any further complication in mode 
competition. 

The most obvious manifestation of the leakage 
was the lighting up of a neon bulb placed near the 
hole in either pole piece, with an intensity in some 
cases comparable to that obtained by placing the 
bulb directly in front of the wave-guide output. A 
second manifestation was that abnormally low 
values for the unloaded Q of a leaky magnetron 
were obtained in ‘‘cold’’ measurements in which 
power was fed into the magnetron from an ex- 
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Fic. 7. Step modification at ends of anode block. 


ternal source. Moreover, bringing metallic ob- 
jects into the neighborhood of the pole piece holes 
affected the measured values of the Q’s. Power 
leakage of this magnitude was not observed from 
the 18-vane magnetrons. It varied in amount 
from tube to tube in the 26- and 38-vane mag- 
netrons. Very bad leakage was observed in one 
38-vane tube in which the [(V/2) —1_] mode was 
nearly coincident with the z-mode, but was not 
observed in the same tube when a cathode of 
slightly smaller diameter was used. In the 
latter case the two modes were more definitely 
separated. 

A comparison of pole piece hole and cathode 
stem dimensions for the 18-vane tube with those 
of the 26- and 38-vane tubes shows that only 
principal coaxial mode propagation is possible 
with the first set of dimensions, whereas second 
coaxial mode propagation as well is permitted by 
the dimensions in the 26- and 38-vane tubes. This 
indication that the leakage occurs because of 
second coaxial mode propagation appears to have 
been confirmed by the fact that it was eliminated, 
on the cathode support end, by the use of a choke 
structure (Fig. 8) whose dimensions were based 
"on the assumption of second coaxial mode propa- 
gation. Furthermore, the leakage was eliminated 
on the opposite end by a reduction, below the 
cut-off value for the second coaxial mode, of the 
diameters of the pole piece hole and of the heat 
conducting stem which protruded from the 
cathode. The diameter of the hole in the iron was 
the same in both pole pieces. The reduction at 
one end was made by the insertion of non- 
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magnetic material in order to preserve symmetry 
in the magnetic field. 

While the evidence is strong that the leakage 
was due to second coaxial mode propagation, the 
fact that it did not invariably occur when the 
dimensions permitted requires further explana- 
tion. This may be obtained by comparing the 
radiofrequency field patterns in the interaction 
space, for the z- and [(N/2)—1] modes with 
those of the first and second coaxial modes in the 
cathode support line. For the principal coaxial 
mode, no strong coupling to any of the modes of 
the resonant system is apparent, even when the 
possibility of asymmetries is taken into account. 
For the second coaxial mode, there is obviously a 
strong coupling between the cathode support line 
and the m = 1 component of the [(V/2) — 1] mode, 
but there should be no coupling to this trans- 
mission line either by the z-mode or by the zero 
component contamination present in it, under 
symmetrical conditions. Some coupling of the 
second coaxial mode to the -mode may, how- 
ever, be expected under asymmetrical conditions, 
the coupling being stronger the smaller the sepa- 
ration between the z- and [(N/2)—1] modes. 
The amount of leakage present may therefore be 
expected to be related to accidental asymmetry 
in cathode location and on the degree of excita- 
tion of the [(N/2)—1] mode. Both of these may 
vary from tube to tube. 

If the coupling to the [(V/2)—1] mode is 
much stronger than to the 7-mode, as suggested 
above, the insertion of a choke in the cathode 
support line has, for the [(N/2)—1] mode, the 
effect of connecting the interaction space to a 
short circuited section of transmission line, in 
which case the frequency of this mode alone 
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should be appreciably affected by the position of 
the choke. Such a selective dependence of the 
[(N/2)—1] mode frequency on choke position 
has been observed in experiments with a movable 
choke. These results have been reported in detail 
by one of us.® It will be necessary, therefore, to 
use some care in the location of a choke intro- 
duced for the purpose of eliminating leakage, in 
order to avoid an unnecessary reduction in mode 
separation. 

Another consequence of the increase in the 
number of cavity resonators is the introduction 
of new problems in anode block construction. 
The. technique of hobbing* had been very suc- 
cessful in forming 18-cavity anode blocks and the 
construction of blocks with the larger number of 
resonators by the same method was attempted. 
The vane thickness and the smallest space be- 
tween vanes was nearly the same as at 18 vanes. 
The latter dimension is of significance because it 
is the smallest dimension of the hob sector which 
forms the resonant cavity in the anode block, 
and is therefore important in determining the 
strength of the hob. The important differences 
between the 18 and larger N anode blocks are 
that the large anode diameter requires the dis- 
placement of a much greater volume of copper in 
the center, and that in the case of the closed-end 
anode blocks, larger cavities are required. 26-vane 
blocks were hobbed with some success, but at 38 
vanes, hob failure occurred on the second or 
third, if not on the first, use of the hob. This 
circumstance led to'the use of the cavity struc- 
ture shown in Fig. 9. The sector structure could 





Fic. 9. Hole and vane construction of anode block. 


* Columbia Radiation Laboratory Progress Report, Oc- 
tober 1, 1946. 


VOLUME 19, JULY, 1948 


“eo, 0 e- 


et oe ry 
24 30 36 42 a8 





Fic. 10. Ratios of cathode diameter, current density, 
and By for N resonators to corresponding values for 18 
resonators. Arrows indicate N values studied. 


be produced by hobbing because of the absence 
of the long sectors of large area. The holes joined 
onto alternate sectors were bored after the 
hobbing operation. Recently it has been found 
possible to produce the complete anode block by 
the hot hobbing technique.* In hot hobbing the 
maintenance of the copper in an annealed condi- 
tion greatly reduces the forces which tend to 
break the hob. 


V. DESIGN AND PERFORMANCE 


In this section we summarize the essential 
design principles and the bases for comparison 
used in evaluating the performance of the various 
tube types. Performance charts and some tube 
dimensions are presented for the more successful 
magnetrons. 

In designing magnetrons of N greater than 18, 
in predicting their advantages, and in evaluating 
their performance, a scaling procedure taking 
into account the change of N and the accom- 
panying change in « was employed, using the 
extensively studied 18-vane design as a starting 
point. The system of reduced variables formu- 
lated by Slater’ was used for this purpose. The 
coordinates V» and By of the point of tangency of 
the Hartree line for the operating mode to the 
cut-off parabola, and the anode current Jo for 
this point, form the basis of the system of reduced 
variables. Performance charts for two properly 
operating magnetrons are essentially the same, 
neglecting differences in circuit losses, when 
plotted in terms of V/Vo, B/Bo, and I/Io, even 
though they differ in wave-length, number of 


‘Columbia Radiation Laboratory Progress Report, 
January 1, 1947. 

7G. B. Collins, Microwave Magnetrons (McGraw-Hill 
Book Company, Inc., New York, 1948) p. 416. 
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resonators, and operating voltage range. For the 
voltage at which operation in a mode just begins, 
the quantities V/V» and B/ By are connected by 
Hartree’s equation 


V/Vo=2(B/ Bo) —-1. 


The beginning of oscillation is obscured by effects 
which spoil the operation at very small currents. 
Since, however, the dynamic impedance of mag- 
netrons is small, the above relation can be taken 
to hold approximately for the operating voltages. 
The quantities Vo, Bo, and J» are related to the 
dimensions and the wave-length of the mag- 
netron as follows: 


Vo =2.52X10°(D, nvd)°* volts, 
21,400 


By=— —— gauss, when \ is in cm, 
(1—o*)mdr 


5.23 x 10° Davi fh 
ji gS 
1 mh dD, 
a-ev(41) 
o 


where / is the anode lengthf and a, is a function 
of ¢, nearly equal to unity in the range of o of 
interest here. A plot of a, as a function of ¢@ is 
given by Slater.’ 

All of the magnetrons constructed and tested 
in these experiments were designed for the same 
values of \ and Vo. Small variations in \ and V5 
resulted from errors in resonator dimensions 
arising from the use of approximate calculations 
and from errors in construction. Operation in the 
x-mode (n= N/2) was the objective in all cases. 
Other values of m enter into the problem when the 
possibility of competition from other modes is 
considered. 

The initial values of o chosen for magnet- 
rons with an untried number of resonators 
was based on Slater’s semi-empirical formula 


o=(n—2)/(n+2). After experimentation it was 


found that the best compromise between effi- 
ciency and stability occurred very nearly at 
ao =0.94[(n—2)/(n+2) ] in all cases. 

The possible advantages that can be derived 


t An effective value of h is used in the computation of J 
for closed-end tubes, in which the r-f electric field intensity 
varies sinusoidally along the length of the anode. The 
complication of the step construction in some of the tube 
types has also been included in the effective values. 
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from the use of a large number of resonators as- 
suming the above design relationship to hold are 
shown in Fig. 10, where the quantities By, D. and 
cathode current density are plotted as a function 
of N for magnetrons of the same A, Vo, and / and 
operating at equal values of J/Io. N values which 
have been actually used are indicated. The in- 
crease in cathode size and the reduction in cur- 
rent density resulting from large N should help to 
extend practical limits on scaling down in wave- 
length, on reducing operating voltage, or on 
operating at higher power input. In the latter 
connection it is of interest that one form of 38- 
vane 1.25-cm magnetron has been driven to 
90-amperes pulsed current, yielding 300 kw of 
pulsed power output for a brief period. 

To evaluate the performance of each of the 
magnetrons tested, an extensive program of 
measurements was carried out. Since the major 
problem was that of obtaining stable 2-mode 
operation, a large part of these measurements 
was concerned with mapping out the regions of 
the performance chart in which such operation 
occurred, and the regions in which disturbances 
appeared. Wave-length measurements on the 
competing modes were made for identification 
purposes. While the unwanted modes made their 
presence known in different ways, depending 
upon particular tube design and on circumstances 
of operation, their most persistent manifestation 
in these experiments was in the form of pulse-to- 
pulse mode shifting. Here operation occurs in the 
mw-mode for the duration of some pulses, and in 
the unwanted mode during other pulses, the 
shifting occurring apparently at random. Ob- 
served symptoms of this kind of mode shifting 
were: (a) the appearance of apparently simul- 
taneous traces on an oscilloscope which records 
current pulses, indicating distinctly different 
values of pulsed current for each mode, (b) the 
random appearance and disappearance of lines in 
the presentation of the z-mode spectrum on the 
screen of a spectrum analyzer, and (c) two dis- 
tinct responses on the wave meter, one for the 
m-mode and one for the unwanted mode. The 
term “‘mode competition,’’ used in Sections II 
and III fits this situation particularly well. The 
competition occurs in the starting of every pulse, 
the outcome being in doubt until the pulse is well 
established. Since the rate of build-up of the 


JOURNAL OF APPLIED PHYSICS 














m-mode, for which the output circuit of the 
magnetron is designed, depends upon the ex- 
ternal loading, the stability of r-mode operation 
is dependent on the external loading as well as 
upon applied voltage and magnetic field. The 
operation over a reasonable range, in magnitude 
and phase, of external load impedance was there- 
fore investigated for stability. 

The experimental observations on the mag- 
netrons constructed show that the mode compe- 
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ANODE CURRENT—AMPERES 


(a) 18-vane magnetron: pulse duration—0.3 ysec.; duty cycle—0.0003; 
wave-length—1.29 cm. 


Fic. 11. Comparison of performance of 18-, 26-, and 38- 
vane rising sun magnetrons. (Electronic efficiencies are 
given in parentheses above the over-all efficiencies. ) 
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tition can be eliminated, by proper design, for 
values of N as great as 38, and there is no indi- 
cation that this is an upper limit. For N = 38 this 
was accomplished only with the closed-end 
design. For 26 vanes stable operation with open 
ends was observed for one particular set of di- 
mensions, but the closed-end design allows more 
latitude in the choice of design parameters, and is 
therefore to be preferred. 

The elimination of mode competition alone 
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(b) 26-vane magnetron: pulse duration—0.2 ysec.; duty cycle—0.0002 
(closed ends); wave-length—1.26 cm. 
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(c) 38-vane magnetron: pulse duration—0.2 ywsec.; duty cycle 
wave-length—1.3 cm. 


0.0002; 
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TABLE |. Design and operating data for best 22-, 26-, 
and 38-vane, 1.25-cm magnetrons, with good 18-vane tube 
for comparison. 











Number of 
resonators 18 22 26 26 38 
resonator ends Open Open Open Closed Closed 
cath. mount axial radial radial axial axial 
D, (in.) 0.165 0.181 0.218 0.234 0.339 
*D,’ (in.) 0.270 0.381 
h (in.) 0.190 0.170 0.170 0.309 0.350 
*h’ (in.) 0.220 0.272 
D, (in.) 0.403 0.410 0.457 fal - 
D, (in.) 0.300 0.311 0.349 0.386 0.479 
ry 1.76 1.76 1.82 2.6 2.68 
o 0.59 0.64 0.70 0.71 0.76 
Vo (volts) 3280 3030 3020 3350 3050 
Bo (gauss) 2840 2780 2620 2630 2060 
Io (amps.) 24 23 27 28 29 
Ax (cm) 1.290 1.207 1.230 1.257 1.300 
External Q 350 500 340 1030 800 


Unloaded Q 750 710 600 860 1000 


Performance at maximum efficiency 


Electronic eff. 


(percent) 43 43 43 44 43 
I (amp.) 11 9 17 12 12 
V/Vo 4.5 5.0 4.9 4.2 4.8 
B/Boy 2.8 3.1 2.8 2.7 2.9 
I/Io 0.47 0.40 0.63 0.44 0.42 


* Step dimensions, see Fig. 7. 
** See Fig. 9. 26-vane tube has Dy =0.096 and Dyc =0.490; 38-vane 
tube has Dy =0.079 and Dyc =0.586. 


may not be sufficient to make available the ad- 
vantages which have been predicted for large N. 
That these are indeed available may be demon- 
strated by a comparison, in terms of reduced 
variables, of performance charts of the large NV 
magnetrons with those of the 18-vane magnetron. 

A typical good 18-vane tube and best indi- 
vidual specimens of the 26- and 38-vane tubes are 
compared on the above basis in the performance 
charts in Fig. 11. Reduced variable scales, as well 
as the ordinary scales appear in the charts, and 
electronic efficiencies as well as ordinary eff- 
ciencies are given. 

Electronic efficiencies are given for comparison 
in the performance charts because of the compli- 
cation arising from variations in loading. Most of 


-the effect of this variation on the observed eff- 


ciency may be removed by comparing electronic 
efficiencies rather than over-all efficiencies. The 
electronic efficiency also varies with the loading, 
but slowly enough that uncorrected values of this 
quantity serve for a fair comparison. It may be 
seen that the three performance charts are quite 
similar, particularly in terms of the reduced 
variables and the electronic efficiency. There is, 
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in fact, no more variation among them than is 
found among samples of presumably identical 
magnetrons at this frequency. It may be con- 
cluded, therefore, that all expectations have been 
fulfilled in this scaling operation, and that opera- 
tion at N values up to 38 can be as good as at 
N=18. 

A criterion of performance that is of special 
interest in connection with scaling to short wave- 
lengths is the low field behavior, since at very 
short wave-lengths the lower values of B/Bo may 
be the only ones within reach. It may be seen 
from the performance charts that, in this respect, 
the 26- and 38-vane magnetrons are, if anything, 
slightly better than the 18-vane magnetron. 

Table I gives design and operating data for 
these three magnetrons and also for open-end 
22-vane and open-end 26-vane magnetrons. The 
performance of the magnetrons in the table for 
which performance charts are not given is well 
represented by the maximum electronic efficiency 
and the coordinates at which it appears on the 
performance chart. In all cases operation was 
stable x-mode operation with good spectrum over 
the region covered by the performance charts 
shown. 


VI. CONCLUSION 


From the results of the experiments reported 
here we may conclude that the rising sun design 
can be applied to the construction of magnetrons 
with a number of cavities considerably greater 
than 18. Successful operation for open cavity 
resonant systems has thus far been limited to 
magnetrons with 26 resonators in the voltage 
range in use for pulsed magnetrons, and can 
perhaps be extended to 30 resonators for low 
voltage CW magnetrons. When closed-end cavity 
resonators are used the number of cavities can be 
increased to at least 38, and there is no reason to 
believe that this number constitutes an upper 
limit. 

The advantages derived from the use of a large 
number of cavities are principally in the increase 
in cathode size and in some reduction in the re- 
quired magnetic field for magnetron operation. 
These advantages are, of course, very attractive 
when one is concerned with the design of high 
frequency generators in the millimeter range. It 
should be pointed out, however, that serious 
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difficulties remain in realizing satisfactory mag- 
netron operation in the millimeter range. In addi- 
tion to the increase in constructional difficulties 
there is the problem of decreasing distance be- 
tween cathode and anode for constant voltage 
difference, as the wave-length is decreased. More- 
over, the anomalously poor operation of the 
magnetron in the low current region may be 
aggravated by the high space charge densities 
that are expected with increasing magnetic field. 
It is hoped, however, that the use of a larger 
number of resonators for magnetrons will facili- 
tate the extension of the present magnetron 
techniques to wave-lengths of a few millimeters. 


The increase in the number of resonators can 
certainly lead to increased output power at 1 cm 
or longer wave-lengths in situations where the 
output power is limited by the cathode area. 
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18-22 American Institute of Electrical Engineers (Midwest General 
Meeting), Milwaukee, Wisconsin 
21-23 Optical Society of America, Detroit, Michigan 
25-28 Technical Association of the Pulp and Paper Industry, Buffalo, 
New York 
25-29 American Society for Metals, Philadelphia, Pennsylvania 
25-29 Society of Motion Picture Engineers, Washington, D. C. 
30 American Mathematical Society, New York, New York 
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Transmission Losses in 2n-Terminal Networks 


VitoLD BELEVITCH 
Bell Telephone Manufacturing Company, Antwerp, Belgium 


(Received October 23, 1947) 


The semipositive character of the dissipation function of 
any physical 2-terminal network imposes some restrictions 
on the realizability of prescribed transmission losses be- 
tween the various pairs of terminals when they are con- 
nected to given impedances. These restrictions are most 
easily expressed in terms of elements of the efficiency 
matrix of the network, which is defined in the first part of 
this paper. Discussions of the application of matrix calculus 
to the solution of reactive and resistive, transformer-type, 
and general 2n-terminal networks are presented. 

A very useful tool in solving telephone transmission 


HIS theory was suggested by an old paper 
by Campbell and Foster,' and a part of 
the following theorems is an extension of their 
results in modern notation using matrix calculus. 


I. DEFINITION OF TRANSMISSION LOSSES 


Assume a 2n-pole network with terminal pairs 
numbered (1,1’), (2,2’), ---, (m,n’). Any pair 
(p,p’) will be across a generator or a receiver 
impedance R,, presumed, as is usual in trans- 
mission theory, to be a pure resistance. The 
network is thus associated with a given set of 
terminating resistances. 

Suppose that a generator of voltage V, be 
connected in series with resistance R,, producing 
a current J, through resistance R,, when all 
other terminal pairs are across their respective 
resistances. We shall let S,, be the square root 
of the ratio between the complex power R,/,? 
delivered to the receiver, and the maximum 
power V,?/4R, available from the generator. 
Thus 

Sup =(2(RpRy)'Ia/ Vy). (1) 


The usual definitions of the effective loss A,, 
in nepers, and of phase shift B,, in radians are 


A gp t+ jByy = — logeSyp. (2) 


In the following we shall use instead the ex- 
pression S,,, and call it the complex efficiency 
factor between terminals p and q. 


1G. A. Campbell and R. M. Foster, “Maximum output 
networks for telephone substation and repeater circuits,” 
Trans. A.1.E.E. (February, 1920), 
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problems is the conventional theory of 4-terminal trans- 
ducers, often called 4-pole networks. Although this theory 
is sufficient for the study of most circuits, some problems 
involve networks with a larger number of terminals, which 
are often referred to as 2n-terminal, or 2n-pole networks. 
This paper is a contribution to their theory. 

The particular point that will be considered in this paper 
is the restriction imposed by physical realizability on the 
arbitrary choice of transmission losses between various 
pairs of terminals connected to a given set of terminating 
(generator or receiver) impedances. 


It will also be convenient to introduce a factor 
Spp, Which will be defined as the reflection coeffi- 
cient at the terminal pair (p,p’) ; 


Spp= (Rp—$p Ry+%p), (3) 


where ¢, is the impedance measured between 
the terminals (p,p’), with all other terminal pairs 
closed on their associated resistances. By writing 


A ppt JB pp = —log.Spp, (4) 


the usual definitions of the return loss and phase 
are obtained. The factor S,, is thus the square 
root of the complex ratio of reflected power to 
maximum power. The condition S,,=0 expresses 
an impedance match at terminals (p,p’). 

The matrix S, composed of elements S,,, will 
be termed the efficiency matrix of the 2n-pole 
network associated with a given set of termi- 
nating resistances. If the main diagonal of this 
matrix is zero, the network is simultaneously 
matched at all its terminal pairs. 

A simultaneous multiplication of the pth row 
and of the pth column by —1 is irrelevant, 
because it merely interchanges the terminals p 


and p’. 
II. CALCULATION OF EFFICIENCY MATRIX 


We start from the usual set of equations of a 
2n-pole network in terms of the self and transfer 
impedances ; 


Vp— Ryply=ZyrlitZypelet: ++ tZonlns (5) 


where p=1, 2, ---, m. In (5), for the sake of 
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venerality, it is supposed that generators were 
present in any terminal pair. 

Introducing normalized currents, voltages, and 
impedances by writing 


I,’ = (R,)'I,, Vp = V,/(Rp)', 
Zvq =Z/(RpR,)', (6) 


the system of equations in (5) becomes 
ey 7 | oe =Z yi Ty! +Z 2’ Te’ + eel +555... (7) 
Since we will always consider normalized 
quantities in the following equations, we may 
drop the primes and write (7) in matrix form: 
V—l=ZI, (8) 


where V and J are one-column matrices, and Z is 
the usual normalized impedance matrix of the 
network. 

Introducing the unit-matrix E, (8) may be 
written 


V=(Z+E)I/. (9) 
Solving for /, 
IT=QV, (10) 
where Q is the inverse matrix of Z+E; 
Q=(Z+E)". (11) 


A generator JV’,, acting alone in the terminal 
pair p, produces in the terminal pair g a current 
¢ as ? 4 * Si i ~ 4 - i . 
I,=Qy»V». Since normalized expressions are 
equivalent to making the terminating resistances 
equal to unity, (1) is reduced to 


Sap = 2Qop- (12) 


Similarly, Q,, is the normalized driving-point 
admittance of mesh p. If we denote by ¢,’ the 
normalized impedance ¢,/R,, (3) becomes 


Spp= (1 fp’) /(I +57’). (13) 
But, by Thévenin’s theorem, 
Opp =1/(14+8,'), (14) 
and, from (13) and (14), 
Sop = 2Qpp— 1. (15) 


Equations (12) and (15) may be combined in 
matrix form: 


S=20-E=2(Z+E)"-E, (16) 
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which may also be written 
S=(E-Z)(E+Z)“. (17) 


This is a formal extension of (13). 

In this calculation, we assumed that both 
matrices Z and Q exist, but this is not necessarily 
true. The degenerate cases will be examined 
later on. 

It is shown in (17) that the matrix S is sym- 
metrical, as is Z (reciprocity theorem). By re- 
placing Z with Z~ (or taking the inverse net- 
work), S is simply changed into —S, and this 
is a form of the duality principle. Inverse net- 
works give the same transmission losses. 


III. REACTIVE 2N-POLE NETWORKS 


In the case of a network composed of only 
reactive elements, the matrix Z is purely imagi- 
nary. Taking the complex conjugate 8S, Z is 
simply changed into —Z; 


S=(E+Z)(E-Z)". (18) 
By multiplying (17) and (18), 
SS=E, (19) 


and the efficiency matrix of a reactive 2n-pole 
network is therefore a complex orthogonal matrix. 
The diagonal terms of (19), written explicitly, are 


exp(—2A pi: +exp(—2A p2)+°:- 
+exp(—2A pn) = 1, (20) 


where p=1, 2, ---, 2. This demonstrates that no 
power can be consumed in a reactive network; 
it is an extension of a familiar relation first 
established by Feldtkeller for reactive 4-pole 
networks. 


IV. RESISTIVE 2N-POLE NETWORKS 


In the case of a purely resistive network, the 
matrix Z is real. It is easily shown that Z and S 
may be simultaneously reduced to diagonal 
matrices (respectively, Z and $) by means of a 
non-singular linear orthogonal transformation. 
Equation (17) is then transformed into 


Sop = (1—2ps)/(14+-Zop)- (21) 


The expression for the power consumed in a 
passive network is known to be of a semipositive 
quadratic form; when the quadratic form is 
transformed into a sum of squares, the coeff- 
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cients are Z,, and must be positive or zero. We 
thus have 


-—1¢5,, <1. (22) 


The case of a degenerate network having no 
impedance matrix may be considered as being 
the limiting case of a network for which at least 
one of the numbers Z,, becomes infinite, so that 
(22) must be extended: 


—1¢§,,<1. (23) 


In this last form the condition obviously ex- 
presses the fact that both matrices E+S and 
E—S are semipositive. This can be written 
explicitly in the usual way in terms of principal 
minors. 

As an example, we will consider the case of a 
hexapole simultaneously matched at all its 
terminal pairs. The efficiency matrix is 


| 0 Sir Sis| 
S= Sis 0 Ses|| . 


By noting that the determinants of E+S are 
non-negative, one obtains the two non-trivial 
conditions 

1 +2S 12893513 — Sis” — So3? — S12? 2 0. (24) 


If the hexapole has the same losses between all 
terminal pairs, i.e., if 


| Si2| = Si3| - | Ses | =, 
one of these conditions gives 
1 —3s?—2s*=(1+5)?(1—2s) 20, 
the other one being trivial. Thus 
s<}, 


and the minimum value of the loss is 0.69 neper 
=6 decibels, or twice as great as could be ex- 
pected from the division of the input energy 
between the two output resistances. Actually, a 
3-decibel loss is obtained in a hybrid coil, but it 
is not symmetrical with respect to all terminal 
pairs. 
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V. IDEAL TRANSFORMERS 


A network composed of only ideal transformers 
may be considered as being the limiting case 
either of a reactive or of a resistive network. As 
its efficiency matrix is real, (19) becomes 


S=E. (25) 


For an orthogonal symmetric matrix, the num- 
bers S,, are either +1 or —1. Since the sum of 
these numbers is equal to the sum of the main 
diagonal elements of S, the sum must be zero for 
a network simultaneously matched at all of its 
terminals, which is obviously impossible if the 
order of the matrix is odd. As a consequence, 
a network matching an odd number of terminal 
pairs and composed only of ideal transformers 
that match an odd number of terminal pairs does 
not exist. 


VI. GENERAL 2N-POLE NETWORKS 


The restrictions imposed on the efficiency 
matrix of a general 2m-pole network are obtained 
by stating that the real part Zea: of the im- 
pedance matrix is semipositive. From (16) or 
(17), 

Z=(E—S)(E+S)1=2(E+S)'-E, 
and 


Z eat =3(Z+Z) =(E+S8)"—E+(E+S)" 
=(E+8)-(E-SS)(E+S)-. 


Solving for E—SS and noting that (16) gives 
E+S=2Q, one obtains 


E—SS8 =4QZ rea. (26) 


Consider now a semipositive quadratic form of 
matrix Zrea: and apply to the variables the non- 
singular complex linear transformation defined 
by the matrix Q; the quadratic form becomes a 
Hermitian form in the new variables and has 
the transformed matrix (26). As the semipositive 
character is not changed by the transformation, 
the conditions of physical realizability can be 
obtained explicitly in terms of the efficiency 
factors by writing that the principal minors of 
E—SS must be positive or zero. The results 
previously obtained for reactive or resistive net- 
works are particular cases of these conditions. 
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Transmission of Monoenergetic Slow Neutrons through Solid Solutions and 
. Mechanical Mixtures of TiC and WC* 
> 
S. S. Sipnu** 
; Argonne National Laboratory, Chicago, Illinois 
(Received December 9, 1947) 
) TiC and WC form continuous substitutional solid solutions at least up to 25 percent TiC+75 
percent WC by weight when heated to 4000°F, TiC being the solvent and WC the solute. 
: Total neutron scattering cross sections of a series of finely powdered solid solutions and corre- 
f sponding mechanical mixtures of the two carbides were determined as a function of the mole 
L percent of WC from the measured transmission of monoenergetic neutrons with energies 
. 0.08138, 0.03479, 0.02056, and 0.01526 ev or wave-lengths \=1.004A, A\=1.536A, A=1.998A 
. and \=2.319A, respectively. The cross sections of both the solutions and the mixtures increase 
. with increasing mole percent of WC. For mixtures they are the same as the additive cross 
sections, each constituent scattering independently of the other. They are analogous to x-ray 
’ scattering for the mixtures in which each constituent gives its own diffraction pattern and the 
] patterns of the two constituents are superimposed. The neutron cross sections of solid solutions 
5 are lower than those of the mixtures at low percentages of WC but increase rapidly with the 
. concentration of the solute molecules. They appear to be sensitive both to the structure of 
the solvent and the concentration of the substituting atoms. 
1. INTRODUCTION tering cross sections of the material for slow 
1 HE transmission of neutrons with wave- "¢Utrons. Therefore, the change in the eae 
; lengths comparable to the interplanar ‘TOSS Section may be used as acer O we 
r spacings in a crystalline substance has been Presence and the degree of such actors. ie 
shown'~* to be modified by the physical state of In the present investigation apne ror: 
the substance. The scattering is due mainly to of monoenergetic slow neutrons t Hue ; ow, 
Bragg reflections from crystalline planes that powdered solid solutions and le wwe ci 
satisfy the usual Bragg relation, n\=2dsin@. The ‘TES of TiC and WC was ane and t cen 
scattering cross sections that are determined rmieet ane ype aa os 
from the measured transmission would conse- ‘U¢ton of the mole percent o : 
‘ quently vary with the number of Bragg reflec- 2. SAMPLES 
, tions. Any factor which would tend to cause a The samples studied were titanium carbide, 
change in the lattice of a crystalline material tungsten carbide, and a series of solid solutions 
) would also cause a change in its scattering cross and mechanical mixtures of the two carbides 
f section. Crystallinity, crystallization, grain size, made in the following percentage by weight: 
order disorder, orientations, strain, presence of oe son 
i- : ae P ‘ (1) 85 percent TiC+15 percent WC, 
d different phases, solute atoms in a solid solution, ; 
as ld h rye peal (2) 75 percent TiC+25 percent WC, 
a € i wou ave a mar effect on the scat- (3) 65 percent TiC+35 percent WC, 
“ * This document is based on work performed under the (4) 50 percent TiC+50 percent WC, 
e auspices of the Atomic Energy Commission at the Argonne (5) 35 percent TiC+65 percent WC, 
: National Laboratory. ( TiC+75 WC 
‘ ** On leave from the University of Pittsburgh. 6) 25 percent TiC+75 percent . 
ve 1M. D. Whitaker and H. G. Beyer, Phys. Rev. 55, 1101 ; ni 
‘ (1939). All the samples were as free from impurities 
f (1940) G. Beyer and M. D. Whitaker, Phys. Rev. 57,976 and moisture as they could be made with usual 
" * Franco Rassetti, Phys. Rev. 58, 321 (1940). care. They were ground to finer than 325 mesh 
t 58, AS (1980). G. Beyer, and J. R. Dunning, Phys. Rev. and packed into standard cylindrical aluminum 
°F. C. Nix and G. F. Clament, Phys. Rev. 68, 159 holders for cross-section measurements. 
(1945). -rav di i imen ob- 
+E. Fermi, W. J. Sturm and R.G. Sachs, Phys. Rev.71, _ /\" X-tay diffraction pattern of a spec . 
589 (1947). tained from each sample was made to determine 
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TABLE I. Total scattering cross sections for \=1.004A 
as a function of mole percent of WC in units of 10°% 
c m?/molecule. 


Mole percent o measured @ measured Additive* cross 


of WC for solution for mixture section 
5.1 13.1 14.5 14.1 
9.3 13.3 14.7 14.3 
14.1 13.4 14.8 14.6 
23.4 14.9 15.1 15.1 
36.2 16.4 16.0 15.9 
47.8 17.6 16.3 


16.6 


* Additive values were computed trom measured values of TiC and 
WC according to the proportion in which each compound was present. 


the crystalline state of the material and also to 
differentiate solutions and me- 
chanical mixtures. The sample 25 percent TiC 
+75 percent WC, which contained the largest 
amount of WC was a complete solid solution and 
well within the limit of solid solubility of WC in 
TiC at 4000°F, the temperature used for making 
the solid solutions. 


between solid 


3. MEASUREMENT OF TOTAL NEUTRON 
CROSS SECTIONS 

The method of measuring relative neutron 
transmission through the samples was essentially 
the same as used by Sturm,’ and by Sawyer, 
Wollan, Bernstein, and Peterson.’ The neutron 
Hux from the heavy water moderated Argonne 
pile was diffracted from the (111) planes of a LiF 
single crystal, and the neutrons satisfying the 
Bragg relation, m\ = 2d siné, where \=/h/ mv, with 
energies 0.08138, 0.03479, 0.02056, and 0.01526 
electron volts or wave-lengths \=1.004A, 
A\=1.536A, \=1.998A, and \=2.319A, respect- 
ively, were used. The transmission of the sample 
for neutrons of above wave-lengths was de- 








termined by setting the crystal at the glancing 
angle @ and the counter at an angle of 2@ to the 
incident beam. The counter used was a pro- 
portional counter filled with BF; enriched in the 
B” isotope. The counting rates were measured 
with and without the sample in the path of the 
diffracted beam first at the glancing angle @ and 
second with the crystal rotated through 1°, while 
the position of the counter was kept fixed. The 
transmission was then calculated from the four 
measurements according to the relation 


T/Ip=(Ri—1r1) (Re—roe), (1) 


in which R, is the counting rate of the beam 
transmitted through the sample and the sample 
container, R, the counting rate through the con- 
tainer only, and 7; and rz the corresponding 
counting rates after the crystal was rotated 
through 1°. The values of r; and r2 gave the back- 
ground due to incoherent scattering which 
reached the counter through the sample. (R2— re) 
=I], was taken as the intensity of the incident 
beam and (R,—7,;)=T7 as the intensity of the 
transmitted beam. 
Total cross sections were calculated from the 
relation 
I/Ig=e"**, (2) 


where o is the total neutron cross section of 
sq. cm per molecule and N is the number in 
molecules per sq. cm of the sample. The value of 
N for a solid solution or a mechanical mixture 


was obtained from the formula 
N=(wt. percent/molecular wt. of TiC 
+wt. percent/molecular wt. of WC) 
6.023 K 10" ‘surface area, (3) 


Fic. 1. X-ray diffraction photographs of TiC and solid solutions of TiC and WC. (A) TiC, 

(B) 50 percent TiC+50 percent WC, and (C) 25 percent TiC+75 percent WC. 
7 William J. Sturm, Phys. Rev. 71, 757 (1947). 
*R. B. Sawyer, E. O. Wollan, S. Bernstein, and K. C. Peterson, Phys. Rev. 72, 109 (1947). 
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Fic. 2. X-ray diffraction photographs of mechanical mixtures of TiC and WC and WC. (A) 50 percent 
TiC+50 percent WC, (B) 25 percent TiC+75 percent WC, and (C) WC. 


and of o from: 


o =(1/N) In(Ip/I). (4) 
4. RESULTS 


X-ray diffraction photographs of TiC and 
some of the solid solutions are given in Fig. 1, and 
those of the corresponding mechanical mixtures 
and WC are given in Fig. 2. The total cross 
sections of the solutions and the mixtures for 
\=1.004A, as a function of the mole percent of 
WC are given in Fig. 3. Similar curves were ob- 
tained for other wave-lengths used here. For 
comparison the measured and the additive cross 
sections are given in Table I. 


5. DISCUSSION OF RESULTS 


Both the x-ray diffraction patterns and the 
total scattering cross sections of the solid solu- 
tions and the mechanical mixtures show several 
points of interest. It is quite apparent from Fig. 1 
that TiC is the solvent and WC the solute, and 
they form continuous substitutional solid solu- 
tions at least up to 25 percent TiC+75 percent 
WC when heated to 4000°F. 

The measured cross sections of both the solu- 
tions and the mixtures, as shown in Fig. 3, 
increase with increasing mole percent of WC. For 
mixtures they are practically the same as the 
additive cross sections. as shown in Table I. It is 
apparent that each constituent scatters inde- 
pendently of the other. This is quite analogous to 
X-ray scattering from the mixtures. As shown in 
Fig. 2 each constitutent gives its own diffraction 
pattern and the patterns of the two constitutents 
are superimposed. 

The scattering cross sections of solid solutions 
are different from those of the mixtures or the 
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additive cross sections. They are lower than those 
of the mixtures at low percentages of WC, but 
increase rapidly with the concentration of the 
solute molecules. This effect was observed for 
neutrons of all four wave-lengths used in this 
investigation. Since all. solid solutions were 
heated to the same temperature and powdered to 
the same mesh this increase could not be ascribed 
to the grain size and the particle size of the 
material. 

The measured scattering cross sections of TiC 
and WC for neutrons of \= 1.004A were 13.8 and 
19.6 barns, respectively. Since for low concen- 
trations of WC in TiC the scattering cross 
sections are lower than those of the corresponding 
mixtures and increase rapidly with the solute 
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Fic. 3. Total cross sections of microcrystalline (A) solid 
solutions and (B) mechanical mixtures of TiC and WC as 
a function of mole percent of WC. Neutron energy 81.38 
10-3 ev or wave-length 1.004A. o in units of 10°** em? 
per molecule. 


molecules, it appears that they are sensitive both 
to the structure of the solvent and the concentra- 
tion of the substituting atoms. For a given 
structure, therefore, the scattering cross section 
may be used as a measure of the solubility of the 
solute. 

Mr. David Meneghetti helped in taking some 
of the readings. 
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Experimental Results on Standing Wave Type Linear Accelerators for Electrons 


E.utiott J. Lawton anp W. C. HAHN 
General Electric Research Laboratory, Schenectady, New York 


(Received December 18, 1947) 


The work that was done in this laboratory on the 
development of a resonant cavity type linear accelerator 
is briefly described. The experimental results demonstrate 
that electrons entering at low energies will accelerate to 
a higher energy in passing through a group of cavities of 
equal 8(Ao/2) spacing and along which appears a standing 
wave field distribution. Starting at 57 kv, an energy of 
0.345 Mev was attained with two cavities. An 8-cavity 
accelerator of cavity length corresponding to 8=0.75 had 
a maximum output of 0.57 Mev. Experimentally deter- 
mined energies were in fair agreement with theoretical 
values at equal axial gradients. Excitation power was 
found to be proportional to [(Mev)?/Length ](1/Q) X coup- 


INTRODUCTION 


HE advent of the high power r-f pulse type 

magnetron equipment and techniques de- 
veloped during the war in the 2000- to 3000- 
megacycle range has been the cause for a renewed 
interest in the linear accelerator principle that 
was first tried by Wideroe' in Germany in 1928 
and later in this country by Sloane and Lawrence? 
and Beams and Snoddy.* The peak power output 
of some of these new magnetrons is in the 
megawatt range and so made energies of the 
order of 1 Mev/ft. appear practical. 

At these frequencies, instead of an individual 
accelerator gap being formed between the ends of 
two cylindrical electrodes, it becomes part of a 
cavity resonator. A series of such resonators 
joined together in line at B(Ao/2) interval form 
the accelerator. 8= V/C, where V is the electron 
velocity for the given gap spacing and C the 
velocity of light. A» is the resonance wave-length 
of the structure. 

Operation depends on establishing a standing 
wave field distribution along the axis of the line 
of cavities such that a complete phase reversal 
takes place over each fundamental cavity length. 
An electron whose velocity matches the phase 


a Wideroe, Archiv. fur Electrotechnik 21, 387-406 
1928). 
(1933) H. Sloan and E. O. Lawrence, Phys. Rev. 38, 2031 
1). 
J. W. Beams and L. B. Snoddy, Phys. Rev. 44, 784 
(1933) and 45, 287(A) (1934). 
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ling coefficient. The power required for an axial gradient 
of 1 Mev/ft., even for the most favorable design, was high 
and would require the full output of one HK7 type mag- 
netron for each foot of accelerator length. Energy spread 
in the accelerated beam was large for the two cavity tubes 
but small for the longer 8-cavity tube. It is believed that 
the apparent low energy spread observed in the 8-cavity 
case resulted from a limited range in entering angle brought 
about by the reduced axial gradient. There was little 
evidence of radial defocusing. A-change in resonance wave- 
length amounting to +0.000225 cm per degree centigrade 
was observed. There was indication that the magnetron 
would follow a resonance change of approximately 5°C. 


velocity of this wave will experience an ac- 
celerating force at each gap or cavity. 

Electrons can be injected at practical velocities, 
8=0.19 to 0.41 corresponding to 10 to 50 kv, 
providing the spacing of the first few gaps is so 
adjusted that the electron receives a net ac- 
celerating force at each gap that tends to syn- 
chronize it with the phase velocity of the exciting 
wave. This spacing must be systematically in- 
creased until the electron has reached a nearly 
constant velocity approaching that of light. It is 
not necessary, however, to increase progressively 
the spacing of each individual gap, as there are 
strong-phase focusing forces in the beginning 
along a group of equally spaced gaps that tend to 
synchronize the velocity of the particle with the 
wave. 

A theoretical study of this phase of the 
problem on the differential analyzer indicated 
that for an axial excitation field of 1 Mev/ft. 
three group spacing changes would be necessary 
to bring the electron up to approach the velocity 
of light ; namely, spacings correspond to 8=90.75, 
0.96, and finally 1.0. Eight of the 8=0.75 and 11 
of the 6 = 0.96 cavities were indicated to bring the 
electrons up to the approximately 2-Mev entering 
velocity for the 8=1.0 final section. For phase 
focusing and bunching, the electron, on entering 
a given cavity, must be in advance of the peak 
accelerating field. Figure 1 shows some typical 
differential analyzer runs for these conditions 
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with 1-Mev/ft. gradient and 50-kv injection. 
Curve A shows the variation of energy of the 
electron as it traverses the accelerator, for a time- 
entering angle 45° in advance of the peak of the 
a.c. gradient. Likewise B corresponds to 90° and 
C to 135° angular advance. From these it is seen 
that the stable phase position corresponds to an 
entering angle between 45° and 90°. A corre- 
sponding study for 0.763 Mev/ft. shows that 
there was no such stable entering angle. 

The holes in the cavity washers set up radial 
components of the electric field, so that an elec- 
tron on passing through the cavity will see a 
focusing field on entering and a defocusing field 
on leaving. Unfortunately, the condition that the 
particle be slightly in advance of the peak ac- 
celerating field for bunching results in its passing 
through a radial field on leaving the cavity that 
is larger than the entering focusing field. The 
conditions for phase focusing and bunching are 
therefore incompatible with those for radial 
focusing. 

The associated magnetic field is of such a direc- 
tion as to produce a focusing action during the 
phase focusing period and therefore counteracts 
the defocusing effect of the electric field. Whether 
the counteracting effect of the magnetic field 
would balance the defocusing effect of the electric 
field is problematical and would depend some- 
what on the final configuration. The experimental 
evidence for the first eight cavities points to a 
favorable balance. 

Since many mutually coupled resonant cavities 
are involved, there will be other undesirable 
modes of oscillation that must be avoided in 
order not to have interference between the 
operating and adjacent modes. This interference 
will set an upper limit to the length of the ac- 
celerator. It has been demonstrated‘ that the 
greatest mode separation is near the center 7/2 
mode of the over-all +1 modes of the mode 
spectrum, making this mode best for operation. 
Mode separation can be further increased by 
symmetrical excitation at the center of an array 
of equal number of cavities where the field dis- 
tribution is maximum for the center mode but 
zero for the two adjacent modes. Choice of the 
™/2 mode and a practical value of coupling 


* Elliott J. Lawton, J. of App. Phys. 19, 534 (1948). 
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K=9.05 percent and a Q=9000 leads to a 
limiting accelerator length of 10 to 20 ft. at 10.7- 
cm wave-length. This length could probably be 
extended in proportion to the square root of the 
wave-length through the use of a longer operating 
wave-length. The power, required to maintain a 
given axial gradient, however, would need be 
increased in proportion to the (A)! because of the 
resultant lower Q. 


EXPERIMENTAL 


The schematic arrangement of an 8-cavity ac- 
celerator tube and associated equipment are 
shown in Fig. 2. The photograph of Fig: 3 shows 
the experimental set-up with a 2-cavity ac- 
celerator section in the test position. 

The resonant cavities of the accelerator are 
formed by joining together alternate thin copper 
washers W, and large diameter ring spacers S at 
regular B(Ao/2) interval (see Fig. 2). The number 
of washers W for each interval establishes the 
desired mode of oscillation while the dimension of 
washer and spacer determines the wave-length 
and coupling K between cavities. K is the wave- 
length difference between the lowest and highest 
modes expressed as a percentage of the wave- 
length of the lowest order mode. To establish the 
desired x/2 operating mode two partitions were 
necessary in each 8(Ao/2) interval. Dimensions 
were usually determined experimentally and, be- 
cause of the high Q’s involved (9000), machining 
tolerances of the order of +0.0005” necessarily 
had to be maintained. 

Electrons entered and left the accelerator 
through terminating cylinders of diameter less 
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Fic. 1. Differential analyzer run for a 1-Mev/ft. gradient 
and 50-kv injection. @= time entering angle. 
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Fic. 2. Schematic diagram 
of an 8-cavity UHF linear 


ae accelerator for electrons. 
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INJECTION GUN ACCELERATOR CAVITIES 
than the “cut-off” diameter for the resonance 
wave-length. By adjusting the amount that the 
end of the terminating cylinder extended into the 
end cavities, the termination could be made to 
appear as a short circuit. The termination could 
be either at a mid-plane through the cavity be- 
tween two partitions or through the partition 
proper. In the case of short tubes (two cavities) 
the partition termination resulted in a higher Q 
than the cavity type because of the one less 
partition required. 

R-f power was supplied from a 1-megawatt 
pulse type HK7, 10-cm tuneable magnetron. A 
pulse width of 3.44 seconds and at a repetition 
rate of 262 cycles per second was used. A mini- 
mum pulse width of 1.54 second was found 
necessary to allow for the build-up time due to 
the high Q of the cavity and for the r-f field to be- 
come steady inside the cavity. 

Power was transmitted from the magnetron to 
the accelerator through 14’’X3" wave guide. 
Excitation of the accelerator cavities was through 
a hole in the outer diameter of one of the cavities. 
The size of this hole was so adjusted that the 
unloaded Q, equalled the external Q.x:, the con- 
dition for match. 

The photographs of Figs. 4 and 5 show, re- 


spectively, a completed 2-cavity accelerator of 


the type terminated at a point midway between 
partitions and of coupling K =17.75 percent and 
a 2-cavity tube terminated at a partition and of 
coupling K =9.05 percent. The cavity length in 
both cases corresponds to 6=0.75. In Fig. 5, the 
wave guide connection has been removed and 
shows the coupling hole, which in this case 
necessarily had to be centered about the partition 
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SYSTEM SCREEN 


MAGNETIC SPECTROMETER 


for this type of feed and termination. In the case 
of the tube in Fig. 4, the coupling hole lies mid- 
way between the two partitions. Figure 4 also 
shows the wave guide connection to the tube and 
resonant type glass window used as the vacuum 
seal between the wave guide and the accelerator 
cavities. 

For stable operation of the magnetron at the 
desired resonance wave-length, it was necessary 
to include a series impedance between magnetron 
and cavity (see Fig. 2). It was necessary to 
dissipate approximately 40 percent of the total 
magnetron power in this series load in order to 
destroy the condition for resonance at more than 
one wave-length between the two oscillating cir- 
cuits, magnetron, and loaded transmission line. 

The above duty cycle and the unavoidable 
amount of power that necessarily had to be 
dissipated in the series stabilizing load limited the 
useful power from the HK7 to approximately 0.3 
megawatts at the tube. 

Injection was with a standard x-ray tube type 
gun. Focusing of the injection beam was ac- 
complished with an external magnetic field. A 
steady d.c. injection beam current of the order of 
0.2 to 15 ma was used. 

The energy of the accelerated beam was meas- 
ured with a simple magnetic spectrometer. The 
deflection of the accelerated beam was deter- 
mined by observing its position on a fluorescent 
screen. This also provided a means of determining 
energy spread in the accelerated beam. The 
spectrometer was calibrated with the d.c. injec- 
tion beam. The screen was observed from a dis- 
tance with a telescope because of the danger of 
x-rays. 
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The over-all tube including gun, accelerator, 
and spectrometer was evacuated to a pressure of 
the order of 2 to 5X10-° mm of Hg with two oil 
diffusion pumps. No extensive preliminary out- 
gassing was found necessary. The vacuum seal at 
the wave guide connection to the tube was with a 
resonant type glass-to-fernico seal window.® Me- 
chanical alignment of gun, accelerator, and 
spectrometer was accomplished optically by 
sighting through sets of peepholes along the axis. 
Final electrical alignment, to compensate for 
stray magnetic fields which were quite objection- 
able, was taken care of by shifting the gun and 
associated focusing coil. 

Surrounding the accelerator section were mag- 
netic field coils for compensating for any radial 
defocusing. 


EXPERIMENTAL RESULTS 


The experimental results are confined to ac- 
celerator sections made up of cavities of length 
corresponding to 8<1 and which represent the 
initial gun stages necessary to bring the electron 
up to the entering velocity of the final 6=1 
stages. 

Four different accelerator sections based on the 
7/2 mode design were tested, namely: No. 1, an 
eight-cavity tube of cavity length corresponding 
to 8=0.75; No. 2, a two-cavity tube of the same 
design and coupling K=17.75 percent as the 
eight-cavity tube; No. 3, a two-cavity tube of the 
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Fic. 3. Experimental set-up with a 2-cavity accelerator in 
the test position. 


5M. D. Fiske, Rev. Sci. Inst. 17, 478 (1946). 
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same cavity length as 1 and 2, but of coupling 
K=9.05 percent and terminated at a partition 
instead of midway between partitions; and No. 4, 
a two-cavity tube of the same coupling and termi- 
nation as 3, but of cavity length corresponding to 
8=0.5. 

The experimental arrangement was such that 
the magnitude of the accelerated beam currents 
could not be measured. 

The experimentally determined energies were 
found to be in fair agreement with the theo- 
retically predicted values where it was possible to 
compare the two at equal axial gradients and in- 
jection voltages. The comparison is shown in the 
Table I. 

In the case of tubes 1 and 2, the maximum 
gradient was limited to the values listed because 
of limited r-f output power from the HK7 
magnetron. 

In the case of the 8-cavity tube 1, whose design 
was based on a 1-Mev/ft. gradient, there was evi- 
dence of a critical phase relation between beam 
and wave at the maximum possible gradient 0.52 
Mev/ft. and a complete lack of phase focusing at 
0.44 Mev/ft. Below this value there was a sharp 
decrease in the output voltage to a value not 
much in excess of the injection voltage. The 
theoretical study indicated an expected value of 
the minimum gradient for phase focusing that 
would be only slightly lower than 0.763 Mev/ft. 
The total output voltage as a function of excita- 
tion power for this tube is shown for different 
injection voltages in Fig. 6. 

The minimum practical gradient for phase 
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Fic. 4. Two-cavity accelerator, showing terminating cylin- 
der and wave guide connection. K = 17.75 percent. 
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Fic. 5. Two-cavity accelerator, showing coupling hole 
between wave guide and cavity. Terminating is at a 
partition. K =9.05 percent. 





focusing and bunching as predicted by the theo- 
retical study was 1 Mev/ft. (see Fig. 1). The 
power required to excite the 1 ft. long, 8-cavity 
tube to 0.52 Mev/ft. corresponded to 0.360 
megawatt/ft. Since excitation power appears as a 
copper loss in the walls of the cavity it should be 
proportional to 


[(Mev)?/Length ] 1/Q. 


The excitation power for the 8-cavity tube 
extrapolated to the value required for 1 Mev/ft. 
would be (1/0.52)?X 0.36 = 1.34 megawatts/ft., a 
prohibitively large value. These results were con- 
firmed by the tests of the 2-cavity tube No. 2 of 
the same coupling and design. The ratio of the 
voltage output of the 8-cavity tube to that of the 
2-cavity tube at the same injection voltage and 
power level was 2.18, which is in agreement with 
the value 2.0 as predicted by the power relation. 

The tests of the 2-cavity tube No. 3 of lower 
coupling K=9.05 percent indicated a larger 
voltage output at a given excitation power. The 
total voltage output vs. total excitation power at 


different injection voltages is shown in Fig. 7. In 


Fig. 8 is shown the voltage per gap, or average 
field, after deducting the injection voltage, as a 
function of the excitation power per gap. For 
comparison, the field results for the higher 
coupling K = 17.75 percent, tube No. 2, are shown 
in Curve A. It is to be pointed out that the 
excitation values given for curve A are corrected 
for Q to compare with the 9000 value of tube 
No. 3. The Q of tube 2 was 6000. It will be noted 
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that the curve for the tube 3 is drawn through the 
group of points for injection voltages between 41 
and 93 kv where the field values are closely 
grouped. Below this range, 30 to 12 kv, the 
entering phase relations were probably such that 
the beam did not pick up enough energy in 
falling through the first gap and so emerged at a 
lower value. This was somewhat borne out by the 
behavior of the beam intensity, as at these lower 
values the injection beam current had to be in- 
creased by a factor of approximately 5 to 10 times 
to produce a visible accelerated beam on the 
fluorescent screen. It can be seen that the excita- 
tion power required for 4 given field, in the case 
of tube No. 3 with lower coupling, is approxi- 
mately one-half of that required for tube No. 2 
which had double the coupling, indicating that 
excitation power is directly proportional to the 
coupling over this coupling range. The excitation 
power required for a 1-Mev/ft. (120 kv/gap) 
gradient for this tube is 0.624 megawatt /ft., and 
would represent the power required for the first 
gun section of a limiting length of accelerator of 
10 to 20 ft. 

The experimental results for the shorter cavity 
(8=0.5) tube No. 4 of same design and coupling 
as tube No. 3 indicate that at a given excitation 
power 0.300 megawatt and 60-kv injection, the 
maximum output was only 0.190 Mev as com- 
pared to 0.390 Mev for the longer cavity higher 
Q tube No. 3. The Q of this tube was 6280 as 
compared to 9000 value for tube 3 and is lower in 
proportion to the decrease in volume due to the 
shorter length cavity. Since the coupling was the 
same in the two cases, the ratio of the two output 
voltages, at constant power level, would be ex- 
pected to be the same as the square root of the Q 
times length ratio as predicted by the power re- 
lation. The agreement is only fair and is 2.05 for 
the voltage ratio as compared to 1.47 for the QL 
ratio. The discrepancy is not apparent, but might 
be associated with slightly different electron 
dynamics in the two cases. The results are suffi- 
ciently close, however, to allow an extrapolation 
of the data on a QL basis, and so lead to a de- 
termination of the approximate power require- 
ments for the intermediate 8=0.96 gun section 
and final 8=1.0 stages of the accelerator. 

The approximate excitation for these stages at 
1-Mev/ft. gradient and extrapolated from the 
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0.624-megawatt/ft. value for tube No. 3 are 0.306 
and 0.270 megawatt per ft. for the intermediate 
8=0.96 and final 6 = 1.0 stages, respectively. The 
corresponding expected Q’s would be 11,700 and 
12,000, respectively. 

For a 10-ft. accelerator, which length would 
approach the limiting value for the 7/2 mode 
design and of coupling K=9.05 percent, the 
average excitation power required to maintain a 
1-Mev/ft. gradient would be 0.319 megawatt/ft. 
In terms of the number of HK7 magnetrons re- 
quired, this would amount to approximately one 
for each foot of accelerator. 

While it has been demonstrated by the M.I.T. 
Accelerator Group® that stable parallel operation 
of magnetrons up to three magnetrons is possible, 
the usefulness of this type of accelerator as a 
generator of high energy particles and x-rays is 
greatly complicated by the large number of 
magnetrons and associated equipment involved. 
The successful outcome of the development of 
higher power magnetrons now in progress would 
greatly alter this picture. 

The energy spread in the case of the two cavity 
tubes tested was large and of the order of 3 to } of 
the total output at high values of injection cur- 
rent. The apparent energy spread decreased as 
the injection beam current was decreased until at 
low values the spot was no larger than the injec- 
tion beam spot, indicating more current coming 
through at the optimum phase. Although the 
energy spread for the 8-cavity tube was ex- 
tremely low (accelerated spot same size as in- 
jection beam spot), it was felt that this was due in 
part to a limited range in entering angle brought 
about by the reduced gradient, rather than a 
direct indication of improved energy spread re- 
sulting from an increased number of cavities. The 
spread, however, should decrease with increasing 
number of cavities because of the greater time 
allowed for bunching. The theoretical study indi- 
cated that at 1-Mev/ft. gradient an energy 
spread of the order of 10 percent to 20 percent 
would be expected after the electron had trav- 
ersed the first 3 ft. of tube. It is true that electrons 
starting at different points along the length would 
come through at varying energies, however, only 


__* Winston Bostick, Edgar Everhart, and Melvin Labitt, 
lech. Report No. 14 of Research Laboratory of Electronics, 
Massachusetts Institute of Technology. 
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TABLE I. Comparison of experimental with theoretical 
accelerated voltages at given axial gradients. 








Experimental Theoretical 
_ Gradient Max. Ex- 
Inj. attained Mev Inj. Gradient 
Accelerator section kv Mev/ft. output kv Mev/ft. Mev 


(1) 8cavity, 8=0.75 length 57 = 0.52 0.57 50 = 0.763 = (0.625 
K=17.75 percent 





(2) 2-cavity, 8=0.75 length 65 0.72 0.234 50 0.763 80.190 
K=17.75 percent 


(3) 2-cavity,° 8=0.75 length 58 0.763 0.250 50 0.763 0.190 
K =9.05 percent 58 1.00 0.315 50 1.00 0.250 
58 1.33 0.390 8650 1.33 0.345 











°=terminations at a partition 


those starting near the first few cavities would 
have a chance to synchronize with the wave. 

A mild form of r-f breakdown inside of the 
cavity occurred at a field of approximately 0.06 
Mev/gap. This was observed as a general back- 
ground illumination on the fluorescent screen. 
Both the background intensity and the x-ray 
output increased with excitation power. Also a 
slight increase in pressure was observed during 
the initial operating stages but disappeared after 
a short period of operation. The effect was usually 
encouraged by an axial magnetic field, which was 
probably due to a restriction of electron motion 
to the region of highest electric field. There was 
no indication that this breakdown was detri- 
mental, as its occurrence produced no change in 
the r-f envelope shape. 

In the case of tube 3, a well defined spot of 
energy approximately corresponding to that for 
the given excitation accompanied the breakdown, 
and so provided 4 means of determining the 
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FiG. 6. 8-cavity, cavity terminated, 8=0.75 tube total 
Mev vs. total megawatts input at different injection 
voltages. Ao = 10.692 cm (vac) (20°C); @=8000; K =17.75 
percent. 
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Fic. 7. Two-cavity, partition terminated, 8=0.75 tube 
total Mev vs. total megawatts input at different injection 
voltages. 


voltage at which complete breakdown occurred 
at different pressures. At a given cavity field, as 
the pressure was allowed to increase slowly, no 
change in spot deflection or r-f envelope shape 
was observed until very near the critical pressure 
for complete breakdown, which probably indi- 
cated that there was no predominant glow dis- 
charge period preceding the critical breakdown 
pressure. At the critical pressure the accelerated 
spot disappeared, the x-ray output decreased to 
zero, and the build-up and decay periods of the 
r-f envelope both collapsed. The critical pres- 
sures for two values of r-f field 0.09 Mev/gap and 
0.14 Mev/gap (1.15 Mev/ft.) were 4.2 and 2.0 
microns, respectively. 

The resonance wave-length of the accelerator 
was observed to change with temperature and by 
an amount corresponding to +0.000225 cm per 
degree centigrade in the case of the 8-cavity tube 
No. 1. For a Q of 10,000 and A» = 10.7 cm, a 4.8°C 
change should reduce the output to one-half. 
Actually the output remained essentially con- 
stant over approximately this range as_ the 
magnetron appeared to follow the resonance 
change. Temperature changes, in the case of 
short tubes such as would be used for x-ray appli- 
cation, should not be a serious consideration. In 
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Fic. 8. Two-cavity, partition terminated, 8=0.75 tube 
Mev vs. megawatt input for one cavity. \o= 10.701 cm (vac) 
(20°C); Q@=9000; d= 1.366"; K =9.05 percent. 


the case of long tubes, however, where it is im- 
portant to maintain a constant cavity length so 
as to maintain proper phase relation between 
beam and exciting wave, temperature changes 
along the length of the accelerator would become 
a serious consideration. 


CONCLUSIONS 


The main conclusion is that although electrons 
can be accelerated to energies corresponding to a 
gradient of 1 Mev/ft. by this means, the power 
required to maintain this gradient is sufficiently 
high as to require an alarmingly large number of 
present day r-f power supplies such as the HK7 
type magnetron. This viewpoint will change as 
large supplies of r-f power become more readily 
available. 

Although the limiting length, set by mode 
interference, appears to be between 10 and 20 ft. 
at 10 cm, it is felt that this could be extended in 
proportion to the square root of the wave-length 
used. However, in doing so, the power required 
would be increased in proportion to the (A)}. 
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A Method for Measuring the Complex Dielectric Constant of Gases at Microwave 
Frequencies by Using a Resonant Cavity* 


C. K. JEN 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 
(Received January 8, 1948) 


A sensitive method has been developed for the measurement of the complex dielectric 
constant of low pressure gases at X-band microwave frequencies. The gas, whose complex 
dielectric properties are to be determined, is contained in a resonant cavity which is a part 
of the microwave circuit. With the introduction of the gas, the real part of the dielectric 
constant changes the resonant frequency of the cavity, while the imaginary part changes the 
amplitude and breadth of the cavity response curve. By a rapid variation of the frequency 
across the cavity resonance, the real and imaginary parts can be conveniently and accurately 
determined from measurements on an oscilloscope. Tests on such a system have been performed 
with satisfactory results. Sample data are presented for methylchloride and deuterated 
ammonia. This method can be used for the measurement of the resonant dispersion and 
absorption of microwaves by gas molecules. Also, the same method can be used for the measure- 


CS 





ment of loaded and unloaded Q-values of a cavity. 


INTRODUCTION 


HE experimental study of the complex 

dielectric properties of materials as a func- 
tion of frequency has had a long history. There- 
fore, it is only natural that such measurements 
have now been extended to the microwave fre- 
quencies as a result of the rapid development of 
these very high frequency techniques in recent 
years. Both wave guides and resonant cavities 
containing the dielectric media have been com- 
monly used for this purpose.'~* Until quite re- 
cently most of the attention has been directed to 
solid dielectrics and comparatively little to gases, 
but the tremendous progress made in microwave 
spectroscopy with gases is exciting considerable 
interest in measurements in gaseous media. 
Notable examples are found in the measurements 
of the absorption coefficient of O2 by Beringer,° 
of NH; by Bleaney and Penrose,* and of NH; and 
many other gases by Hershberger.’ 


* The research reported in this document was made 
possible through support extended Cruft Laboratory, 
Harvard University, jointly by the Navy Department 
(Office of Naval Research) and the Signal Corps, U. S. 
Army, under ONR Contract N5ori-76, Task e. I. 
The contents of this paper have been reported earlier in 
Progress Report No. 2, for Contract NS5ori-76, Cruft 
Laboratory, January 1, 1947. 

!W. R. Maclean, J. te 7 Phys. 17, 558 (1946). 

2 T. W. Dakin and C. N. Works, J. App. Phys. 18, 789 
(1947). 

3C. N. Works, J. we Phys. 18, 605 (1947). 

4S. Roberts and A. V. Hippel, J. App. Phys. 17, 610 
(1946). 

5 R. E. Beringer, M.I.T. Radiation Lab. Report No. 684. 
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The object of the present work is to develop a 
more sensitive and accurate method of measuring 
at microwave frequencies the complex dielectric 
constant of gases at fairly low pressures. Such a 
method should be usable both for the non- 
resonant and the resonant dielectric properties. 
The X-band (8500-9750 Mc) is chosen partly be- 
cause of the availability of apparatus and partly 
because of the anticipation that in the near future 
the X-band will also be a fruitful field for 
microwave spectroscopy. The method is, how- 
ever, equally applicable to other frequency 
ranges. A resonant cavity has been chosen for this 
study because it is equivalent to a very long wave 
guide at its resonant frequency. 


THEORY OF MEASUREMENT 


A resonant cavity behaves in principle like a 
low frequency resonant circuit. Its properties are 
usually characterized by a loss factor (in this 
paper it will be called 6, which is the reciprocal of 
the customary Q-factor) and a resonant fre- 
quency. When a resonant cavity is coupled to a 
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Fic. 1. Coupled cavity. 


°B. Bleaney and R. P. Penrose, Proc. Roy. Soc. 189, 
358 (1947). 


7™W. D. Hershberger, J. App. Phys. 17, 495 (1946). 
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eu Hected wave from the cavity to that of the inci- 
"| dent wave, we have 
=t}-8 ia c¢  cawiTy do 2 f—fo 
— oe rt 
OUTPUT Y¢ soe Y Cc 61 5; fo 
: ake r=— = —__. -- -———-, (3) 
Fic. 3. Cavity coupled to magic tee. Yet+ Yc 5p f-—fo 
(+: 5, j ) 
wave guide, as shown in Fig. 1, the situation is 6 db fo 


again analogous to an ordinary coupled circuit at , 
low frequencies. We can then derive by this The vectorial plot of I is a circle as shown in 
analogy the admittance of the cavity (a “reflec- Fig. 2. In order to translate the amplitude of the 


tion” type in this case) as seen from the wave reflection factor into the usual form of a resonant 


guide, with the following result: response, we now adopt a so-called magic-tee 
arrangement as shown in Fig. 3. Here, the inci- 

bo 1 f?—fo’ dent wave enters the magic tee at A, and then 

Yce= | 5 ft Ff, | (1) divides equally into branches B and C. At the end 


of branch B, there is a shorting plunger which 








where always has a reflection factor [=—1 in the 
Y¢=admittance of the cavity as seen from the immediate vicinity of the short. Depending upon 
wave guide, ; the position of the plunger, this reflected wave 
Ye=characteristic admittance of the wave can assume any phase difference relative to the 
guide, wave reflected from the cavity when the two re- 
energy loss in cavity wall flected waves combine at D. Furthermore, by the 
plus loss in medium per radian __ property of a magic tee, the vector difference of 
— es energy stored the two reflected waves passes through branch D, 
energy escaping from the while the vector sum returns through branch A. 
; input window per radian _ Referring to the incident wave present in either 
= - coc” 


energy stored branch B or branch C as unit amplitude, let 

fo=resonant frequency of the gas-filled 
cavity, 

f=frequency of the incident wave. 


F=vector difference of the two reflected waves 
when branches B and C are equal in 
microwave phase, 

In the neighborhood of the cavity resonance, i.e., G = vector difference of the two reflected waves 

. f*=fo, we have when branches B and C differ by 7m in 

microwave phase ; 


do / 2 f—fo 
Yo= Yej —+j— ~- (2) then we have 
61 6; fo 
F=-1-T, (4) 
Defining I as the reflection factor which denotes G=+1-T. 








650 JOURNAL OF APPLIED PHYSICS 





eS 


in 





It follows that 
1 
[1+ (50/81) 2 +L(2/51)(f—fo/ fo)? 
[ (60/81) P+0(2/81)(f—fo/ fo) 
[1+ (50/61) 2+[(2/61)(f—fo/ fo) 


At f=fo, |F\? is a maximum while |G|? is a 
minimum. Thus, we have 





|Fit=4 





| F| max 1G| min = 51 bo. (5) 


Referring to Fig. 4, at | F/ Fmax|?= 3, the “‘half- 
power” frequency difference is related to the 6- 
values as follows: 


(f" —f')/fo= bot 41. (6) 


Now, the left-hand terms of Eqs. (5) and (6) are 
two directly measurable quantities. For brevity, 
let 


x= (f" — f'/fo), 
te | F| max/|G| min; 


then we can determine 49 directly from the 
following expression : 


bo = (x/1+y). (7) 

Thus far we have not considered the differential 

effect of the gas. Suppose, to start with, the 
cavity is empty. Then in Eq. (2), we have 


fo=fe, 
bo=8., (8) 


To D.C. Repeller 
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where the subscript e refers to quantities in an 
empty cavity. When the cavity is filled with gas, 
we must consider the changes resulting from the 
presence of the complex dielectric constant. Let 
the dielectric constant be expressed as follows: 


e=e' — je’, (9) 


where, for most gases, e’’<<e’ = 1. The admittance 
of a gas-filled cavity may then be written as 




















_f 
& 2 (e)! 
Ye= Ye] —+j-— —— 
51 61 4 
{ (e)* J 
— x fe 4 
E+” 2 (e’)! 
=Y, eee ee ams (10) 
61 61 Is 
i (’)) J 








Comparing Eq. (10) with Eq. (2), we see that 
with a gas-filled cavity, 


= f./(e')! 
a 


We now conclude from Eq. (8) and Eq. (11), that 
for e’ very close to unity 


2 
‘’=1{-— (fo) gas — (fo)empty J, (12 
: Gouna J - : ' 
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Fic. 5. Complete schematic diagram and CRO pattern. 
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and 


” ies = ( 50) ens — ( do)empty e ( 13) 


From Eqs. (12) and (13) we see that the real part 
of the dielectric constant can be calculated 
readily from the shift of the cavity resonance fre- 
quency and the imaginary part from the change 


in the loss factor as defined in Eq. (7). It should | 


be remarked that the procedure of determining 
e”’ in Eq. (13) with the aid of Eq. (7) is essen- 
tially different from the apparently simpler pro- 
cedure of using the following expression from 
keq. (6): 


"i = (d0+ 51) eas po (do0+ 5: )empty- 


The latter statement would be true if 6; remains 
sensibly constant. Since the coupling factor 6, 
may vary with the introduction of gas, the loss 
factor 59 must be isolated from the total effect, 
thus necessitating the general 
veloped in Eqs. (4) to (7). 


treatment de- 


EXPERIMENTAL ARRANGEMENT 


The experimental arrangement for carrying out 
the proposed measurements should consist es- 
sentially of one part which yields an amplitude- 
frequency display of the cavity resonance, and 
another part which serves to measure frequency 
differences on the same display. These objectives 
can be accomplished in a variety of ways, de- 
pending upon the degree of sensitivity required 
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Fic. 6. e'—1 and e’”’ values for methyl-chloride 
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and the availability of instruments. It was found 
convenient for the present experiment to adopt a 
scheme as shown in Fig. 5. 

The cavity tested here is of a type ordinarily 
used as an X-band wave meter, having an inner 
diameter of 4.952 cm. The cavity is tunable in the 
range 8500-9900 Mc. It has an unloaded Q-value 
of about 25,000. For evacuation and gas-filling 
purposes, the whole cavity is enclosed in a 
cylindrical brass housing, which is sealed on to 
the wave guide leading into the cavity. A sealed 
glass tube projecting from the same enclosure 
allows the cavity to be tuned by a magnetic 
micrometer head by means of an external magnet. 
The cavity is sealed from the wave guide by a 
window consisting of a quartz disk of about one- 
mm thickness with a silver-chromium evaporated 
film bordering the rectangular opening, which is 
soldered with Wood’s metal to a special choke 
ange. To cancel the effect of the wave reflected 
by the quartz window, a similar window is placed 
symmetrically with respect to the magic-tee so 
that the reflected wave will not appear at the 
output. The microwave source is provided by a 
2K25 reflex klystron, whose frequency can con- 
veniently be swept in any desired manner. 

For the measurement of frequency differences, 
use is made of another 2K25 klystron, frequency- 
stabilized by Pound’s method.* The output is 
applied to a push-pull modulator, which mixes 
the microwave with a radiofrequency wave of a 
few Mc, thus producing two sideband frequencies 
with the carrier almost completely suppressed. 
One of the sideband frequencies can be con- 
veniently used for measuring frequency differ- 
ences since its position can be continuously 
changed relative to the stabilized carrier. 

Measurement of the voltage ratios, either for 
determining the ‘“half-power”’ points or the ratio 

F | max/\|G\ min, is accomplished by comparison 
with a known voltage ratio in a slotted section of 
wave guide, in which the field distribution may 
be assumed to be sinusoidal. 

Both the signal and frequency-measuring 
waves are then detected, amplified, and finally 
displayed on an oscilloscope. An X-band spec- 
trum-analyzer,’ designed by M.1I.T. Radiation 


8 R. V. Pound, Proc. I.R.E. 35, 1405 (1947). 
* For an alternative scheme, see L. B. Young and L. G. 
Mann, Internal Report 55.5, M.I.T. Radiation Lab. (1945). 


JOURNAL OF APPLIED PHYSICS 





sly 


for 
tio 
son 
| of 
lay 


ing 
lly 
eC 
ion 


45). 


ICS 





Laboratory, has been used for this purpose. Since 
the frequency of the local klystron of the spec- 
trum analyzer is already swept at 60 cycles per 
second, it is necessary to sweep the signal fre- 
quency at a rate which is very much lower or 
higher than this. We have actually used a 
frequency of about 50 kc, which frequency 
modulates the klystron, as shown in Fig. 5. This 
operation results in an oscilloscopic pattern of 
closely spaced lines whose envelope represents the 
amplitude of the input wave. The frequency- 
measuring side frequency is superimposed on this 
pattern and is movable with respect to it by 
adjustment of the r-f oscillator. The whole pat- 
tern is also shown in Fig. 5. 

In the measurement of the complex dielectric 
constant, two similar operations are performed, 
one with gas in the cavity and the other with the 
cavity evacuated. In each case, the center of the 
resonance curve is matched with one of the side- 
band frequencies and then the frequency differ- 
ence of the “half-power’”’ points is measured. The 
ratio | F| max/|G| min is measured by adjusting the 
shorting plunger in the gas cavity assembly to 
show alternately a | F|?-curve and a |G|?-curve. 
Then, | F| max? and |G| min? are equalized by an 
attenuator at the input to the spectrum analyzer. 
The attenuator is calibrated from the standing- 
wave pattern in the slotted wave guide section of 
the stabilized assembly. The final results are then 
calculated from Eqs. (7), (12), and (13). 


SAMPLE RESULTS OF MEASUREMENT 


The experimental result for the sensitivity of 
the arrangement described above may be ex- 
pressed in terms of the minimum detectable 
change for e’ — 1 and e’’. This value is estimated to 
be roughly 10~° for e’—1 and 2X10~-* for e’’. 

As a sample of a number of measurements, re- 
sults for the complex dielectric constant of methyl- 
chloride (CH;Cl) and of a mixture of deuterated 
ammonia (NH;;: 12.5 percent, NDH2: 37.5 per- 
cent, ND2H: 37.5 percent, ND ;: 12.5 percent, by 
volume) are presented in Fig. 6. 
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DISCUSSION 


The above experimental results show that the 
proposed method can be successfully employed 
for the measurement of the complex dielectric 
constant of gases. From the magnitude of the ab- 
sorption factor (e’’) and its dependence upon 
pressure, it can be concluded that the resonant 
absorption frequencies of deuterated ammonia 
are much nearer to the working frequency than 
those of methyl-chloride. It is well known that 
NH; exhibits ‘‘inversion’’ absorption for a band 
of frequencies around 24,000 Mc, and ND; is ex- 
pected to absorb somewhere around 6000 Mc 
according to Manning’s calculations.’ It is 
reasonable to assume that the intermediate 
species have absorption frequencies between these 
limits. What we get here is the over-all effect of a 
number of selective absorptions by the mixture. 
For methyl-chloride the measured result may be 
due to resonant absorption’ at frequencies some- 
what higher than the normal K-band and partly 
due to non-resonant absorption (the zero-fre- 
quency contribution). 

It can be seen that the present method is 
equally applicable to the measurement of reso- 
nant absorption and dispersion and, further, it 
should be able to predict the general location of 
the absorbing center from the dependence of the 
complex dielectric constant upon pressure and 
frequency. 

From the way in which this experiment is per- 
formed, the loaded and unloaded Q-values are 
necessarily known quantities. Thus the method is 
also useful for the determination of such values, if 
desired. 
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Carbon Film Formation and Commutator Brush-Wear As Revealed by the 
Electron Microscope 
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The faces of carbon or graphite brushes, polished by 
rubbing against a rotating copper disk and subsequently 
examined with the electron microscope by means of surface 
replicas, are found to consist of many tilted “‘fingers,”’ 
which are the projecting portions of the large and irregular 
graphite plates composing the brush. These graphite pro- 
jections show a directional overlap, which is determined 
by the direction of motion of the disk, and when viewed 
stereoscopically they are found to be inclined at various 
angles up to nearly 90° with reference to the plane of the 
rotating base. They account for a previously observed 
“directional memory” of graphited surfaces and they show 
that such surfaces in general cannot be regarded as con- 
sisting simply of graphite plates oriented with the hex- 
agonal planes flat and parallel to the plane of motion. 

The brush track on the rotating copper base is found to 
contain small islands of carbon deposited from the brush. 


INTRODUCTION 


HE carbon film which forms on collector- 

ring or commutator surfaces has been 
isolated and identified chemically,' and its im- 
portance is now well known. However, its struc- 
ture—the arrangement of the crystals within the 
film—has not been revealed. The film has been 
presumed to consist of graphite platelets oriented 
with their slip-planes flat on the surfaces and 
parallel to the direction of motion. This view has 
been drawn from analogy with the established 
structure of “graphoid films’’ on bearing sur- 
faces.” 

Lately it has been suggested* that this view 
might not be correct but that the platelets may 
overlap in the direction of stroke to form a 
“shingled” layer. The initial evidence for this 
new viewpoint was based upon observations of 
steel ball bearings lubricated with graphite (con- 
taining small amounts of Apiezon grease) and 
operated in a sliding contact apparatus* in 
vacuum. It was found that the bearings operated 
quietly, with low friction, when revolving in the 








''C, Van Brunt, unpublished reports (1937); also C. Van 
Brunt and R. H. Savage, Gen. Elec. Rev. 47, 16 (1944). 

2 R. O. Jenkins, Phil. Mag. 17, 457 (1934). 
*R. H. Savage, J. App. Phys. 19, 1 (1948). 
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These islands consist of large numbers of graphite wear- 
fragments, which are individually of very small size in 
comparison with the projections of the brush face. These 
fragments are clearly derived from the projections by a 
tearing apart of their edges and are packed tightly and with 
random arrangement into the tool marks of the copper. 
smoothing over the irregularities of the latter. It is not 
determined whether this ‘‘carbon film’”’ of the brush track 
is continuous. . 

The wear-dust of unlubricated graphite is found to con- 
sist of particles which are very small in average size com- 
pared with those from which brushes are manufactured. 
The size of these wear-fragments extends below the re- 
solving power of the electron microscope, and this is 
consistent with recent data on similar particles of a size 
determined independently by adsorption measurements. 


usual direction of rotation (clockwise) but that 
after reversal of rotation the bearings temporarily 
were noisy and showed high friction (at least 
3 Xnormal as estimated from deceleration meas- 
urements). In the course of a few minutes of 
operation in the new direction the bearing opera- 
tion returned to normal. It was possible to repeat 
the high friction effect merely by reversing the 
rotational direction. The apparatus showed a 
directional memory and was even sensitive to a. 
single slight displacement (e.g., a fraction of one 
revolution) in the direction producing frictional 
disturbance. 

It was possible to understand this curious effect 
by assuming the overlapping already referred to. 
Subsequent experiments with carbor brushes 
operating upon slip-rings have shown supporting 
evidence of this directional effect. For example, 
the films formed under graphite brushes con- 
taining small amounts of inorganic salts, and 
operating with low friction under high altitude 
conditions in which less than the adequate con- 
centration of water vapor was present,' were 
suddenly torn up simply by a reversal of rotation, 
the brushes then wearing out within a few minutes 
(unless given additional water vapor.) 


4R. H. Savage, Gen. Elec. Rev. 48, 13 (1945). 
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The practical and fundamental importance of 
the film structure thus became evident, and work 
has been done to learn more about the associated 
phenomena using the electron microscope. During 
this work it was found essential to investigate 
also the structure of the faces of brushes which 
were used to form the films, and the carbon dust 
formed by wear of brushes in vacuum equipment. 
The initial report on this work follows. 

PREPARATION OF CARBON SAMPLES 

Representative collector-films and _ polished 

brush faces were made by operating graphite 








Fic. 1. Carbon 
motor-brush replica. 
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brushes in a friction apparatus on a copper disk 
rotating at 1750 r.p.m. This disk had been turned 
in a lathe to as smooth a face as could be con- 
veniently machined in the laboratory. Abrasives 
and polishing powders were excluded, and pre- 
cautions were taken to prevent contamination of 
the face after cleaning. The brushes were small 
rods of square cross section, machined from 
electrographitized carbon blocks made from a 
coke base and having the following physical 
characteristics: resistivity 1800 microhmcentime- 
ter, scleroscope hardness 50, density 1.75. They 
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bore normally against the disk with apparent 
bearing areas of 0.06 cm® under spring loads of 
0.092 kg. They did not carry current. The atmos- 
phere was water vapor at the pressure of 7 mm 
(mercury), and the wear rate was observed to be 
nearly zero during the entire run, which was 
limited to one hour in order to limit the film 
thickness to partial transparency. Following this 


run the brush surfaces showed characteristic 


specular polish, and the unmagnified track was 
light gray against pale copper in general appear- 
ance. The brushes and disk were removed from 


Fic. 2. Stereo of carbon motor-brush replica. 


the apparatus with considerable care so as to 
avoid any reversing motion upon the contact 
surfaces, which would tend to disturb the surface 
arrangement of the crystals and particles. Im- 
mediately following, these samples were trans- 
ferred to covered petri dishes for examination. 


ELECTRON MICROGRAPHS 


To employ the high resolution and great depth 
of field of the electron microscope for this exami- 
nation of the surfaces that had been in contact 
with one another, it was necessary to use surface 
























Fic. 3. Stereo of carbon film on copper surface. 
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replicas. These replicas were prepared in the 
following manner: a 10 percent solution of puri- 
fied nitrocellulose in amyl acetate was placed on 
the surfaces to be examined and the solvent 
allowed to evaporate slowly. When the nitrocellu- 
lose film was hard, it was stripped mechanically 
from the surfaces, removing with it that material 
not firmly attached to the surface and repro- 
ducing in replica all other structures. Onto this 
surface a thin film of silica was evaporated in a 
vacuum by heating 1 mg of silica in a tungsten 
coil filament at a distance of 6.5 cm from the 
resin replica. The desired area of the replica was 
then selected under a light microscope and an 
electron microscope specimen screen fastened 
over it with the thin silica film in contact with the 
screen wires. The original replica was then dis- 
solved away and the silica replica washed by the 
condensation of amyl acetate vapors directly on 
them. This operation leaves a thin silica replica 
over the openings of the specimen screen, re- 
taining without displacement or alteration what- 
ever material has been removed from the original 
surface by the nitrocellulose. Care must be 
exercised in the interpretation of an electron 


micrograph of such a combination of graphite and 
replica-of-graphite structures. 

Figure 1 is a replica of a representative small 
area of one brush face. The magnification may be 
calculated from the superimposed micron scale. 
In this micrograph we see that plate-like particles 
comprising a portion of the surface protrude 
above the plane of the surface and overlap in the 
direction of motion of the copper, indicated by 
the arrow, thus forming many overlapping 
tongues or projections. These vary in diameter, in 
respect to surface, from less than 0.02 micron to 
particles larger than one micron, some of which 
are torn off by the original replica. Their number 
and distribution may vary over wide limits, as 
indicated by other pictures. The average area of 
the projections may be estimated at about 0.5 
square micron although the great majority are 
much smaller with areas of the order 0.01 square 
micron. The edges of the projections are ragged 
and suggest that the process of normal friction is 
accompanied by continual tearing away of frag- 
ments from the edges. 

In view of the extreme raggedness of the 
projections on the brush face, the term “‘shingles”’ 




















Fic. 4. Stereo of carbon film on copper surface. 
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is now misleading, but it was originally used be- 
fore any real conception of surface structure had 
been obtained, simply to suggest a directional 
overlapping. 

Figure 2 is a stereoscopic electron micrograph 
of another brush-face replica. It is similar to 
Fig. 1, but stereoscopically one may see that the 
overlapping tongues are inclined to all angles up 
to nearly 90 degrees from the base plane, with 
most at less than 45 degrees. From this the sug- 
gestion is made that the areas of contact consist 
of individual surfaces formed by these tongues 
which are isolated tilted planes scattered over the 
entire face. 

Electron diffraction patterns, contrary to these 
observations, have generally been interpreted as 
_ indicating that the crystals lie flat and parallel to 
the surface. However, such patterns do not reveal 
the detailed structure at the surface. 

Figure 3 is a stereoscopic electron micrograph 
of a replica of the carbon film on the copper 
surface from which some loose carbon particles as 
aggregates have been transferred to the thin 
replica. The replica itself shows the tool marks in 
the copper and also reproduces the structures 
formed by the carbon particles that are tightly 
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Fic. 5. Stereo of carbon motor-brush replica. 


packed into these marks over scattered areas. 
Many irregularities of the copper surface evi- 
dently collect the carbon particles and are quickly 
encased by firmly packed particles. These pro- 
trusions then become relatively smooth, after 
which they remove little more material from the 
brushes. Some of the protrusions are no doubt 
broken, and rebuilt by the desposition of other 
particles. In the study of this type of brush no 
indications have been found that the copper 
contributes to this smoothing process by a 
localized flow, as might occur under the influence 
of increased temperature or pressure. 

The direction of motion of the copper may be 
determined by the shape of some of the pro- 
trusions and by the preferred location of the 
loose carbon particles on the trailing side of these, 
as may be observed in the stereoscopic pictures. 
The direction of motion of the copper is indicated 
by the arrow in Fig. 3. 

The pioneering work of Van Brunt' on carbon 
films clearly indicated a porous* but nevertheless 
continuous structure under the optical microscope 
at 100 diameters’ magnification. Relatively large 


* Pores were not actually observed, but the reagents 
quickly penetrated the film and attacked the underlying 
base. 
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sheets were isolated intact, and were found to be 
transparent over small thin areas, but were 
streaked and mottled in general appearance be- 
cause of the opacity of thicker areas which were 
themselves discontinuous although contained 
within the whole film. The electron micrograph at 
high magnification (Fig. 3) was at first interpreted 
as indicating a discontinuous film structure, the 
only apparent film consisting of replicas of the 
isolated carbon-particle patches and clusters of 
carbon particles themselves. (The fine lines 
clearly represent the roughness elevations and 
depressions, or secondary tool marks in the 
copper.) However, such an interpretation of this 
micrograph is difficult to reconcile with the 
optical observations already referred to. Further- 
more, it does not account for the observed fact 
that the nitrocellulose is much more easily 
stripped from the brush track than from the 
untouched copper. Accordingly, two alternative 
possibilities should be recorded. (a) The carbon 
film over the copper surface may be incompletely 
developed in the present instance. Its rate of 
formation may have decreased enormously after 
the first few traverses of the brush during which 
time enough islands were formed to support the 











Fic. 6. Graphite stock powder No. 1107. 
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brush face continuously. The film which had been 
found to be continuous under the optical micro- 
scope! had been prepared by a much longer run of 
72 hours compared with the present film of 1 
hour. Or (b) an extremely thin carbon film may 
actually cover these tool marks, conforming to 
them in such a way as to indicate no superim- 
posed structure such as is shown by the large and 
semi-opaque graphite crystals occasionally noted. 

Figure 4 is another stereoscopic electron micro- 
graph of a replica of a brush track on the copper 
surface and shows in this instance a few graphite 
crystals which are very large and undistorted 
compared to the mass of particles on the surface. 
Here again the great majority of the particles 
have the same size and structure as the wear dust 
(see later) and were apparently produced by a 
breaking away from the edges of the tongues on 
the brush face, since their size is comparable with 
the irregularities of the tongues. 

None of the micrographs of the brush track has 
shown evidence of even partial orientation of the 
graphite crystals relative to the base plane. The 
indicated crystalline arrangement is completely 
random wherever visible. 

From the foregoing electron micrographs we 
are led to conclude that the brush face normally 





Fic. 7. Wear dust from No. 1107 brush. 
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contains a myriad of individually isolated graph- 
ite projections (tongues) stroked in the direction 
of motion so as to lie at varying obtuse angles to 
the base plane. Thesg projections evidently pro- 
vide the surfaces of true contact. They are 
undoubtedly flexible and elastic, and should sup- 
port the load as a unified system which is proba- 
bly unique in its resiliency and adaptability when 
in contact with irregular metal surfaces rotating 
at high speed. 

The process of wear is expected to be a normal 
accompaniment of the friction between the brush- 
face projections and the collector surface and may 
occur intermittently as a gradual erosion of these 
projections until, as a limit, flat areas have de- 
veloped. We should then expect renewed deposi- 
tion of graphite over these flat and presumably 
inelastic areas or a tearing of them, under the 
excessive pressures of nonconforming surfaces, to 
uncover fresh projections. Film formation and 
wear are thus fundamental counterprocesses ac- 
companying normal friction. Evidently, upon 
reversal of rotation a major rearrangement of the 
projections is required, which accounts for the 
directional memory of graphited surfaces as 
observed. 


SUPPLEMENTARY OBSERVATIONS 


An additional study has been made to check 


and to extend the foregoing observations. In this 


a second copper disk and new brushes were used. 
These were operated first in vacuum with high 
wear rates (track 1) in order to fit the rather 
irregular brush faces to the disk surface, then on a 
fresh area in a water vapor atmosphere at 6-mm 
(mercury) pressure for 6 minutes (track 2) to 
form a relatively thin film, and, finally, in the 
water vapor for 4 hours (track 3) to form a ‘“‘well 
developed”’ film. The electron micrographs of the 
surface replicas made from these new samples 
confirmed the initial observations in all of their 
fundamental aspects. Small tilted projections on 
the brush face were again noted, and these were 
characterized by the directional overlap. Inci- 
dentally it was noted that the 6-minute film was 
similar to the 4-hour film, indicating that the 
time factor within these limits had not been an 
important variable. One difference was noted, 
however, between the first micrographs and the 
later ones: in the latter the projections were 
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seemingly more rounded (Fig. 5). The suggestion 
is made tentatively that this difference may have 
resulted from slight structural or particulate 
differences between the different carbons which 
were used, or even between different successive 
cross sections of the same carbon, or, on the 
other hand, from different evaporation rates of 
silica from the filament during the formation of 
the replicas. Much more work is needed in order 
to clarify points such as this, and it is suggested 
that one approach which should also yield much 
other useful information will be through the use 
of carbons which differ greatly in their apparent 
physical properties, such as electrographitized 
black-carbon in comparison with crystalline 
natural graphite. 


WEAR DUST 


In order to understand more clearly the nature 
of the wearing process in relation to film forma- 
tion and to film structure, attempts were made to 
determine the extent to which the individual 
particles in the brush are subdivided as a result of 
the wear. For this purpose it was essential to have 
reference micrographs of the particles initially 
making up the brushes and comparison micro- 
graphs of the dust produced by the wear of these 
same brushes. It was not possible to obtain 
satisfactory reference micrographs for the brushes 
which were used in the film study reported in the 
preceeding discussion because those brushes con- 
tained a large proportion of graphite derived 
from the coked binder; this had been crystallized 
from organic material during the final firing 
process so that the determination of its structure 
and form was greatly complicated. Therefore in 
place of the electrographitized brush a study was 
made of a natural-graphite type. This type had 
been made by firing a pressed crystalline natural- 
graphite powder (No. 1107) with about 2 percent 
coked binder, to about 1350°C. 

Figure 6 is an electron micrograph of the No. 
1107 graphite stock powder before the brush was 
pressed from this material. The sample was 
mounted by dispersing the particles in a nitro- 
cellulose solution and drawing down the sus- 
pension on glass with a spatula.’ The particles 
are seen to be predominantly very large in area, 


a M. C. Schuster and E. F. Fullam, Ind. Eng. Chem. 
Anal. Ed. 18, 653 (1946). 
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and so thin in many instances as to be penetrated 
by the electron beam. The characteristic hex- 
agonal plate-like crystal structure may be easily 
identified. 

Figure 7 is the comparison micrograph and 
shows the dust produced by wear of the No. 1107 
brushes in the vacuum equipment. This dust 
shows a few large particles that have been chipped 
from the brush, and a large number of fine par- 
ticles occurring in a size range which evidently 
extends below the resolving power of the electron 
microscope. It is thus evident that the process of 
wear and removal of debris results in a subdivision 
of the original graphite particles to very small 
fragments. 

Electron micrographs of the wear dust made 
from the electrographitized brushes used in the 
film study were very similar to Fig. 7, showing a 
few chips of large size and a great number of 
small fragments. Many of the latter were of a 
size and thickness corresponding to edge irregu- 
larities and thickness of the projections of Fig. 1. 
From these observations it may be concluded 
that the process of wear involves the continual 
tearing away of the tongues of the brush face into 


many smaller fragments. Many of these frag- 
ments may be broken further into still smaller 
particles before being thrown from the copper 
disk. (Such an abrasive and crushing action as 
might be expected may partly separate at the 
edges the sheets of graphite comprising the 
particles so as to expose a large surface area.). 
Specific surface measurements® of the dust by 
nitrogen adsorption at — 195°C indicate a surface 
area of about 400 sq. m/gram. 
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The Radiation Patterns of Dielectric Rods—Experiment and Theory* 
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A theory is developed describing the electromagnetic radiation from dielectric rods used as 
terminations to rectangular wave guides. The radiation patterns are computed from the fields 
due to equivalent magnetic and electric currents on the surfaces of the rods. These equivalent 
currents are determined by use of Schelkunoff's equivalence principle applied to approximate 
fields within the dielectric. Specific calculations are carried out for rods whose average dimen- 
sions are 0.9’ 0.25" and whose lengths vary from three to ten wave-lengths. Experimental 
patterns are obtained for eight polystyrene rods, some of which are variable in length. Agree- 
ment with theory is found for rods from three to six wave-lengths long; differences between 


experiment and theory are discussed. 


INTRODUCTION 
A THEORY has been developed describing 


the electromagnetic radiation from dielec- 
tric rods used as terminations to rectangular 


-*The work described in this paper was done at the 
Defense Research Laboratory under the sponsorship of the 
Bureau of Ordnance, Navy Department, Contract NOrd- 
9195. 
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metallic wave guides. In particular, the theory 
leads to expressions for the radiation patterns 
for certain dielectric rods, which are in agree- 
ment with experimental observations for rods 
not longer than six wave-lengths. 

The theory is an application to dielectric 
rods of the equivalence principle stated by 
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Schelkunoff.! This principle leads to the estab- 
lishment of equivalent magnetic and electric 
currents over the surface of a particular rod; 
from these currents the radiation pattern may 
be computed. The equivalent currents are calcu- 
lated from the tangential components of the 
electric and magnetic vectors at the surface of 
the rod. In order to carry out this calculation 
exactly, the electromagnetic field in the di- 
electric must be known. Although the problem 
of finding this field has not been solved, reason- 
able approximations may be made for the field. 
The radiation patterns thus computed agree with 
experiment. Specifically, patterns are calculated 
for a dielectric rod whose average dimensions 
are 0.9’ 0.25” and whose length varies from 
three to ten wave-lengths. This dielectric rod is 
inserted in a wave guide which is excited in the 
T Eo, ; mode. 

Experimental patterns are obtained for eight 
different polystyrene rods of rectangular cross 
section. These patterns agree well with the 
corresponding ones based on theory; detailed 
analysis is made of the differences in terms of 
the theoretical approximations and experimental 
variations. 


THEORY 


The equivalence principle states (in part) that 
the electromagnetic field inside a surface = due 
to sources outside the surface can be produced by 
sheet electric currents J and sheet magnetic 
currents M over > given by 


J= —n XH’, 
M=nxXE’, 


(1) 
(2) 


where n is the unit normal vector directed out- 
ward from , X denotes the vector product, and 
E’ and H® are the values of the electric and 
magnetic vectors E and H on the surface 2. 
In terms of these equivalent currents, the mag- 
netic vector potential A and the electric vector 


Ex \\ 


Fic. 1. Dielectric rod as termination for 
rectangular wave guide. 











Firetree ceasing 


1 Bell Sys. Tech. J. 15, 93 (1936). 
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potential F at a point P within the surface are 


given by: 
1 exp(jkr’) 
na (Ps 
4nr Jy 


yr’ 


1 exp(jkr’) 
pot fa, 
4dr  - r’ 


where r’=distance from d> to the point P, 
k=2x/r, \=wave-length in free space, and 
j=(—1)}. In terms of these potentials, the electric 
and magnetic fields at the point P are: 


(4) 


E = jwpA— (1/ jwe) grad divA—curlF, (5) 
H = jweF —(1/jwu) grad divF+curlA, (6) 
where w=angular frequency, ¢=dielectric con- 


stant, w= permeability. These six equations allow 
calculation of the radiation pattern provided 
that E® and H® are known. 

In the application of these general statements 
to the problem of the radiation from a dielectric 
rod it is necessary to examine the detailed nature 
of the fields E®, H®. The dielectric rod serves as a 
termination for a rectangular metallic wave 
guide, as shown in Fig. 1. The dimensions of 
the rod are the same as the wave guide dimen- 
sions at the junction point and are assumed to 
vary only slowly with distance from the end of 
the wave guide. The surface corresponding to ¥ 
consists of the outside surfaces of the wave guide 
and dielectric rod and is closed by an infinitely 
large sphere about the origin. The sources are 
within the wave guide. It is assumed that for 
short distances (a few wave-lengths) outside of 
the wave guide, the field within the dielectric 
is not far different from the field in the guide. 
This is reasonable, for the dielectric is inserted 
within the guide a considerable distance. A ta- 
pered termination is included to prevent ex- 
cessive reflections and generation of harmonic 
waves at the junction between air and dielectric. 
Also, the shape of the dielectric is sensibly the 
same as that of the guide. Hence the field in 
the portion of the dielectric within the guide is 


CE TET LEE tat fae ea 








Fic. 2. Variation in electric field along 
length of dielectric rod. 
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accurately the same as the field in the empty 
guide; the considerable difference in dielectric 
constant between the dielectric and air should 
prevent excessive changes in the field in the di- 
electric outside of the guide, at least for short 
distances. 

The tangential components of the electri¢ and 
magnetic vectors at the surface Y are the only 
quantities giving rise to the sheet currents, as is 
indicated by Eqs. (1) and (2). The tangential 
components over the surface of the wave guide 
will be assumed zero. This is a very good approxi- 
mation for the electric vector, since the wave 
guide is a good conductor. The fields are further 
assumed to vanish at large distances, so that the 
only part of the surface Y on which the tangential 
components of the field are significant is that 
portion in contact with the surface of the 
dielectric. 


APPLICATION TO TE,,, MODE OF EXCITATION 


The transverse electric wave in the wave 
guide designated TE», is of especial interest. 
The electric vector has the form 


wy 
E.1= i(cos - ) exp( — j(wt— k’s)), (7) 


b 


where k’=27/\,, \y=wave-length inside the 
guide, i=unit vector in the x direction, and 
b=width of wave guide, and also of dielectric 
rod. The magnetic vector has the form 


ko omy «Kk ry 
He 1 = 1 j— cos—+—— gn— 
‘ ray b jp 


Xexp(— j(wt—k’z)), (8) 


where j and k are unit vectors in the y and z 
directions, respectively. 

The field at the surface of the dielectric will 
produce (in general) five categories of equivalent 
currents. The first three of these will be con- 
sidered in detail. 


(a) Magnetic currents on top and bottom (wide sur- 
faces). 

(b) Electric currents on the sides (narrow surfaces). 

(c) Electric currents on the top and bottom. 

(d) Magnetic currents on the sides. These currents arise 
from the tangential component of the electric vector at 
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Fic. 3. Coordinate systems used in pattern computations. 


these surfaces. The electric vector vanishes at the sides in 
the wave guide, and therefore cannot be of appreciable 
magnitude in the dielectric. These currents are then 
assumed negligible. 

(e) Electric and magnetic currents on the end. Because 
of the small radiation area of the end face compared to the 
other surface areas of the dielectric rod, the effect of any 
equivalent currents on the end is neglected. 


(a) Magnetic Currents on Top and Bottom 


The electric field inside the wave guide has 
zero tangential components on these surfaces; 
if the field in the dielectric were exactly the same 
as in the wave guide, the magnetic currents 
would yanish. There is undoubtedly a flow of 
energy through these surfaces, and hence it is 
concluded that the electric vector departs from 
the normal and has a tangential component 
parallel to the z axis. A suggested variation of the 
field is indicated by Fig. 2. The tangential com- 
ponents at top and bottom are in opposite direc- 
tions, but so are the corresponding vectors n; 
hence the magnetic currents M are in the same 
direction. The resultant magnetic currents are 
assumed to vary across the width of the guide 
as does the x component of the electric field, 
since it is the distortion of this component that 
leads to the currents. Hence 


M =j Mo cos(ry/b) exp(— jwt) sink’z. (9) 


Here k’ =22/dw, where \. is the wave-length in 
the dielectric. These currents are distributed 
over the length and width of the dielectric rod. 
The general form for the variation in the z di- 
rection, exp(jk’z), is particularized here to sink’z 
in order that the currents shall vanish at the 
opening of the wave guide as required by the 
absence of a tangential component within the 
guide. Since only the magnetic currents on 
the surfaces of the dielectric are to be considered, 
the electric currents are taken as zero, and, 
correspondingly, A from Eq. (3) is zero. The 
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value of F is given by Eq. (4) as: 


j wy 
F= - fto( 00s ) sinks 
4dr b> b 
exp( jkr’) 


Xexp( — jwt) : d>. 


r 


4 


(10) 


In evaluating this integral it is convenient to 
designate the field point P by polar coordinates 
(r,0,6), while the location of the surface element 
dy on the dielectric rod is specified by the 
Cartesian coordinates (x,y,z). The relationship 
between these coordinate systems is illustrated 
in Fig. 3. It is easily seen that when r is very 
much larger than the length of the rod, by 
neglecting higher order terms 


r’ =r—x sin@ cosd—y sin8é sing—z cosé. 


(11) 


The variation of r’ over the length of the di- 
electric rod is insignificant in determining the 
value of the integral in Eq. (10), except for 
the exponential term exp(jkr’). Thus r’ is put 
equal to r in the denominator, and this term may 
be removed from the integral together with the 
term exp(jkr) resulting from substitution of Eq. 
(11) for 7’ in the exponential. Hence Eq. (10) 
becomes 

exp( jkr — jut) 
F=jP,(6,o) , 


4nrr 


| 
wy 

(6,0) = [ Ma( cos =~) (sinks) (12) 
> i] 


Xexp(— jk(z cosé+-y sin8 sing 





+x sin@ cos@) )d>. 


This integral is a function of the variables (6,9), 
and of the parameters, Mo, k, k’, and the dimen- 
sions of the dielectric rod ; the integration is to 
be carried out over the surfaces of the dielectric. 

Next, the vector operations of Eqs. (5) and (6) 
must be carried out to find E and H at the field 
point P. Since A is zero, Eqs. (5) and (6) reduce 
to the following : 


= —curlF, (13) 


1 
H = jweF ——— grad div F. 
Jwp 


(14) 


These operations are to be taken with respect to 
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the coordinates (7r,6,¢), and hence P,(@,¢) must 
be explicitly considered. However, if 7 is taken 
large, the terms involving derivatives of P(6,¢) 
become negligible. Hence P, will be taken as a 
scalar constant in the evaluation of Eqs. (13) 
and (14). Neglecting terms in (1/7’), 


1 | jk 
i( -- cost ) 
4x| r 
jk 
-n(= sin@ cose) | 
r 


Xexp( jkr — jwt) P1(0,0). 


curlF = 


(15) 


Now E=-—curlF, which ‘is the negative of Eq. 


_ (15); when kK and i are converted to polar coordi- 


nates, then (the term exp(jkr— jwt)/4ar which 
represents a spherical diverging wave is omitted 
throughout the remainder of this section) 

E,=0, 

E= (jk cosd) P (8,9), 

E,=(— jk sing cos@) P\(0,0). 


(16) 


A complete solution requires the detailed 
evaluation of P;(0,¢). As is shown in Fig. 1, the 
rod is supposed to be of rectangular cross section, 
of length /, average width a and depth 6. Accord- 
ingly, the limits on the integral over Y run from 
0 to / in the gz direction, and from -~—(06/2) to 
+(b/2) in the y direction; but x has only two 
values, +(a/2) and —(a/2). The integrand 
arises from the evaluation of (exp(jkr’)/r’). Since 
x takes on only the values +(a/2), x may be 
replaced immediately by these values, so that 
two terms arise in the integrand. Also, d= may 
be replaced by dS=dydz. Thus: 


a 
P,(6,o) = f af exr( “— sin@ cose ) 


a Ty 
+exp( ik. sin@ cose ) cos— - (sink’z) 
2 b 


Xexp(— jk(z cosé+y sin@ sing))dydz. (17) 
Consolidating terms gives: 
ka 
P\(0,6) = Mo coo sind cose 
wy 
f (cos *) (sink’z) 
s b 
Xexp(— jk(z cosé+y sin@ sing))dydz. (18) 
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Separating the y and z portions of the integral 
then gives: 


ka | 

P\(0,6) = Mo co sind cose | -Ty-1s, | 
| 

| 


+(b/2) Ty 
=f (cos “) 
(b/2) b (19) 


Xexp(— jky sin@ sing)dy, | 
| 
l 
n= f (sink’s) exp(— jkz cos@)dz. 

0 


An evaluation of J, and J. leads to: 


° 


r\? kb 
(*) cos(— siné sing ) 
2b \2 2 


;= Se ee ene ee ere 
°° 


x {T\?* kb 2 
( ) —{ — siné sind ) 
2 2 


1 
I,=—(A —jB), 
2k 


1 L (20) 
A =————[1 —cos}{ (n —cos6)kl} ] 


n—cos@ 


1 
+——[1—cos{(n+cos6)k/} ], 
n-+cos@ 


sin}(m—cos@)kl} sin} (n+cos@)kR/} 


n—cos@ n-+-cos@ 


n=k'/k. 





From Eq. (16) it is seen that for large distances 
from the dielectric rod the electric vector forms 
a plane wave. The magnetic vector may be 
found directly without recourse to Eq. (14), as: 

H,=0, 
k . . 
H,= —-—E, = j— sing cosOP (0,9), 
AYN AYN (21) 
k Rk 
H,=—E4 = j— cos¢P (6,9). 
ANT wu 

The terms in Eq. (16), as amplified by Eqs. 
(19) and (20), have the following explanation : 

(a) The term cos($ka cos@ sin@) in Eq. (19) 
arises because of the fact that the magnetic 
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currents are present on two plane parallel sur- 
faces separated by a distance a. 

(b) The sing and cos¢ terms in Eq. (16) arise 
because the electric field is polarized. 

(c) The term cos@ in E, arises because the 
magnetic currents in a surface element d.S behave 
like dipole oscillators parallel to the y axis. 
(Later it appears that the factor cos@ occurs in 
Es when the source is electric currents on the 
side.) 

(d) The integral J; in Eq. (19) represents the 
part of the pattern contributed by the finite 
width 6 of the dielectric rod. The term cos(zy/b) 
in the integrand of J; represents the shading of 
the magnetic currents in accordance with the 
T Eo,, excitation. 

(e) The integral J; in Eq. (19) represents the 
effect of the length of the dielectric rod. Both 
k and k’ enter into the expression ; the significant 
quantity appears to be k’/k, which depends on 
the thickness of the dielectric rod. 

Experimental patterns are usually confined to 
certain planes: specifically, measurements of the 
E-vector in the planes ¢=0 and ¢=27/2. The 
equations for E simplify: for the E-plane, ¢= 0, 


2b 
Iyw=—, 


T 
I. remains the same, 


ka 2b r ae 
P (0,0) = Mo coo sind] hanes Ts, 


T 





E4= jk: P,(8,0), 
E,=0. | 


Likewise for the /-plane, ¢= 7/2, 


r\? kb 
( ) cos( - sin) 
26 \2 2 


lL (2/2) =— — oe — 


'@O-G 


’ 





(23) 
I, remains the same, 
P(0,9/2) = Mo: Ti¢2/2)* Ie, 
E.=0, 
E,=— jk cosé: P,(0,2/2). 
665 
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Fic. 4. Patterns in E-plane for dielectric rod 
with /=4.22\, n=1.10. 


(b) Electric Currents on the Sides 


A similar theoretical analysis may be carried 
out for the electric currents on the sides of the 
dielectric rod. These currents are given by Eq. 
(1); for the TE»; mode, the field at the surface 
should have a tangential component correspond- 
ing to the component of the field in the z direction 
inside the wave guide: 


HT,=B sin(ry/b) exp(jk’z — jot), (24) 


where B is an arbitrary constant. Evaluating 
the term in y for the sides of the rod, the surface 
electric currents are: 


J =iJ» sink’z exp(— jot). (25) 


In this equation, the term exp(jk’z) has been 
replaced by the term sink’z, indicating that the 
currents in the neighborhood of the opening of 
the wave guide are insignificant in contributing 
to the radiation pattern. This choice of the varia- 
tion with z is based primarily on experimental 
evidence. Making use of Eq. (11) and putting in 
the appropriate values of r’ as before leads to 
the value of the vector potential A: 
exp( jkr — jut) 
=i—- P(0,), 
4rr 
+ (26) 
P.(0,o) -{ J» sink’z exp( — jk(z cos@ 





+y sin@ sing+ x sin@ cos@))d>. | 


Considering only the contributions from these 
currents, and therefore putting F=0, Eq. (6) 
simplifies to: H=curlA. Ignoring the action of 
the vector operations on P2(6,@) in Eq. (26) and 
neglecting terms of the order of 1/r’, the com- 
ponents of H are found to be (omitting the 
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Fic. 5. Patterns in E-plane for dielectric rod 
with /=4.22d, n =1.30. 


factor exp(jkr— jwt)/4ar throughout the rest of 
this section) : 


H,=0, | 


I1,= jk(sing) P2(8,¢), (27) 


IT, = jk(cosé cos) P2(8,¢). | 


These equations are recognized as those for a 
plane wave when r is large, and therefore instead 
of using Eq. (5) for E, more directly: 


E,=0, } 


Eo = (wu/k)H, = (jwp)(cosé@ cosd) P2(8,¢), > (28) 


E,= — (wu/k) Hs = — (jou) (sind) P2(6,¢). 


As before, the integrand in P2(6,@) arises in part 
from the evaluation of (exp(jkr’)/r’). Since here 
y takes on only the values +43), two terms arise 
in the integrand: these two terms become con- 
stants so far as the integration is concerned, and 
the surface element reduces to dxdz. Thus, 


P»(0,6) = —Jo-2 cos{ (kb/2) sin@ sing | | 
Ts Ts) 


+(a/2) 
I,= f exp(— jkx sin8@ cosd)dx, . 
~(a/2) 


| 


(29) 


| 


i= f sink’s exp( — jkz cos@)dz. | 


0 ) 


It is seen that J, is the same as 72 in Eq. (19); 

this integral is evaluated in Eq. (20). J; has the 
value: 

2 sin{ (ka/2) sin@ cos¢ | 

calcd — . 


k sin@ sing 





(30) 


Evaluating the equations for the £- and H-planes 
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of 





gives: E-plane: ¢=0 
2 sin[ (ka/2) siné | 


jg ase 


k sin@ 


I,=I_ remains unchanged, | 
{ 
| 
| 





2 sin[{ (ka/2) sin@ } (31) 
P (0,0) = —2J9— ——_—_——_- 1, 
k sind 
Ee = jwpu(cos@) - P2(6,0), 
Ez =(. 


H-plane: ¢=7/2: 
I 3¢2/2) =a, } 


I, remains unchanged, 


T kb 
p(s, ") = —2J ya coo — sind, - (32) 
2 


Es=0, 
Es = — jopP[6,(4/2) ]. J 





The expressions for the patterns for the E- 
and H-planes for magnetic currents on top and 
bottom and electric currents on the sides as 
summarized in Eqs. (22), (23), (31), and (32) are 
seen to have certain similarities and differences. 
For electric currents the pattern in the E-plane 
indicates a null at 6=2/2 due to the factor cosé. 
This null is absent in the pattern for the mag- 
netic currents. Equation (22), however, indicates 
a null in this latter pattern if ($ka sin@) = 7/2. 
This means sin@=2/ka; but for ka=1.96, the 
experimental value, r/ka>1 and the null is not 
possible. Again, for electric currents the pattern 
in the H-plane indicates a null if (326 sin@) = 2/2, 
or sin@ = x/kb; the experimental value of kb =4.5 
gives +/kb=0.7, whence 6=45°. The pattern 
for magnetic currents does not lead to sucha 
null, but contains a factor of cos@ leading to 
a null at @= 2/2. Other nulls for both cases may 








i \_ _ - _~ 
fa ee wey | 
| A\| 7 x \ | 
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FiG. 6. Experimental pattern in E-plane for dielectric rod 
with /=10.2d, n=1.10. 
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arise for smaller angles nearer 6=0, because of 
the term J; which is common to both pattern 
expressions. Such nulls are found in computed 
values for J2, and are also observed experi- 
mentally. 


(c) Electric Currents on Top and Bottom 


These currents arise from the tangential com- 
ponent of the magnetic vector on the top and 
bottom surfaces. The magnetic vector has the 
components: 


wy 
H,=B, — exp(jk’z— jut), 


my 
H,=B, wal” exp(jk’z— jt). 


These expressions have the same values on the 
two faces, since they do not vary with x; the 
evaluation of the electric currents then gives 
currents which are everywhere parallel, but of 
opposite direction on the bottom compared to 
the top. Hence these currents can produce no 
contribution in the field at large distances in the 
z direction (@=0). Some contribution should be 
expected at moderate angles, but no contribution 
at right angles to the rod (6 = 42). The resultant 
pattern should be somewhat on the order of the 
pattern due to a quadrupole. No detailed evalua- 
tion is given here of this pattern, since it does 
not appear to give contributions near 6=0, the 
direction of primary interest. 


PATTERN COMPUTATIONS 


The fields given by Eqs. (22) and (23) have 
been calculated for specific values of the param- 
eters. These expressions were modified in the 
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Fic. 7. Half-width of major lobe as a function of 
length of dielectric rod. 
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Fic. 8. Angular position of minor lobes as a function 
of length of dielectric rod. 


calculation to find the magnitude of the complex 
expressions involved, and to produce normalized 
values for the fields. The maximum value of the 
magnitude was used as a normalizing factor, 
whether this maximum occurred at @=0 or 
elsewhere. Unfortunately, the expressions for A 
and B cannot be evaluated in terms of a general 
parameter, but must be computed for each choice 
of m and ki. To allow comparison with experi- 
ment, several values for m were chosen, namely, 
1.00, 1.05, 1.10, 1.20, 1.30, and 1.50, for ki = 26.4. 
Typical computed patterns are shown in Fig. 4 
for n=1.10, k1=26.4, and in Fig. 5 for n=1.30, 
kl = 26.4. After comparison with experiment, the 
value of 1.10 was chosen for further computa- 
tions, using a range of values of k&/ from 19.0 
to 64.1. A summary of data for these several 
patterns computed for the E-plane is shown in 
Figs. 7, 8, and 9, together with comparable ex- 
perimental data. The parameter » varies with 
the thickness of the dielectric rod. In very thin 
rods k’~k, and n is very nearly unity. In thick 
rods k’>k, and n becomes nearly equal to the 
refractive index for the dielectric. In experi- 
mental practice it is found desirable to taper the 
rod (see Fig. 1) so that the thickness a decreases 
along the length. This procedure reduces reflec- 
tion effects and improves the radiation pattern. 
For such a situation, k’ would actually be a 
function of z in the equations. This variation is 
allowed for by choosing a mean value of k’ in 
correspondence with experimental evidence. 

In the same way, the fields given by Eqs. (31) 
and (32) were calculated, but for only one or 
two values of the parameters. Because of the 
common term J», essential agreement was ob- 
served between these patterns and those for 
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Fic. 9. Relative height of first minor lobe as a 
function of length of dielectric rod. 


Eqs. (22) and (23), the only differences arising 
in the heights of the minor lobes. These differ- 
ences were not considered significant, and hence 
no further values were computed, it being con- 
cluded that the patterns calculated from the 
effects of the magnetic currents were sufficient 
to describe the fields. 


EXPERIMENT 


Experimental Arrangement 


Experimental patterns for a number of poly- 
styrene rods as terminations for a wave guide 
were obtained by rotation of the rods in a plane 
polarized electromagnetic field. This field was 
produced by a dipole antenna mounted in a 
24-inch parabolic reflector and connected to a 
magnetron operating at X-band frequency. This 
source was placed on the roof of one building 
and directed at another building 400 feet away, 
on the roof of which was mounted the receiving 
equipment. Such an arrangement produced a 
plane polarized field at the receiving position and 
reduced reflections from the ground and other 
surroundings to a minimum. A rotating stand 
coupled to a motor drive was used for rotation 
of the rods in the field. The signal received by a 
rod actuated a suitable receiver whose output 
was recorded on the chart of an Esterline-Angus 
recording milliammeter. A marking relay on this 
recorder allowed interpretation of the record in 
terms of the rotation of the rod on its stand. 

The rods were all of rectangular cross section, 
designed to fit the wave guide snugly near one 
end, and ending within the wave guide in a 
tapered section designed to prevent reflections 
at the interface between dielectric and air. The 
outer section was tapered uniformly to give an 
average thickness of slightly more than half the 
thickness of the wave guide. This section had 
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different lengths for the several rods. A typical 
shape is illustrated in Fig. 1. 


Pattern Measurements 


Patterns in E- and H-planes were taken for 
eight different polystyrene rods of various lengths 
and thicknesses. To obtain a pattern in the 
E-plane, the rod and associated wave guide were 
rotated about an axis perpendicular to the axis 
of the rod and the electric field direction. To 
obtain a pattern in the H-plane, the rod and 
associated wave guide were rotated about an 
axis again perpendicular to the axis of the rod 
but parallel to the electric field direction. In each 
case the orientation of the field from the source 
was such as to produce the 7E»,; mode in the 
receiving wave guide. 

Seven of these rods had constant width equal 
to the inside dimension of the wave guide, but 
were tapered in the E-plane so that the average 
thickness was the same for each. As is shown in 
Fig. 1, the taper does not begin exactly at the 
junction with the wave guide; the taper itself 
may differ somewhat from rod to rod. Hence 
the several average thicknesses are not precisely 
the same. The eighth rod had a length com- 
parable to one of the first seven, but a consider- 
ably different thickness. Typical patterns are 
shown in Figs. 4, 5, and 6. The pattern in Fig. 4 
is for the E-plane for a rod whose length is 
/=4.22d and whose average thickness is 0.1952. 
On this same figure is also shown the theoretical 
pattern. It is observed that agreement is very 
good out to about 45°. The general shapes of the 
curves are the same, the widths of the major 
lobes at the —10-db points are the same, and 
the positions of minima agree well. The heights 
of the first minor lobes are very nearly the same, 
but the second minor lobe is only suggested by a 
slight dip on one side of the experimental curve. 
The theoretical curve shows a minimum about 
45°, but at this same angle it was shown in the 
section on theory that other effects there neg- 
lected might be of some importance. In addition, 
reflections from surrounding objects are more 
likely to affect the low level fields represented 
by minor lobes than the strong signal of the 
major lobe. Hence the deviation of the experi- 
mental curve is not unexpected. The pattern in 
Fig. 5 is for the E-plane for a rod of length 
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1=4.22 and thickness of 0.390. The agree- 
ment between the experimental and theoretical 
patterns is here not as good; the experimental 
curve has some variations near 0° not evident 
from theory, and the positions of the nulls do 
not coincide. A small variation in the choice of n 
in the calculation would improve the agreement 
of the nulls. The pattern in Fig. 6 is for the 
E-plane for a rod of length 10.2A, and average 
thickness 0.195. It is seen that the major lobe 
has coalesced with one or more minor lobes, while 
other minor lobes are considerably smaller. The 
theory developed above indicates for such lengths 
very prominent first minor lobes, but no major 
lobe at all. This discrepancy is attributed to the 
failure of the continuing assumption that the fields 
in the dielectric rod are similar to the fields within 
the wave guide throughout the length of the rod. 

The seven rods of comparable thickness each 
had a portion of uniform rectangular cross section 
about the junction between the rod and wave 
guide; by inserting more or less of the rod into 
the wave guide a certain variation of length was 
available for each rod without appreciable change 
in the average thickness. About 30 different pat- 
terns could thus be obtained for the seven rods, 
with some overlapping of values for length. 
Figure 7 shows the resulting values of half- 
width of the major lobe, illustrating the scatter- 
ing of the experimental points. A smooth curve 
shows the general tendency of these points; this 
curve is not far different from the curve obtained 
from theoretical patterns, for lengths as great as 
about 6d. For rods of longer length, coalescing 
of the major and first minor lobes begins, whereas 
the theory indicates a continuing narrowing of 
the major lobe. A set of data is shown in Fig. 8 
for the angular position of the minor lobes as a 
function of length of rod. The experimental 
points have been omitted here for purposes of 
clarity ; the scattering is much the same as that 
shown in Fig. 7. Agreement with theory is again 
quite good up to /=6,, but is less satisfactory 
thereafter; the minor lobes tend to disappear or 
to coalesce with the major lobe. Agreement be- 
tween theory and experiment for angles greater 
than 45° or so is not good, that is, for further 
minor lobes. For these large angles the theoretical 
assumptions are not expected to hold accurately, 
and reflections from surrounding objects may 
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well upset experimental conditions. Similar data 
for the nulls of the patterns show equally good 
agreement between theory and experiment as 
for the positions of minor lobes. 

The relative height of the first minor lobe ob- 
tained from experiment as a function of length 
is compared to the theoretical curve in Fig. 9. 
Here there is little or no agreement except for 
order of magnitude; the theory indicates in- 
creasing minor lobe height, while experiment 
leads to decreasing height. The heights of minor 
lobes are likely greatly affected by the approxi- 
mations of the theory and by the experimental 
conditions, whereas the positions of lobes and 
nulls are less affected by these factors. Thus the 
good agreement of positions of minor lobes and 
nulls is not accompanied by good agreement in 
minor lobe height. 

A comparable set of data was compiled for the 
H-plane patterns. The same agreements and 
deviations were noted as discussed above for 
the E-plane, without producing any fundamental 
differences between the results for the two planes. 





SUMMARY 


A theory is developed for computing radiation 
patterns for a dielectric rod used as a termination 
for a rectangular wave guide. Subject to some 
restrictions, the theory is applied to the TEo,, 
mode of excitation in the wave guide, and 
specific expressions are worked out for the radia- 
tion patterns. The expressions contain as param- 
eters the length, width, and thickness of the 
dielectric rod, and its apparent index of refrac- 
tion. Suitable choice of the index of refraction 
produces computed patterns agreeing in con- 
siderable detail with patterns for polystyrene 
rods whose lengths vary from 3d to 6, and whose 
average thicknesses and widths are the same. 
Also, a comparable computed pattern is obtained 
for a considerably thicker rod in the same length 
range. 

It is concluded that, given the physical param- 
eters, the analysis may be used to compute 
radiation patterns for the TE ,, excitation of 
rectangular dielectric rods. 





The Attenuation of Spherical Shock Waves in Air* 


RicHarD G. STONER AND WALKER BLEAKNEY 
Palmer Physical Laboratory; Princeton University, Princeton, New Jersey 


(Received February 13, 1948) 


The peak pressure of spherical shock waves from small 
explosive charges has been determined as a function of 
distance by measuring the velocity of propagation and 
applying the velocity-pressure relation derived from the 
Rankine-Hugoniot equations. The charges were composed 
of TNT or pentolite (TNT/PETN 50/50 by weight), 
weighed from 1.45 to 8 lb., and were spherical or cylindrical 
in shape, except for one series of small rectangular blocks. 
Measurements on the non-spherical charges were made in 
the plane through the center of the charge, perpendicular to 
the axis. 

The pressure-distance relations for the four principle 
charge types are given by the following fitted equations, in 
which II represents excess peak pressure in atmospheres, 
' and the distance, scaled according to charge weight, is 
given by the non-dimensional variable Z = R/(pr)! where R 
is the distance from a charge of volume r and specific 
gravity p: 


I. INTRODUCTION 


T has long been recognized that pressure 
waves of finite amplitude in mediums like air 


* This work was supported in part by the Office of 
Naval Research, Contract No. N6ori-105 Task II. 
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}-lb. rectangular blocks, TNT, 

Ig = 13.50/Z —769.9/Z* +-36280/Z?; 
8-lb. cylinders, Pentolite, 

I4 = 10.49/Z — 135.5/Z?+21070/Z3; 
4-lb. cylinders, TNT, 

IIg = 11.34/Z —185.9/Z*+19210/Z3; 
3.75-lb. spheres, Pentolite, 

Ie =8.63/Z +295.1/Z?+7823/Z'. 

These equations are valid for values of Z between approxi- 
mately 18 and 110. The indicated probable error of the 
fitted curves is of the order of one percent for intermediate 
distances, increasing to from two to seven percent at the 
extremes of the distance range covered. 

The curve for spherical charges is in agreement with the 
Kirkwood theory, and the results for cylinders having vari- 
ous length/diameter ratios indicate large dependence of 
pressure on charge shape for values of Z up to 50 at least. 


propagate at greater than sound velocities, the 
velocity increasing with amplitude. For this 
reason, a continuous wave of finite amplitude 
changes its form with time, a pressure front be- 
coming steeper until limited by the viscosity 
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and heat conductivity of the medium. When 
this relatively stable form is reached, the pressure 
rise occurs within an extremely narrow region, 
and for most purposes may be considered a 
perfect discontinuity. Such a disturbance is called 
a shock wave. The conditions governing the 
propagation of such discontinuities were formu- 
lated by Rankine and Hugoniot, who gave ex- 
pressions for the conservation of mass, mo- 
mentum, and energy across the shock front." 

The spherical wave produced by the detona- 
tion in free air of an explosive charge is a wave of 
this type, having the characteristic form shown 
in Fig. 1 for a pressure jump of the order of 0.5 
atmosphere.? The initial rapid rise to a peak 
pressure P is followed by a more gradual decay, 
a shallow suction phase, and a return to atmos- 
pheric pressure Po, either gradual, or including a 
small secondary shock. The peak pressure decays 
as the wave propagates outward for two reasons: 
(1) because the energy must be distributed over a 
constantly increasing area, and (2) because there 
is a finite dissipation of energy in the transition 
through the shock front. Exact solution of the 
hydrodynamical equations, subject to initial con- 
ditions describing the release of energy in the 
explosive process and to the Rankine-Hugoniot 
conditions at the shock front, is not possible be- 
cause of the non-linear nature of the equations 
of motion, and direct numerical integration is, 
in general, prohibitively laborious. In a recently 
published paper,’ Kirkwood and Brinkley, by 
imposing a similarity restraint on the energy- 
time curve of the shock and using the second 
law of thermodynamics to determine, at an 
arbitrary distance, the distribution of the initial 
energy, have reduced the problem to two ordi- 
nary differential equations which are amenable 
to numerical integration. Their computations 
then provide a theoretical peak pressure vs. 
distance relation. 

Experimental determination of the peak pres- 
sure of shock waves has also proved to be diffi- 
cult. The most common method and that which 


'W. J. M. Rankine, Phil. Trans. A160, 277 (1870). 
H. Hugoniot, J. de l’ecole polyt. 57, 3 (1887) ; 58, 1 (1888). 

? The form of very weak shock waves is somewhat 
Starens, and has been investigated by J. W. M. DuMond, 
E. R. Cohen, W. K. H. Panofsky, and E. Deeds, J. Acous. 
Soc. Am. 18, 97 (1946). 

7S. R. Brinkley, Jr.and J. G. Kirkwood, Phys. Rev. 71, 
606 (1947). 
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was for a long time considered the most promising, 
has been the use of piezoelectric crystals, which 
produce a charge proportional to the applied 
pressure. The fundamental objection to such 
gauges, which has only recently become clear, 
results from the disturbance of the flow caused 
by the presence of the gauge itself. This effect 
makes the apparent pressure as read by the gauge 
less than the true pressure (i.e., the pressure that 
would exist in the absence of the gauge). The true 
pressure, moreover, is not a linear function of the 
gauge reading. Calibration of the gauges there- 
fore becomes a difficult matter. These difficulties 
can be reduced considerably by suitable stream- 
lining of the gauge and mount, but a complete 
solution requires a method of dynamic calibra- 
tion exactly duplicating the conditions under 
which the gauge is used, which has not so far 
been achieved. 

A technique for pressure measurements which 
is not subject to these objections has been sug- 
gested in the application of the relation between 
the peak pressure and the velocity of propagation 
of a shock. If gauges are used simply to time the 
passage of the shock over a known distance, the 
need for calibration is eliminated, and the local 
effect due to the presence of the gauges is un- 
important so long as the over-all speed of propa- 
gation of the wave is not reduced. While the 
experiment requires careful planning to obtain 
satisfactory precision, the elimination of sys- 
tematic errors should assure more reliable results 
than have been obtained previously. 

The derivation of the pressure-velocity relation 
has been given in numerous places.‘ Using the 
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Fic. 1. Typical shock wave from a small charge. (8 feet 
from 0.5 lb. of TNT. P—Po~0.5 atmosphere.) 


*G. I. Taylor and J. W. Maccoll, Thery ede of — 
pressible Fluids,’ from Aerodynamic T: te yy 
Durand (Julius Springer Verlag, Berlin” 1935) ol. iI 
Div. H; H. W. Liepmann and E. E. Puckett, Aerodynamics 


671 














e———1.—__—er" + 





CaP l -—4f— -—- —--}-—_ - +44 x 
° initia —i ~aQwo< 0---0-- O--- 
/ a ; 
a oe CHARGE GAUGE LINE 

POSITIONS 

ales SCALE 

N some 
Qo 10 20 
FEET 


Fic. 2. Plan of experimental set-up. Ten gauges were mounted uniformly, 18 inches apart, along each gauge line. 


Rankine-Hugoniot relations and calculating the 
internal energy of the gas from the ideal gas law, 
we obtain 


(P— Po) Py=[2y yt DIL(V? a*) —1], (1) 


where P is the absolute pressure immediately 
behind the shock front, Po» is the atmospheric 
pressure ahead of the shock, V is the velocity of 
the shock relative to the medium ahead, a is the 
velocity of sound in this medium, and y is the 
ratio of specific heats Cp/Cy. At high pressures 
the ideal gas approximation with constant 
in the energy equation is not correct. Comparison 
of Eq. (1), using y=1.40 for air, with more 
exact calculations® indicates, however, that for 
an excess pressure P — Py of 100 lb. /in.? the error 
in P—P, is less than 0.7 percent and that this 
error diminishes rapidly for lower pressures. 
Examination of Eq. (1) shows that the ratio 
of the percent error in excess pressure P — Po to 
the percent uncertainty in the velocity ratio V/a 
is never less than two, and is very large for pres- 
sures only slightly greater than Po. In the pressure 
region of interest in this investigation (~0.1 to 


10 atmos. excess) the error ratio is about ten in 


order of magnitude. Thus in order to determine 
the excess pressure within 1 percent, velocity 
measurements must be good to about 0.1 percent. 
This was taken as the design criterion, and the 
experiment was planned so that this precision 
could be expected in all the measured quantities 
involved in determining the velocity ratio. 
Il. EXPERIMENTAL PROCEDURE 

To measure shock wave velocity, quartz piezo- 
electric gauges were chosen as being best suited 
of a Compressible Fluid (John Wiley and Sons, Inc., New 
York, 1947). 


5 |. G. Kirkwood, S. R. Brinkley, Jr. and J. M. Richard- 
son, unpublished work. 
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to the requirements of the problem. These con- 
sist of an x-cut quartz disk, 0.9 inch in diameter 
and 0.1 inch thick, mounted in a sealed brass 
housing with a thin dural face plate which com- 
municates the external pressure to the x face of 
the crystal. Used “‘face-on”’ (the blast incident 
normally on the sensitive face) these gauges 
provide an extremely fast response to pressure 
changes, and the flat face furnishes a definite 
point from which to measure the baseline dis- 
tance. Ten such gauges were mounted rigidly at 
intervals of 18 inches and connected in parallel 
on a single channel of a multichannel cathode-ray 
oscilloscope unit. This was part of a mobile unit 
developed especially for blast recording, and 
provided amplifiers with good high frequency 
response and exceptionally well focused traces so 
that the time of arrival of a shock wave at a gauge 
could be determined within about a micro-second. 

Two lines of ten gauges each were set up on 
opposite sides of the charge, with a distance of 40 
feet between the first gauge of each line, as shown 
in Fig. 2. With a spacing of 18 inches between 
gauges, a total distance range of about 7 feet to 
47 feet could be covered by moving the charge 
position back and forth between the two gauge 
lines, 18 pressure points being obtained from the 
two gauge lines on each shot. Charge positions 
were located at the point midway between the 
gauge lines, and symmetrically on each side of 
this point. In this way the two independent sets 
of gauges could be checked against each other 
directly (charge in center position) and in effect 
interchanged in covering distances close to and 
far from the charge by taking shots at corre- 
sponding positions on opposite sides of the center. 
Excellent agreement was observed between the 
results of the two groups of gauges. 
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To reduce the risk of the nearer gauges blanket- 
ing those behind, each gauge line was rotated at 
an angle of about 9° from the line joining the 
farthest gauge of each line. A test showed that 
with the gauges staggered to this extent there 
was no observable attenuation of the shock wave 
due to its passage along the gauge line. Careful 
measurements of the angle and distance of each 
gauge line with respect to each charge position 
were made, as well as of the intergauge distances 
along each gauge line, so that the differences in 
the radial distances of successive gauges from the 
charge could be calculated. These represent the 
true baselines over which the spherical shocks 
were timed. When the charge position was 
changed each gauge was rotated about a vertical 
line through the center of its sensitive face so as 
to remain normal to the radial line from the 
charge. 

Short pulses of about five microsecond dura- 
tion were superimposed at a rate of 10,000 per 
second on each recording channel to facilitate 
timing the intervals between the arrival of the 
shock at successive gauges. These pulses were 
produced by a crystal controlled oscillator 
checked against the standard carrier frequency 
of station WWV. 

The charges used were either spherical or 
cylindrical in shape, weighing between 3 and 8 
pounds. The cylinders were hung with their axes 
vertical, so that all the gauges were in the 


equatorial plane of the charge. In measuring 
“free air’’ peak pressure, the charge and gauges 
must be sufficiently far above the ground to 
eliminate interference with the incident wave 
arising from reflection phenomena at the ground. 
In the present experiments the gauges and 
charges were located in a horizontal plane 10 
feet above the ground, a height sufficient to 
guarantee that the effect of the ground could 
not be felt until the shock wave had passed the 
farthest gauge. 

The velocity of sound in the medium must 
also be known in applying the pressure-velocity 
relation. Because of the difficulty of controlling 
or measuring random local variations in the 
temperature, humidity, local velocity, etc., of 
the outdoor atmosphere, which affect the velocity 
of both sound and shock waves, the velocity of 
sound was measured directly and as nearly as 
possible at the same time and place as were 
the shock waves themselves. Ordinary telephone 
receivers were set up as indicated in Fig. 2 to 
provide two independent baselines of 30 and 42 
feet over which the velocity of the wave pro- 
duced by the detonation of a blasting cap could 
be measured. The caps were located at a suffi- 
cient distance for the velocity of the pressure 
wave to have decayed to sound velocity within 
the limits of observation. Two caps were deto- 
nated in close succession at positions N and S 
(Fig. 2) to provide determinations of sound 


TABLE I. Charges used in experimental program.* 








’ Number 

Average» Diametere Height of pres- 

Desig- weight Height* or width diameter Number _ sure ob- 

nation Composition (pounds) Shape (inches) (inches) ratio of shots servations 
Q Pressed TNT 0.5 Rectangular block 3.50 1.75 — 17 117 
A Cast Pentolite* . 8.00 Cylinder 8.75 4.50 1.94 5 88 
B Cast TNT® 4.108 Cylinder 7 3.62 1.93 8 144 
Cc Cast Pentolite* 3.75 Sphere — 5 -- 8 142 
D Cast TNTe 2.89% Cylinder 5 3.62 1.38 S 90 
E Cast TNT® 2.89" Cylinder 12 351 5.20 5 90 
F Cast TNT® 1.45' Cylinder 2.5 3.62 0.69 1 18 
G Cast TNT® 1.45) Cylinder 6 2.31 2.60 3 52 
H Cast TNT® 1.45! Cylinder 8 2.00 4.0 3 45 
K Cast TNT® 1.45! Cylinder 12 1.62 7.4 1 18 








® The spherical pentolite charges were prepared by the Woods Hole group of Division 2, NDRC. All cylindrical charges were supplied by the 


Picatinny Arsenal. 


> The weights of individual charges varied by as much as +5 percent in some cases from the average values listed. Accurately measured weights 


were used in scaling to reduce results to a common basis. 


* Dimensions given are nominal values. Variations of the order of +5 percent occur in some cases. 
450 percent TNT and 50 percent PETN, by weight. Specific gravity ~1.60. 


* Specific gravity of cast TNT ~1.55. 

‘ Includes 0.031-Ib. tetryl pellet booster. 
* Includes 0.096-Ib. tetry! pellet booster. 
» Includes 0.067-Ib. tetryl pellet booster. 
‘ Includes 0.030-lb. tetryl pellet booster. 
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velocity in opposite directions, from which the 
velocity of sound relative to still air and the 
component of wind velocity in the direction 
parallel to the gauge lines could be determined. 
The effect of cross wind velocity (although it 
was not measured) is entirely negligible since it 
is a second-order effect, and, by good fortune, 
the prevailing wind direction during the period 
of shooting was very nearly parallel to the gauge 
lines, so that the cross-wind component was in 


5, 





any case only a fraction of the parallel wind 
component which was measured. The two caps 
and the main charge were set off automatically 
at intervals of about 0.1 second to eliminate as 
far as possible the danger of fluctuations in the 
atmosphere between the time of the sound 
measurements and the shock measurements. 
The proper wind velocity correction was added 
to the shock velocity measured over each in- 
terval of the gauge line, and the corrected shock 
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Fic. 3. Experimental pressure vs. distance results for four types of charges. Excess shock wave pressure in non-dimen- 
sional units, I1=(P—Pp»)/Po, is plotted against the non-dimensional distance Z= R/(pr)!. For specific gravity p= 1.60, 
Z=1.06R/D, where D is the diameter of a sphere having the same volume as the charge. Thus Z represents approximately 
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the distance in charge diameters, for the explosives used. (For TNT, p=1.55; for pentolite, p= 1.60.) 
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TABLE II. Fitted curves giving best values of experimental pressure measurements on four series of charges. 








Distance range covered 
by experimental data 


Fitted curve for pressure 


Standard deviation* of 1 
(expressed in percent) 





Type of charge [Z =R/(pr)}] =(P—Po)/Po o1 min o 
TNT Eng’r Blocks (Q) 21 <Z<107 Ilg= 13.50/Z —769.9/Z?+-36280/Z3 +3.9% +2.1% +4.9% 
Pentolite Cylinders (A) 14.7<Z< 91 Il4 = 10.49/Z — 135.5/Z?+21070/Z3 7 0.9 1.4 
TNT Cylinders (B) 18.6<Z<i1ll1 IIp = 11.34/Z —185.9/Z?+19210/Z3 6.5 0.8 aon 
Pentolite Spheres (C) 18.6<Z<119 Ie= 8.63/Z2+295.1/Z?+ 7823/23 2.1 0.6 1.4 











* o1 is the standard deviation at the smallest Z indicated in column 2 and o2, at the largest Z. omin represents the minimum value of standard 
deviation reached usually at about the midpoint of the Z range. The curve for @ is generally flat in the neighborhood of omin, and the standard 


deviation is roughly <1.Semin in the region 40 <Z <90. 


velocity V was inserted in Eq. (1), along with 
the measured sound velocity a, and the corre- 
sponding pressure calculated. Since the shock 
velocity changes appreciably with distance, the 
pressure calculated from the measured average 
velocity over the interval corresponds to the 
instantaneous pressure at a point within the 
interval determined by the relation 


Rm+ 
V(Ry) =(1/Rag) J V(R)\dR, (2) 
R 


m— 


where Ry is the distance from the charge at 
which the instantaneous pressure equals the 
pressure computed from the average velocity, 
R,, is the distance from the charge to the mid- 
point of the interval over which the velocity is 
measured, g is the ratio of the length of the inter- 
val to Rn, Rm+=Rn(1+q/2), Rn—=Rmn(1—q/2), 
and V is the function expressing shock velocity 
in terms of distance from the charge. 

An approximate expression for Ry may be 
obtained by assuming that the pressure over a 
short interval may be represented by an ex- 
pression of the form P—P,»=AR", where A 
and are constants, and expanding the resulting 
expression for V(R) obtained by substituting 
this in Eq. (1). The expression valid for low 
pressures was first given by Carr, Schwarzschild, 


and Weiss, working at Aberdeen, and has the 
form 


Ry =Rni1—(14+n/2)q?/24}. (3) 


A slightly different expression is obtained when 
an expansion valid at high pressures is used. 
A more detailed derivation of Eq. (3) shows 
that it actually holds for moderately high pres- 
sures, and that in the most unfavorable case in 
the experiments reported here the error in Ry is 
only 0.18 percent. 
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III. RESULTS 


The physical characteristics of the types of 
charges investigated are given in Table I. In 
order to compare results from charges of different 
weights, the data have been reduced to a common 
basis by scaling the distances according to the 
principle of similitude. For similar charges of 
equal density, equal pressures are expected at 
distances proportional to the linear dimensions 
of the charges, and hence to the cube root of 
the charge weight. This scaling law has been 
shown experimentally to be valid over a large 
range of charge weights. The scaled distances are 
given in terms of the non-dimensional variable 
Z=R/(pr)', where R is the actual distance from 
the charge, having a volume 7 and specific 
gravity p relative to water. The pressure is 
expressed non-dimensionally in units (P — P»)/Po, 
denoted hereafter by II. 

The charges designated A, B, and C in Table I 
formed the major part of the experimental pro- 
gram, and the complete sets of pressure-distance 
points for these three series and for the pre- 





®*It has been customary in the past to represent the 
scaled distance by the dimensional variable z = R/w! where 
R is the distance in feet, and w is the charge weight in 
pounds. The quantity z thus has dimensions ft./(lb.)!. 
The relation between z and the non-dimensional variable 
Z=R/(pr)', which we have introduced, is obtained by 
replacing w by 62.43pr. It follows that Z=(62.43)!z 
= 3.9672. Thus a distance in Z units when divided by 4 
represents roughly the distance in feet from a one pound 
charge. 

The convenient relationship between the Z units and 
charge diameters may be seen by substituting V=2xD*/6 
where D is charge diameter. We then have Z=(xp/6)~! 
<(R/D). Thus Z represents exactly the number of charge 
diameters (R/D) when (xp/6)~ = 1; i.e., for p=6/r = 1.910. 
Because of the cube root involved, Z also represents 
approximately the distance in charge diameters for a 
considerable range of values of p. For p= 1.50, for example, 
(rp/6)~' = 1.084, and Z differs from the number of charge 
diameters by 8.4 percent. 

The “density of loading,’’ as the term is customarily 
used, is the same as the specific gravity p used here. Hence 
R will be in feet if 7 is in cu. ft., in cm if 7 is in cc, ete. 
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liminary series (Q) using U.S. Corps of Engineers 
Demolition Blocks of pressed TNT are shown as 
small dots in Fig. 3, A, B, C, and Q, respectively. 

The curves shown represent least square fits 
of an expression of the form Il=AZ“'+BZ~ 
+CZ-*, introduced by Reynolds,’ which has been 
found to give a satisfactory representation of the 
data over the range of Z involved. Asymptotic 
solutions of the propagation equations show a 
decay in pressure proportional at large distances 
to Z~' except for a slowly varying logarithmic 
factor, and experimental measurements in general 
point to Z~* as a possible limiting value for the 
dependency at small distances. The Z~ term 
provides a third parameter necessary for fitting 
the data successfully. Since the pressure in the 
explosive is known to be finite, an expression of 
this form cannot hold for very small Z. For this 
reason the fitted curves should be extrapolated 
only with great caution beyond the range of the 
experimental data. The equations of the fitted 
curves for the first four groups of charges in 
Table | are given in Table I 1. The limiting values 
of Z indicate the range covered by the experi- 
mental data. The quantity o is an estimate of 
the reliability of the fitted curve, expressed as a 
standard deviation; that is, if the experiment 





were repeated, in 2 of the cases the value of Il 
I 3 


given by the fitted curve would lie within +o of 
the “‘true’”’ curve which would be obtained from 
a very large number of observations in the same 
kind of experiment. The values of o given in 
Table Il are expressed as a percentage of the 
value of II as given by the fitted curve, which 
represents the best value of the pressure as 
indicated by the experimental data. The stand- 
ard deviation of an individual point is ~5 times 
the o for the curve. The statistical treatment by 
which these o’s were calculated took into account 
the non-independence (in the statistical sense) of 
data from adjacent gauge line intervals.* 

The relationship of the various experimental 
determinations to each other and to the theory 
are shown in Fig. 4. To emphasize the comparison 
of the various curves, the pressure curve for 
pentolite spheres (II, in Table Il) has been 
chosen arbitrarily as standard, and the differ- 
ence of each of the other curves from it, ex- 
pressed as a percentage of excess pressure, is 
plotted as ordinate. The ordinate is thus the quan- 
tity 100(1I—TII,)/1, (equivalent to 100(P —P,) 
(P.—P»)). The abscissa is the non-dimensional 
scaled distance Z. Thus the curve for pentolite 
spheres is represented by the horizontal axis. 











Fic. 4. Comparison of results. 
Each curve represents the per- 
cent deviation of excess pressure, 
100(1l—Ic)/Ic, from the value 
for pentolite spheres IIc, plotted 
as a function of the non-dimen- 





sional distance Z. The experi- 








mental curves are designated by 















letters as in Table I. Curve Th 
represents Kirkwood and Brink- 
ley’s theoretical prediction for 
TNT spheres. 








Zs » & 
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7G. T. Reynolds, “Small charge air blast experiments,’"” OSRD Report No. 1518 (1943). 
® We are indebted to Dr. H. Sheffe for developing a form of statistical treatment appropriate to this problem. 
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(The choice of this curve as standard is only for 
convenience in making the greatest number of 
meaningful comparisons and, of course, implies 
no fundamental distinction in the reliability of 
this curve relative to the others.) The shaded 
region indicates the standard deviation band +o 
(expressed in percent). Bands of similar width 
apply to the other experimental curves. The 
various curves are identified according to the 
designations used in Table I. 

Curve Th represents the theoretical curve for 
spherical charges of TNT computed by Kirk- 
wood and Brinkley on the basis of their recently 
reported theory.’ A curve for pentolite, which 
could be compared directly with the experimental 
results, is not available, but experimental work 
at Woods Hole indicates that 5 percent probably 
represents the superiority in pressure level of 
pentolite over TNT, averaged over the whole 
distance range involved. Thus the approximate 
constancy of the difference between curve Th 
and curve C, with an average value of from five 
to seven percent, represents good agreement 
between the theory and the experimental data. 

The effect of charge shape may be seen directly 
from curves C and A which represent pentolite 
spheres and cylinders, respectively. Curve B is 
for TNT cylinders having the same length, 
diameter ratio as the pentolite cylinders. Its 
approximate similarity to curve A is a further 
confirmation that the observed difference be- 
tween cylinders and spheres is real, but the two 
curves raise some question whether (at least in 
the case of cylindrical charges) the ratio of the 
pressures from TNT and pentolite is a constant, 
independent of distance. The shape of curves A 
and B may be explained qualitatively by the 
observation that the shock wave from the sides 
of a cylindrical charge should at first be more 
nearly cylindrical than spherical, and thus show 
a slower decay and relatively higher pressure at 
small distances from the charge. At some point, 
the rarefaction propagating in from the ends 
would begin to take effect, and the wave could 
then be expected to decay more rapidly than the 
corresponding spherical wave. It is somewhat 
surprising that this effect is apparently so pro- 
nounced at distances greater than 20 charge 
diameters, but in view of the indicated precision 
of the experimental determinations, the observed 
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effect must be considered real. The explanation 
suggested is of course only qualitative. 

The remaining types of charges in Table I, 
consisting of smaller cylinders of TNT having a 
larger range of length/diameter ratios, were not 
used in sufficient numbers to make it advisable 
to treat the data statistically, but curves could 
be fitted by eye with reasonable assurance. Be- 
cause of the smaller size of the charges and the 
smaller number of observations, the results are 
less reliable, but are almost certainly good to 
better than 10 percent, except at the extremes of 
the distance range. Curves D and E are for 
2.9-lb. TNT cylinders with length/diameter 
ratios of 1.38 and 5.20, respectively. The former 
is somewhat less than the ratio for the charges 
of curve B, the latter considerably larger. The 
relative values of the three curves at small 
distances are consistent with the argument sug- 
gested above, and the three curves approach 
agreement at large distances, which is to be 
expected. The data from 1.45-lb. cylindrical TNT 
charges with four different length /diameter ratios 
varying from 0.69 to 7.4 also showed similar 
behavior, but the number of shots was so small 
(only 2 shots, providing 18 pressure points each, 
with no overlap in distance, in the case of the 
smallest and the largest L/D ratios) that an 
accurate comparison is of little significance, and 
the curves have not been included in Fig. 4. 

An extreme example of the effect of charge 
shape is seen in curve Q, representing measure- 
ments made with standard Corps of Engineers 
Demolition Blocks of pressed TNT, whose meas- 
urements are given in Table I. These blocks have 
a square cross section, with slightly rounded 
corners, and measurements were made in the 
equatorial plane in a direction normal to a flat 


face. Thus the effect of both corners and ends 


may be expected, and a large divergence from 
the pressure-distance law for spherical charges 
is observed, considerably larger than in the case 
of even the most elongated cylinders investigated. 

There have been published, as far as the 
authors are aware, no experimental results with 
which the data in this paper can be compared. 
They are aware, however, of the existence of 
other results to which reference cannot be made 
at this time because they have not been pub- 
lished. 








IV. CONCLUSION 


The velocity method of measuring shock wave 
pressures which has been developed appears to 
offer a means of determining absolute peak pres- 
sure which is (1) free from the systematic error 
inherent in direct gauge measurements due to 
the physical presence of the gauge in the flow, 
and (2) capable of precision of better than i 
percent in the excess pressure when a sufficiently 
large number of measurements is made. A funda- 
mental limitation appears to be the nonuni- 
formity of the atmosphere under even the best 





field conditions. On the basis of such measure- 
ments, Kirkwood’s theory has been checked 
experimentally, and a difference in the decay of 


‘pressure with distance arising from the shape of 


the charge has been observed which persists at 
relatively large distances from the charge. 

The authors wish to take this opportunity to 
express their appreciation to Professor A. H. 
Taub for many helpful consultations during the 
early development of the program, and to Messrs. 
L. Christiansen, H. Krieg, R. W. Symons, and 
R. W. Weston, who constituted the field crew. 





Chromatic Aberration and Resolving Power in Electron Microscopy 
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Contrast difficulties in the observation of fine detail in thin specimens may be overcome by 
employing dark-field observation. This introduces the added complication, however, that the 
imaging electrons are now inhomogeneous, having experienced a most probable energy loss of 
the order of 20 volts. The present calculations of the effect of chromatic aberration and diffrac- 
tion show that with an objective aperture smaller than the “optimum aperture,”’ and centered 
optics, such inhomogeneities do not appreciably affect the resolution, even in the extreme case 
that the electrons are uniformly distributed in energy. 


HE recognition of fine detail attainable in 

electron micrographs depends on the ulti- 
mate resolving power of the instrument, as 
determined by lens aberrations! and diffraction 
effects, on the one hand, and by the image con- 
trast under the given operating conditions on 
the other. The condition generally regarded as 
favorable for attaining the highest resolutions— 
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Fic. 1. Arrangement of illumination for bright-field 
and dark-field observation. 








' The effects of field asymmetries as well as voltage 
fluctuations (effective because of the chromatic aberrations 
of the lenses) are included herein. 
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namely, bright-field observation of very thin ob- 
jects with a physical objective aperture at least 
equal to the “optimum aperture’’—inherently 
leads, unfortunately, to relatively low contrasts. 
Apart from contrast limitations, the resolution 
is here limited by spherical aberration and 
diffraction, leading, according to Glaser’s wave- 
mechanical theory,” to the expression 


dinin =0.56(Cf/X)! (1) 


‘ 


for an “optimum aperture’”’ 


opt = 1.13(Cf/d)-*. (2) 


Here dmin is the least resolvable separation of 
two object points, C the dimensionless spherical 
aberration constant of the objective, f its focal 
length, and \ the de Broglie wave-length of the 
imaging electrons. Equation (1) represents a 
lower limit to the resolution obtained with the 








*W. Glaser, “Image formation and resolving power of 
the electron microscope from the standpoint of wave 
mechanics,” Zeits. f. Physik 121, 647 (1943). 
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electron microscope by any means. As already 
indicated, unfavorable contrast conditions pre- 
vent a close approach to this limit for the 
majority of specimens. 

An obvious remedy for this contrast difficulty 
consists in the employment of dark-field obser- 
vation ; here the background brightness is greatly 
reduced since only radiation scattered by the 
specimen is utilized for the formation of the 
image. If the illuminating beam is made to im- 
pinge on the specimen at an angle and the 
objective aperture is centered as for bright-field 
work. (Fig. 1), the geometric aberrations and 
diffraction effects may, in principle, be made 
approximately the same as for bright field. 
However, most of the scattered electrons, unlike 
those which pass through the specimen without 
deflection, experience an appreciable loss in 
kinetic energy. These energy losses lead to 
enhanced chromatic aberration effects in the 
image which, in the past, have been considered 
responsible for the reduction in resolving power 
always observed when an electron microscope is 
used as a dark-field instrument. Since this in- 
crease in chromatic aberration appears to be an 
essential concomitant in the passage from bright 
field to dark field it becomes of interest to in- 
vestigate the velocity distribution in the scattered 
electrons used for dark field imaging and to 
examine to what extent chromatic aberration 
reduces resolving power. 

Figure 2 shows the energy spectra of 50-kv 
electrons which have been scattered by a 100- 
A.U. collodion foil through the small angles 
indicated. These curves, obtained with the micro- 
analyzer arranged as shown in Fig. 3, have been 
corrected for the sensitometric properties of the 
photographic plate, but not for the finite slit 
width of the instrument. Hence they represent 
simply the electron distribution obtained on the 
analyzer plate. The geometry of the apparatus 
was such that the electrons passing through any 
portion of the entrance slit had been scattered 
by the object over an angular range 0.0008 
radian in width. With the slit adjusted so as to 
coincide with the electron probe (Fig. 4) the 
angular deflection of the electrons could be 
computed from the displacement of a particular 
point of the spectrum from the overexposed 
image of the probe; it is obtained simply by 
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Fic. 2. Electron distributions in velocity spectrograms 
formed by electrons scattered by a thin collodion film 
through the indicated small angles. 


dividing this displacement by the length of the 
electron path from the object to the plate. Quite 
similar results were obtained, of course, from 
spectra obtained by displacing the slit; here the 
angular deflection is given by the ratio of the slit 
displacement to the object-slit separation. 

The curves in Fig. 2 show a predominance of 
the inelastic scattering at very small angles, of 
elastic scattering (zero energy loss), at greater 
angles. This is in accord with theoretical expecta- 
tion.? The most probable energy loss of the 
inelastically scattered electrons is seen to be of 
the order of 20 electron volts, which presumably 
represents an averaged energy of excitation from 
the uppermost occupied to the nearest unoccupied 
energy bands of the nitrocellulose lattice. 

In order to compute the effect of chromatic 
aberration resulting from such velocity distribu- 
tions the simplifying assumption was made that 
the physical aperture of the objective corre- 
sponded to an angular aperture appreciably 
smaller than the optimum aperture given by 


3 See, e.g., L. Marton and L. |. Schiff, ‘‘Determination 
of object thickness in electron microscope,”’ J. App. Phys. 
12, 759 (1941). 
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Fic. 3. Arrangement of 
microanalyzer for the deter- 
mination of velocity distri- 
bution of electrons scattered 
through a definite angle. 





Eq. (2). Under these circumstances the effect of 
spherical aberration is negligible, so that only 
chromatic aberration and diffraction need be 
considered in determining the fundamental re- 
solving power of the microscope. Since prelimi- 
nary calculations for the velocity distributions in 
question indicated that the effect of chromatic 
aberration, also, was very small, it seemed of 
interest to treat the most extreme case of chro- 
matic aberration—namely, that obtained with 
uniform energy distribution of the electrons. 

If fo is the focal length of the objective for 
Vo-volt electrons, the displacement of the point 
of convergence of the wave in image space corre- 
sponding to a change eAV of the kinetic energy 
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Fic. 4. Velocity spectrogram of 50-kv electrons which have 
passed through a thin collodion film. 
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of the electrons is 


2>= M?K(A V/ Vo) fo. 


(3) 


Here M is the magnification and K a factor of 
the order of unity. Hence, if u(z,r) is the ampli- 
tude of the wave at an axial separation z and a 
radial separation r from the point of convergence, 
the intensity (or concentration of electrons) at a 
point a distance r from the axis is, for uniform 
energy distribution of the electrons, given ap- 
proximately by 


1= f u(z,r) | ds. 


—2 


(4) 


This relation implies the validity of Eq. (3) 
throughout the range of kinetic energies, and neg- 
lects the effect of the variation of the de Broglie 
wave-length A with kinetic energy as well as the 
variation in the solid angle intercepted by the 
objective aperture at the point of convergence. 
It applies, thus, most accurately for small values 
of r, where electrons with small values of AV 
contribute most to the value of the integral. 

The value of u(z,r) for a converging spherical 
wave limited by a distant diaphragm (the ob- 
jective aperture) to a cone of half-angle a is 
given, according to Debye,’ by 


tk a 2r 
u(z,r) =— f f exp(ik(z cosdo 
2n 0 0 


+r sind» COS¢go)) sind odd od go 


-ik f exp(tkz costo) 


<X Jo(kr sind) sindoddy. (5) 


Here k=27/X and a, the half-angle of the wave 
cone, will be regarded as constants, corresponding 
to the values for the wave with its convergence 
Jo(x) is the Bessel function of 
order zero. The intensity expression in Eq. (4) 


point at z=0. 


may then be evaluated with the aid of a series 
of transformations and a change in the order of 
integration as follows: 


*P. Debye, “Light near a focal point or line,” Ann. d. 
Physik 30, 755 (1909). The calculation given applies for 
an unpolarized wave. 
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sin? (a/2) 
I(r) =4ak f Jo2(2kr(y—y*) dy 


0 


1 x sin*a 
=2k f J?(krx sina)- dx 
0 


(1—x? sin*a)! 





2a : Jo?(u)-u-du 
br? 0 (1 —(u2/k*r2)) 


The last expression is readily integrated if the 
quantity u*/(kr)* is neglected as compared with 
unity. This is quite generally permissible, since, 
even if the pattern is referred to the object 
plane, u?/(kr)? has a maximum value less than 
0.0001. We then find: 


I(r) =k sin’a{ Jo*(kr sina) +J,*(kr sina)}, (6) 


where J;(x) is the Bessel function of the first 
order. For very large r this approaches the 
distribution to be expected for chromatic aberra- 
tion alone : 


2 sina 
I(”)———— 
r 
Jo(x)—(9#x/2)-! cos(x — 4/4); 


J \(x)—(9x/2)-! sin(x — 2/4). 


since 


Close to the Gaussian image point, on the other 
hand, the distribution is closely approximated by 
the distribution for a monokinetic imaging pencil 
in the plane of sharp focus: 


Io(r) = (sin’a/r?) J7(kr sina), (7) 


since for small values of the argument x the 
functions Jo?(x)+J,7(x) and 4J/,;?(x)/x? are both 
given by 1—(x?/4)+ terms of the fourth and 
higher orders in x. 

Figure 5 compares the intensity distribution 
calculated for uniform energy distribution (Eq. 
_(6)) with that obtained for monokinetic electrons 
(Eq. (7)). Inspection indicates that the resolving 
power obtainable with the two distributions is 
practically the same. This is further confirmed by 


5 In this case 








8 sin 
J =4x sin??-2 f dz 2 tana_ 2a 
=47 s =e = ° = 
2 i/tana rz? tanta = tan®a@ r r 
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the distributions obtained for two scattering 
centers separated by small distances, plotted in 
Fig. 6. If the least resolved distance is defined 
as that for which a minimum appears between 
the two maxima corresponding to the individual 
scattering centers, a slightly smaller separation 
may actually be resolved with the uniform energy 
distribution than with the monokinetic radiation. 
In any case, the energy distribution of the 
electrons reduces the contrast. This contrast 
reduction is, however, not so great as to impair 
materially the discrimination of particles in dark- 
field images. It would be more serious in bright- 
field observations. ; 

The conclusions regarding the effect of chro- 
matic aberration on resolving power which have 
here been reached apply throughout the portion 
of the image field in which axial chromatic 
aberration exceeds chromatic change in magnifi- 
cation and rotation. Ina magnetic electron micro- 
scope with centered lenses this is commonly true 
for the entire observed field.® 

These calculations indicate that objectives 
with physical apertures somewhat smaller than 
the optimum for the lens should be capable of 
approximately equal resolving powers in both 
bright and dark fields. The reduced resolution 
demonstrated in the published dark-field images 
may reasonably be attributed to the fact that 
sufficiently small apertures were not generally 
used. In this laboratory some dark-field images 
have been obtained with apertures of 4X 10-* rad. 
and 10~* rad. in which the resolution obtained 
was nearly equal to that in the bright-field 
images. However, in neither type of image did 
the resolution approach the theoretical value. 
This last discrepancy has been ascribed to the 
effects of the charging of contamination layers on 
what were originally clean apertures. 

It now appears that lack of resolution in 
dark-field images has been due to technical 
rather than fundamental limitations. This means 
that properly directed development should make 
it possible to exploit the high contrast possi- 
bilities of dark-field images. 


® See, e.g., Zworykin, Morton, Ramberg, Hillier, and 
Vance, Electron Optics and the Electron Microscope (John 
Wiley and Sons, Inc., New York, 1945), p. 646. 
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The Young’s Modulus and Strain Coefficient of Resistivity of 
Some Bismuth Rich Alloys* 


G. C. KuczyNskI 
Sylvania Electric Products, Inc., Bayside, New York 


AND 


J. T. Norton 
Department of Metallurgy, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received February 24, 1948) 


The Young’s modulus and strain coefficient of electrical resistivity were determined for some 
Bi-Pb, Bi-Sn, Bi-Se, and Bi-Ni alloys. It was found that the modulus of elasticity decreases 
rapidly in the regions of solid solution contrary to the predictions of Guillet. The strain coeffi- 
cient of resistivity was found to increase rapidly within the limits of solid solubility while the 
stress coefficient of electrical resistivity remained constant. 





HE alloys of bismuth with lead and tin 

exhibit unusually rapid increases of specific 
electrical resistivities with small concentration 
of these elements.':? 

On the other hand, the work of Nishayama’® 
and Guillet* seems to indicate that in the range 
of terminal solid solutions, in general, no un- 
usual variations of Young’s modulus should be 
expected. Therefore, it seemed of interest to 
investigate the Young’s modulus and the in- 
fluence of strain upon electrical resistivity of 
these particular alloys. 


METHOD AND MATERIALS USED 


The static method of Young’s modulus deter- 
mination was used because it facilitated the 
measurement of elastic and electrical properties 
at the same time. The apparatus designed for 
this purpose in the Metallurgical Department of 
Massachusetts Institute of Technology is de- 
scribed in another place.’ The increments of 
stress were plotted against the corresponding 
increments of strains. The slope of the linear 
part of the curve thus obtained was accepted 
as the value of Young’s modulus of a given metal 
or alloy. The temperature was maintained at 





* This paper is based — the ScD thesis of G. C. 
Kuczynski, Massachusetts Institute of Technology, 1946. 

!N. Thompson, Proc. Roy. Soc. A155, 111 (1936). 

2 A. Thomas and G. Evans, Phil. Mag. 16, 329 (1903). 

*Z. Nishayama, Sci. Rep. Tohoku Imp. Univ. 18, 359- 
400 (1930). 

4L. Guillet, Bull. Tech du Bur, Veritas. 

5J. T. Norton and G. C. Kuczynski, Instruments, 21 
328 (1948). 
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28°C and controlled within the range +0.1°C. 
The total systematic error in Young’s modulus 
determination was estimated not to exceed 5 
percent. 

At the same time, using the same apparatus 
the electrical resistance at increasing strain was 
measured and its change AR divided by original 
resistance. Ro was plotted against increment of 
strain Ae. The slope of the obtained curve at 
Ae=0 (extrapolated) is defined as strain coeff- 
cient of electrical resistance AR/RoAe. The rela- 
tionship between this coefficient and the corre- 
sponding coefficient Ap/pAe of electrical resistivity 
p in the elastic region can be expressed by the 
following equation: 














Ap/pde= (AR/RAe) —(1+2y), (1) 
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Fic. 1. Young’s modulus as a function of 
composition for Bi-Pb alloys. 
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Fic. 2. Young’s modulus as a function of 
composition for Bi-Sn alloys. 
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Fic. 3. Young’s modulus as a function of 
composition for Bi-Se alloys. 


where yu is the Poisson’s ratio which, in the case 
of bismuth and its diluted alloys, is equal to 
0.33. The changes in length were measured with 
a cathetometer, down to 0.0005 cm. The perma- 
nent elongation never exceeded 0.01 percent. 
The total systematic error in the determination 
of this factor was estimated not to exceed 8 
percent. 

Some Bi-Pb, Bi-Sn, Bi-Se, and Bi-Ni were 
investigated. Bismuth metal was supplied by the 
American Smelting and Refining Co. and con- 
tained no more than 0.0023 percent Pb, 0.0013 
percent Cu, and 0.0005 percent Fe. Lead used 
was chemically pure containing more than 99.9 
percent of lead. Tin was of the Bureau of Stand- 
ard purity. Selenium containing more than 99.9 

percent of selenium used for preparation of bis- 
 muth-selenium alloys was supplied by the Ameri- 
can Smelting and Refining Company. Nickel used 
for making bismuth nickel alloy was Mond Nickel 
containing more than 99.9 percent of Ni. 

The alloys were prepared by melting in vacuum 
then extruded to 0.004-inch diameter wires. 
These wires were examined microscopically and 


6 Smithsonian Physical Tables (1934), Eighth edition, 
p. 137. 
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Fic. 4. Strain coefficient of resistivity for Bi-Pb alloys. 


only those with smooth surfaces were selected 
for investigation. From these only pieces which 
exhibited uniform cross section to 3 percent were 
tested. All wires were vacuum annealed at 100°C 
for 15 minutes. The modulus of elasticity and 
strain coefficient of resistivity for each metal or 
alloy was determined at least from five speci- 
mens. The current passed through the wire 
never exceeded few milliamperes. The electrical 
resistivities of the investigated wires indicated 
that the axis of the wire was parallel to the 
trigonal axis of the bismuth unit cell. All factors 
measured refer to this direction only. 


EXPERIMENTAL RESULTS 


The following Bi-Pb alloys were investigated: 
0.1, 0.2, 0.4, 0.5, 1.0, and 2.5 atomic percent of 











TABLE I. 
Young's 
Concentration Modulus X10" Ap/pAe X10-" 
in atom percent (dyne /cm?) Ap /pAe (cm?/dyne) 

Bi-Pb alloys 
0.00 2.36 —11.11 —4.10 
0.10 1.46 — 6.00 —4.10 
0.20 1.43 — 6.30 —4.30 
0.40 1.41 — 6.10 —4.25 
0.50 1.66 — 7.00 — 4.30 
1.00 1.64 — 8.66 —5.30 
2.50 1.53 — 9,98 — 6.50 

Bi-Sn alloys 
0.17 2.25 — 9.60 —4.2 
0.87 1.60 — 0.32 
1.77 1.62 — 5.82 
8.45 1.70 — 2.32 

Bi-Se alloys 
0.27 1.81 — 9.16 — 5.05 
0.54 1.76 — 8.68 — 5.00 
1.33 1.70 — 8.06 — 5.05 
2.63 1.70 — 7.68 — 4.80 

Bi-Ni alloy 

— 2.08 
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lead. The results are represented in Fig. 1 and 
in the Table I. The smallest addition (about 0.1 
atomic percent) of lead lowers Young’s modulus 
by about 40 percent. After the minimum point 
at about 0.3 atom percent of lead is reached, 
Young’s modulus increases rapidly. After 0.5 
atom percent of lead the alloys exhibit only slight 
changes in Young’s modulus. The rapid decrease 
of modulus of elasticity followed by its rapid 
increase takes place in the region of solid solu- 
bility of lead in bismuth. The solubility limit 
of lead in bismuth was estimated by Guertler 
and other investigators’*® at 0.5 atom percent 
of lead. 

The bismuth-tin alloys behave initially similar 
to bismuth-lead alloys (Fig. 2, Table I). They 
exhibit a shallow minimum at about 1.3 atom 
percent of tin. According to M. Hansen*® the 
solubility limit of tin in bismuth is about 1.5 
atom percent. 

B' smuth-selenium alloys (Fig. 3, Table 1) 
show the same rapid decreases of Young’s 
modulus at low concentrations of selenium, but 
the modulus of elasticity of the alloys containing 
more than 0.27 atom percent of selenium falls 
down slowly and finally assumes a constant 
value. The initial decrease due to the addition of 
one atom percent of selenium is about 30 per- 
cent. The solubility of selenium in bismuth is 
not known but most probably is very low, of the 
order of a tenth of atom percent. 


Ap 


Only one alloy of bismuth-nickel containing 
13 atom percent of nickel was investigated. Its 
Young’s modulus was found to be 2.2510" 
dynes/cm? almost equal to that of pure bismuth 
(2.3 X10" dynes/cm?). 

The above observations seem to contradict 
Guillet’s contention‘ that in the range of terminal 
solid solutions, Young’s modulus does not change 
appreciably. 

The strain coefficients of electrical resistivity 
are given in the third column of Table I and 
Figs. 4-6. The comparison of these curves with 
those of Young’s modulus for the correspond- 
ing alloys reveal striking relationship. Abrupt 
changes of both properties occur within the same 
regions of concentrations. In all cases Young’s 
modulus decreases rapidly as strain coefficient 
of electrical resistivity increases (algebraically). 
In the case of Bi-Pb alloys Young’s modulus 
reaches its minimum at 3 atom percent of lead 
and at the same concentration Ap/pAe has its 
maximum. In the case of Bi-Sn alloys the in- 
crease of strain coefficient of resistivity is so 
rapid that at about 0.8 atom percent of tin it 
reaches value zero and then becomes positive 
and as high as 6.0 at 1.75 atom percent. After 
solubility limit of tin in bismuth (about 1.5 atom 
percent) is exceeded, Ap/pAe falls slowly to the 
value of 2 for 8.45 atom percent of tin. 

In the case of Bi-Se alloys the change of strain 
coefficient does not seem to be so abrupt as in 
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Fic. 5. Strain coefficient of re- 
sistivity for Bi-Sn alloys. 
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* M. Hansen, Aufbau der Zweistofflegierungen (Verlag Julius Springer, Berlin, 1936). 


VOLUME 19, JULY, 1948 


685 














SMTi» Sereda 


Ap 
pie 
—6 
af } 


-g - ———— ied 


—Q : 
—10 
—11 


—12 











j | | | | 
0.5 1.0 1.5 2.0 2.5 
Atom percent Se 





Fic. 6. Strain coefficient of resistivity for Bi-Se alloys. 


the previous cases. This is most probably due to 
the fact that the solubility of selenium in bismuth 
is well below 0.1 atom percent of selenium, so 
the line joining the value of Ap/pdAe for pure 
bismuth and that for 0.27 atom percent of 
selenium should be much steeper (like the dotted 
lines on Figs. 3 and 6). The same refers to 
Young’s modulus curve described above. 

The strain coefficient of resistivity for Bi-Ni 
alloy with 13.0 atom percent nickel was also 
determined and found to be — 4.68. 

The above results show that there seems to 
be a definite relationship between Ap/pAe and 
Young’s modulus within solubility region. The 
shape of the curves suggests that the ratio 
(Ap/pAe)/E where E is Young’s modulus should 
remain constant. In view of EAe= Ao where Ac is 
the increase of stress, we can express our rela- 
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tionship in the following equation: 
Ap/pAo =constant. 


In other words the relative change of the elec- 
trical resistivity of a terminal solid solution of 
Bi-Pb and Bi-Sn alloys is proportional only to 
the stress and independent of either its elastic 
deformation produced or the initial resistivity or 
the composition. Therefore, it is the same as in 
the pure bismuth, namely, 4.10. 

Indeed, the fourth column of Table I shows 
that this factor remains reasonably constant for 
Bi-Pb alloys for concentrations lower than solu- 
bility limit (about 0.5 atom percent of Pb). 
After this limit is exceeded Ap/pAo decreases. 
Of course no such analysis is possible for Bi-Sn 
alloys because some have zero or even positive 
strain coefficients of resistivity. Bi-Se alloys 
exhibit low value of the stress coefficient of 
resistivity (around —5.0X10-" (cm?/dyne) in- 
stead of —4.1X10-" (cm?/dyne)), indicating 
once more that the solubility limit of selenium 
in bismuth is well below 0.27 atom percent of 
selenium, Ap/pAo for Bi-Ni alloy is unusually 
high, — 2.08 X10-" (cm?/dyne). 
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Electron Optical “Schlieren” Effect 


L. MARTON 
National Bureau of Standards, Washington, D. C. 
April 19, 1948 


EVERAL months ago the problem came up of meas- 
uring the magnetization produced on ferromagnetic 
wires in a wire-recording system. Several possibilities were 
explored and finally an electron optical method was chosen 
as one of the possible methods for measuring the fringe field 
extending from the recording wire. This method was con- 
ceived in close analogy to the light optical “schlieren” 
method (often also called Tépler’s shadow method). It 
consists in forming an image of a light source on a con- 
venient stop which intercepts all direct rays (Fig. 1). If in 





SOURCE 


OBJECT STOP IMAGE 


Fic. 1. Schematic representation of the optical path for the 
production of a “‘schlieren”’ image. 


the space between the source and lens there is a variation 
of the index of refraction, an image of that inhomogeneity 
will be produced by means of the same lens in a conjugate 
image plane beyond the stop. Thus a dark-field image of the 
variation of the refractive medium is formed. 

This method is readily adaptable to electron optical in- 
vestigations. In fact, the isolated observations of striae 
occur often enough in electron microscopical work when 
either the objective lens aperture is not properly centered or 
the action of the beam on some unexpected impurities in the 
system produces a charge on some part of the instrument. 
The patterns thus produced in electron microscopes are 
usually very disturbing and every effort has been made to 
avoid them. In our work here we deliberately produce con- 
ditions under which these patterns are formed and use them 
for the observation of the fields making these patterns. 

The most convenient application to electron optical sys- 
tems is that of a beam-producing system like that of an 
electron microscope with a sufficient number of lenses to 
give the required flexibility in concentrating the trans- 
mitted beam and a convenient stop for the production of 
the dark-field image. The first experiments of this nature 
had for object a determination of the magnetization pro- 
duced by a recording head on a steel or other wire used for 
recording purposes. For producing a “‘schlieren”’ image of a 
magnetized wire, the wire was placed at a distance of about 
5 cm from a magnetic lens with a free bore of about 16 mm 
and irradiated with a 40-ekv electron beam. An approxi- 
mately 6-to-1 image was formed on a fluorescent screen. 
Sharp focusing is facilitated by the fact that the image of 
the magnetized wire appears straight only when in focus, 
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whereas out-of-tocus images show the magnetized regions 
out of line with non-magnetic parts of the wire. Thus the 
out-of-focus image of a straight magnetized wire appears as 
a zigzag line. 

After the image of the wire is properly focused, a center 
stop of about 1-mm diameter is used to intercept the 
direct rays. Figure 2 shows how, by placing the center stop 
partially across the visual field, a partly bright-field and 
partly dark-field image of such magnetized wires can be 
obtained. The ‘‘schlieren”’ image of the magnetic fields ex- 
tending on both sides of the wire is readily recognizable. 

Once images of fields are produced by this method the 
intensity distribution in the image can be used for calcu- 
lating the actual field distribution at the sample. Such 
calculations have been carried out recently in this labora- 
tory by Mr. S. H. Lachenbruch for the case of an ideal 
magnetic dipole. In the case of a magnetized wire which was 
our first example, it is reasonable to assume that the field 
around the wire behaves as if it were produced by an ideal 
dipole. By means of this assumption we can calculate from 
the intensity distribution in the image the actual field 
strength in the immediate vicinity of the wire. 

This electron optical method can be made into a very 
useful tool for the observation of electrostatic or magnetic 
fields where such fields are not susceptible to any other type 
of investigation. As an example, the field distribution in 
electromagnetic wave guides may be cited where the field 
distribution in elbows is rather incompletely known. Ex- 
periments have been undertaken here by Mr. D. Reverdin 
to produce such “‘schlieren”’ images by using standing waves 
in a straight portion of a wave guide. It is also an interesting 
tool for the exploration of space charge fields, or any other 
phenomena which do not always allow an analytical treat- 
ment because of their complicated geometries. 

Further applications include a study of contact po- 
tentials, mosaic structure of crystals, ferromagnetic do- 
mains, etc. A preliminary study of ferromagnetic domains 





Fic. 2. ‘“‘Schlieren"’ image of short magnetic pulses on steel wire. 
Left side: dark field; right side: bright field. 
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has been undertaken and is being published elsewhere.! It 
may be worth mentioning in this respect also an early work 
of Boersch, who observed dark-field images produced by a 
vapor jet.? Strictly speaking, his observation does not 
belong here because the primary phenomenon is scattering 
instead of deflection produced by a variation of the electro- 
magnetic index of refraction. 


1L. Marton, Phys. Rev. 73, (1948). 
2H. Boersch, Zeits. f. Physik. 107, 493 (1937). 





Here and There 


Personnel 








J. Howard Dellinger, Chief of the Central Radio 
Propagation Laboratory of the National Bureau of Stand- 
ards, retired on April 30 after 40 years of government 
service. He is now acting as a radio consultant and adviser. 

Isidor Fankuchen is the first incumbent of the new Chair 
of Applied Physics at Polytechnic Institute of Brooklyn. 

Christopher E. Barthel, Jr., has been named Chairman 
of Physics Research at Armour Research Foundation of 
Illinois Institute of Technology. 

Carroll Stansbury has been appointed to the Electronics 
Division of the National Bureau of Standards, where he 
will do research in the engineering electronics laboratory. 

T. Keith Glennan was inaugurated as the fourth Presi- 
dent of Case Institute of Technology, Cleveland, Ohio, on 
May 20 and 21. The program featured a symposium with 
the theme: The Increasing Responsibility of the Scientist 
and Engineer in Public Affairs. 


Lehigh Valley Physics Club 


The Physics Club of the Lehigh Valley was recently 
elected to the status of an Affiliated Society of the American 
Institute of Physics. The Club at present has a membership 
of about one hundred scientists and engineers. 


Eisenhart Elected Fellow of British Society 


Churchill Eisenhart, Chief of the Statistical Engineering 
Laboratory, National Bureau of Standards, has been 
elected a fellow of the Royal Statistical Society of Great 
Britain 


North American Philips Sponsors Course 


A group of nationally known authorities on x-ray 
diffraction lectured on their specialties before the third 
annual spring course in x-ray diffraction, spectrometry, and 
fluorescence analysis, held recently in the New York offices 
of North American Philips Company, Inc. Discussion of 
medical applications of x-ray diffraction was included in the 
course for the first time. Seventy persons, from 15 States 
and Canada, enrolled in the week-long course. 





New Booklets 








The Proceedings of the 1947 National Electronics Con- 
ference has been printed and is now available at $4 a copy 
from Dr. R. R. Buss, Secretary N.E.C., Care of Electrical 
Engineering Department, Northwestern Technological In- 
stitute, Evanston, Illinois. 

Lenkurt Electric Company, Inc., 1105 County Road, San 
Carlos, California, has issued its Bulletin 40A, entitled 
Trancors. 24 pages. 

Industrial Photographic Division of Eastman Kodak 
Company, Rochester 4, New York, recently issued two new 
booklets which are available free on request: Magnifying 
Time (12 pages) and Kodak Linagraph Films and Papers, 
for Instrument Trace Recording (4 pages). 

RCA Laboratories Division, Princeton, New Jersey, 
published this spring an index to RCA Technical Papers, 
Volume IT (b), 1947. 24 pages. 
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